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ABSTRAOT 


The effect of changes in strain rate and temperature on the flow stress 
of Mg single crystals was determined in the range of 4:2° to 420°x. The 
present data, along with those obtained previously, indicate that the 
intersection of dislocations is the controlling mechanism.  Cottrell’s 
original model must, however, be modified to include an effect of stress 
on the activation volume. This effect appears to be due primarily to the 
fact that the force between two intersecting dislocations varies with their 
distance of separation rather than a change in the effective forest spacing. 


§ 1. INTRODUCTION 


Tue deformation of metal crystals is in many cases governed by thermally- 
activated dislocation mechanisms. If a single activation process controls 
the deformation, the strain rate y is given by 

y=vexp(—G/kT) Bed Phan SU agus es at) 
where v is a frequency factor which includes essentially the frequency of 
vibration, the number of places where the activation can take place, and 
the strain per successful fluctuation. G=HW—T'S is the Gibbs free energy 
of activation where 7 is the enthalpy of activation (commonly termed the 
‘activation energy’) and S is the entropy of activation. Conrad and 
Wiedersich (1960) have shown that the value of H can be obtained from 
flow stress data by the following equation : 


ae) _ 5 de r (et) 
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where 7 is the applied shear stress and 7, is the long-range internal stress 
which is proportional to the shear modulus p. +*=7—7, is the excess of 
the applied stress over the internal stress. It is called the effective stress 
or thermal component of the flow stress. 

Two mechanisms have been proposed which control the movement of 


dislocations in their glide plane in metal crystals with a close-packed 
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structure (c.p.h. and f.c.c.): (a) the mutual intersection of dislocations 
(Cottrell 1952, Nabarro 1952, Seeger 1954) and (6) the overcoming of the 
Peierls stress (Conrad 1958, Lytton et al. 1958, Weertman 1958). One 
major difference between these two mechanisms is the effect of stress on 
the activation energy H. For the intersection mechanism Cottrell (1952) 
and Seeger (1954) have suggested that the activation energy is decreased. 
linearly by the effective stress, 1.e. 


H=H,—1ge 2 ee ee 


where v, is termed the ‘activation volume’. v=d lb where d, is the inter- 
secting width (approximately twice the width of the gliding dislocation) ; 
1, is the average spacing between the dislocations threading the glide 
plane (termed the ‘forest’); and b is the Burgers vector. Modifications 
of eqn. (3) have been proposed which consider that the activation volume is 
itself a function of stress (Friedel 1956, Basinski 1959) and temperature 
(Thornton and Hirsch 1958). For the Peierls mechanism Seeger (1956) has 
calculated the activation energy to be 


H=Hy (1+ 4log~ ) Ree ee) ot: 
TT 

for the limit r*<t,°. H, is the energy ofa kink and 7,,° is the value of the 

Peierls stress at 0°K. 

For the deformation of f.c.c. metal crystals in Stages I] and IIT (multiple 
glide) there is now considerable evidence to support the intersection of 
dislocations as the rate-controlling mechanism at low temperatures 
(Friedel 1956, Seeger 1957, Thornton and Hirsch 1958, Basinski 1959). 
For Stage I deformation of f.c.c. crystals or the deformation of ¢.p.h. 
crystals (single glide or so-called easy -glide), the determination of the rate- 
controllmg mechanism has been somewhat more elusive. The major 
problem has been that the experimentally determined flow parameter 
B=(d0lny/dr), is independent of temperature (Conrad and Robertson 
1957, 1958, Conrad 1958) rather than proportional to the reciprocal of the 
temperature as required for the intersection mechanism when eqn. (3) 
describes the activation energy. More recent work on Mg crystals by 
Conrad et al. (1960) has, however, indicated that this apparent discrepancy 
is probably a reflection of the fact that the activation volume is itself a 
function of the stress and that in order to obtain a constant strain rate the 
stress must increase with decrease in temperature. 

Since this latter investigation (Conrad et al. 1960) was limited in the 
temperature range covered (78° to 364°K) and in the number of temperatures 
considered (four), it appeared that additional work over a larger temperature 
range should be conducted to positively identify the controlling mechanism 
during single-glide (easy-glide). Consequently, in the present investiga - 
tion the deformation behaviour of Mg single crystals was investigated in 
considerable detail over the range of 4:2° to 420°K. Of special interest 
was the determination of the effect of temperature on the flow parameters 
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7 and B. From such measurements the effect of temperature and stress 
on the activation energy H and the activation volume v was determined. 


§ 2. EXPERIMENTAL PROCEDURE 

The Mg crystals employed in the present study were approximately 1 cm 
diameter and 15cm long. They were grown from the melt in a graphite 
mould under a static atmosphere of He from redistilled Mg of <0-02% 
impurities. Their orientations were determined by the Laue back-reflec- 
tion technique and are given in the table. Previous to testing, they were 
immersed in 25° by volume HCl to remove a surface layer about 0-1 cm 
thick, and then chemically polished in 20° by volume nitric acid. Fol- 
lowing polishing, the specimens were rinsed successively in tap water, 
distilled water, and acetone, and then allowed to dry in air. 


Initial orientation and critical resolved shear stress 7, of Mg 


crystals 
, °. y 3 
Crystal Xo (1) Ay (2) TK ae ‘ 102 g/mm? 
2 28 40 Yel 1-35 39 
6 23 34 78 1-44 109 
12 43 ol 85 1-15 68 
16 29 38 300 1-70 47 
2] 18 32 4:2 2-50 1 


(1) x)» =Angle between specimen axis and basal plane. 
(2) A,=Angle between specimen axis and (100). 
(3) y=Resolved shear strain rate. 


The experimental parameters of eqn. (2) must be obtained for a constant 
structure. This was accomplished by employing so-called differential 
tests, i.e. tests where the strain-rate or temperature are changed during the 
deformation. Two testing machines were employed: an Instron Model TT-B. 
and a special screw machine built by Areweld Mfg. Co.}. In general, the 
Instron machine was used for tests above 78°K while the Arcweld machine 
was employed for tests between 78° and 4:2°K. Several tests were also 
conducted on the latter machine in the range above 78°K to obtain an oueston 
of data. Only the Arcweld machine was employed for tests below 78°K 
because of the larger space between cross-heads, which permitted the 
i ion of the helium cryostat. 

Barth Ib SR-4 load call connected to the X-axis of an Instron X-Y 
recorder was employed to measure the load in both machines. This gave a, 
stress accuracy of approximately 2%. The cross-head movement was 
measured by a clip gauge employing SR—4 resistance wire gauges connected. 


+ Arcweld Manufacturing Co., Grove City, Pa. 
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to the Y-axis of the recorder, giving a resolved shear strain sensitivity of 
approximately 5x 10->. The strain rate in both machines could be rapidly 
changed by factors up to 1 : 2000 over the range of resolved shear strain rate 
of 10-*min- to 10min“. 

Temperatures between room temperature and 78°k were obtained by 
the method described previously (Conrad 1958), which employs a nitrogen 
vapour spray on the specimen as the refrigerant. To obtain temperatures 
above room temperature, resistance heating elements were installed in the 
low-temperature container. In all cases the temperature was measured 
and controlled by a previously calibrated copper—constantan couple 
attached to the specimen and connected to a L and N Speedomax Type G 
recorder-controller with a D.A.T. anticipatory control. The temperature 
was controlled within +1°K over the length of the specimen from 78° to 
420°K. 

Fig. 1 


[I olet for Lrquid Me 


Stayonory Pull Rod 


OS foo Moving Crossheod 


/ 


~;——— Guter Dewar 


Inner Dewor ———-}— — 
Specimen 


Perforated Tube 


Cold Chomber 


Liquid No 


Schematic of helium cryostat. 


Temperatures between 78° and 4:2°K were obtained by the cryostat 
shown in fig. 1. A temperature of 4:2°K was obtained by filling the inner 
chamber with liquid He, while temperatures between 78° and 4:2°K were 
obtained by controlling the flow of He into the cryostat from a large He 
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dewar in the manner described by Wessel (1956). By this method the 
temperature was maintained constant to within + 1°K. 

To make a temperature change, the specimen was always unloaded to 
about 50% of the flow stress, to prevent creep from occurring during the 
change. A decrease in temperature generally required about 2min, 
while an increase required about 5-10 min. In all cases, the test was con- 
tinued only after the temperature, indicated by the couple attached to the 
specimen, had remained constant for at least 5min. An indication that 
thermal equilibrium had been established was also provided by the fact 
that no further change in dimensions associated with thermal contractions 
or expansions was observed. 

The stresses and strains reported in the present paper are all resolved 
shear stresses and resolved shear strains. 


§ 3. EXPERIMENTAL RESULTS 


Figures 2 and 3 show typical stress-strain curves with temperature and 
strain-rate changes. No marked transients such as those reported by 
Conrad (1958) and Basinski (1959) for f.c.c. metals were observed following 
the changes in temperature and strain rate. 


Fig. 2 
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Effect of temperature and strain rate on the flow stress. Also shown are 
values of the parameter B=(d In y/07)r. 


Also given in figs. 2 and 3 are values of the parameter B=(dln [0r)p 
determined from the changes in strain rate. For temperatures of 40°K and 
above, B was independent of strain and relatively independent of Ronee 
ture, in agreement with earlier results (Conrad and Robertson 1957, 1958, 
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Conrad et al. 1960). However, at 20°K and below, the value of B was larger 
than that for the higher temperatures. Also, in this temperature range B 
decreased with strain. Furthermore, as seen in fig. 3, straining at these 
lower temperatures caused a decrease in B determined subsequently at 
higher temperatures. Thus, at temperatures below 40°K, the effect of 
strain on B for Mg crystals is similar to that observed for Stages I and ITT 
inf.c.c. metals. This feature is being investigated in more detail and will be 
covered in another paper. In the present paper only those values of B will 
be considered which do not include an effect of straining. 


Fig. 3 
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Effect of temperature and strain rate on the flow stress. Also shown are 
values of the parameter B=(0 In y/dr)7. 


Figure 4 shows the variation of 6 with temperature. The increase in B 
below 40°K is clearly evident; also there appears to be a minimum in B 
at ~60°K. The present values of B between 78°K and 300°K lie between 
those obtained previously (Conrad and Robertson 1957, Sheeley e¢ al. 1957, 
Conrad et al. 1960) for Mg crystals at approximately the same strain rate. 
Within the scatter of the data, B tended to decrease with increase in 
effective strain rate} y, over the range of 10-$ min“ to 10 min“ in agree- 
ment with the earlier results (Conrad et al. 1960). However, this feature 
was not studied in detail in the present investigation. 

According to Seeger (1954), the flow stress 7 can be separated into two 
additive components, a thermal component +* which is identical to the 
ee Nae mentioned effective stress and an athermal component Tp80 

at 

dak ile ee em a IER th 
7, is the long-range internal stress and consequently only indirectly 
dependent on the temperature through the shear modulus j.; it is indepen- 
dent of the strain rate. 7*, however, is a function of strain rate and tem- 
perature. 


+ Ye= (12)? where 7, is the initial strain rate and 7, is the final strain rate. 
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To evaluate the activation energy we have to derive r* as a function 
of temperature. One procedure is to determine the critical resolved 
shear stress rt) as a function of temperature with y=const. 7, is then 
obtained from 7) at high temperatures where +*=0. However, this 
requires testing a number of crystals, each at a different temperature, and 
the results are always subject to specimen variations. A better way to 
obtain 7* is from the change in flow stress which occurs with a change in 
temperature. This method seemed especially favourable in the present 
investigation since the previous work on Mg crystals had indicated that 
only z, is increased by work hardening, whereas r* is independent of strain. 
Therefore the separation of the measured flow stress into r* and 7, 8 possible 
without any doubtful assumptions. 


Fig. 4 
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Variation of B=(9 In y/d7)7 with temperature for Mg crystals. 


7, can be calculated at any point on a stress-strain curve such as is shown 
in figs. 2 and 3 by adding the stress increases due to the work hardening 
Az, (after correction for the change in shear modulus), to the initial internal 

7 


stress Ges 


n At, 
T,)=tm, (2+ 3—f sO 
ih 7 br, Lo mie 
whereby the summation is carried out over the pieces of the stress-strain 
curve obtained at the temperatures 7';; 7* is then obtained from eqns. (5) 


and (6): | roe 
(Di) =a(D,9)— wl P) | + SH. i 2 hue) 


In the present investigation 7, was determined directly for only one 
erystal (Mg-2) by obtaining the critical resolved shear stress at 377 
which is greater than the temperature at which 7* becomes zero (~ 330°K 
for y=10-? min“). 7, for the other crystals was then chosen so that their 
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7* at one temperature (generally ~80°K) coincided with 7* of ee 
The values of 7, for these crystals then ranged from 37 to 47 g/mm*, 
which is the order of spread one generally obtains in the critical resolved 
shear stress for different crystals at temperatures above 330°K. 

The values of +* (for y ~ 10-? min) obtained by eqn. (7) are given as a 
function of temperature in fig. 5+. There is good agreement between the 
various crystals and between the values obtained in this manner and 


Fig. 5 
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(0) 100 200 300 400 
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Effect of temperature on the thermal component of the flow stress 7*. 


those derived from the critical resolved shear stress values designatedz,. Of 
significance with regard to the curve of fig. 5 is that an inflexion occurs 
at approximately 30°K, i.e. at the temperature below which the parameter 
B increases strongly. Also one should note that 7*=0 at 7’ >330°K. 


§ 4. Discussion 


Employing eqn. (2) one obtains from the data given in figs. 4 and 5 the 
variation of H with temperature shown in fig. 6; i.e. the activation energy 
for plastic flow of Mg crystals is proportional to the temperature. Taking 
the logarithm of eqn. (1) and rearranging one obtains: 


Hee kT n'y }pa es 


where v'=vexp(S/k). The proportionality between H and temperature in 
fig. 6 thus indicates that v’ is a constant over the stress and temperature 
range considered. Hence, either v’ is independent of temperature and 
stress or the combined effect of temperature and stress on v’ are such to 
yield a constant value of v’ at a constant strainrate. To a first approxima- 
tion it was assumed in the present work that v’ is independent of stress and 
temperature. This approximation is reasonable if the effect of stress and 


SE EE nee Sena EE See e NE Mr EA ey SEE bn, 
+ The shear modulus values (C,,) were taken from the paper by Slutksy and 
Garland (1957). 


Thermally-activated Glide in Magnesium Crystals 185 


temperature on the strain rate through v’ is small compared with their effect 
through the Boltzmann exponential. The agreement between theory and 
experiment to be discussed below lends support to this assumption. 

Since the previous work (Conrad e¢ al. 1960) indicated that H for the 
plastic flow of Mg crystals is primarily a function of stress alone (i.e. £f is 
relatively independent of temperature per se), one can obtain the effect of 
stress on H by substituting the stress from fig. 5 for the temperature in fig. 6. 


Fig. 6 
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Variation of the activation energy H for plastic flow of Mg crystals 
with temperature. 


Fig. 7 
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Effect of stress on the activation energy H for plastic flow of Mg crystals. 


This has been done in fig. 7. ‘The dependence of H on stress observed here 
supports the mutual intersection of dislocations as the rate-controlling 
mechanism over the Peierls mechanism. Furthermore, the value Ot. 
(16-3 keal/g atom) is in reasonable agreement with what one estimates for 
the energy to form a jog in Mg (14-16 kcal/g atom). For further comparison, 
it is desirable to determine the values of the frequency factor v, the forest 
spacing 1, the intersection width d and the force-distance relationship for 
the intersection mechanism from the experimental data. 
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It can be shown (Conrad et al. 1960) that the physically-significant 
activation volume v* for the intersection of dislocations is given byt 


vet =ldb=BkT ee eee te 


where J and d may depend on the stress 7*. An effect of stress on d will 
occur if the force between two intersecting dislocations varies with their 
relative distance of separation; see for example Basinski (1959). That l 
may decrease with stress due to the bowing-out of a dislocation in the slip 
plane has been suggested by Friedel (1960). The variation of v* with stress 
obtained from eqn. (9) by using the values of B given in fig. 4 and the 
values of 7 as a function of r* from fig. 5 is given in fig. 8. A strong effect 


Fig. 8 


t* g/mm* 


Effect of stress on the activation volume v* for plastic flow of Mg crystals. 


of 7* on v* is clearly evident. From vg=v,* (r*=0) one can obtain an 
estimate of the average forest spacing J): 


b= & = =2x104em |)... BP ee 


which is in agreement with measurements of densities of dislocations in 
metal crystals (Cottrell 1953). 

Knowing the forest spacing one can estimate the value of v’ and thus 
compare the estimated value with that given by the slope of fig. 6. Accord- 
ing to Seeger (1954) v’ is given by 


y = NA by exp(S[kh oe. ee sees ay 


where JN is the number of dislocation segments per unit volume contributing 
to the plastic flow, A is the area swept-out per successful fluctuation, b 
is the Burgers vector and v, is the frequency of vibration of a dislocation 
segment. Taking V=1,-*, A =1,? and v9=vab/2lo, where va is the Debye 


nn ee eae rien RNs Sees nine Ne te Pe Ls 

} v*7* is the work done by the applied stress during the thermal activation. 
In general v* is not equal to vp; v* =v9 only for 7*=0 or if the force—distance 
relationship is a square-type function. 
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frequency (10% sec"), and neglecting the entropy contribution one obtains 
v'=1x 10sec. This is somewhat lower than the value obtained directly 
from the slope of fig. 6 (107 sec“'); however, the agreement can be con- 
sidered satisfactory in view of the assumptions made. 


Fig. 9 
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| Some indication of the form of the force-distance relationship for inter- 
section is provided by the plot of +* (= /’/bl) vs. v* (= bid) in fig. 8. How- 
ever, such a plot only gives the true force-distance curve (aside from the 
Sector bl) when two conditions are fulfilled : 
(a) The effective forest spacing / does not vary with stress. 
(b) The force-distance curve either decreases to zero force immediately 
following Fmax or is symmetrical. LG pate: 
If the first Sree satisfied, the data for crystals which differ in initial 
forest spacing J, should superimpose when TU, ( = dy) is ploMed versus 
the ratio v*/vy(=d/do). dy represents the ‘activation distance’ at zero 


188 On Thermally-activated Glide in Magnesium Crystals 


effective stress (r*=0). Figure 9 gives such a plot for three separate 
investigations on Mg crystals. It is here seen that the curves superimpose, 
suggesting that / is not significantly affected by the stress. 

In regard to the actual shape of the force-distance curve, one cannot 
make a decision based on the experimental data. Also, the theory is not 
developed in sufficient detail to be of help. Since there is no reason to 
expect the force-distance curves for intersection to have the special forms 
listed, one must assume that the curves of fig. 9 only represents the force- 
activation distance relationship. Similarly, the ‘force-distance’ curves 
derived by Basinski (1959) for the deformation of Ag, Al and Cu single 
crystals in Stages Il and III are also actually force-activation distance 
curves. If one subtracts the long range stress component from Basinski’s 
curves, their form is similar to that given in fig. 9. 


§ 5, CONCLUSIONS 


In conclusion the present and previous experimental data on the plastic 
flow of Mg crystals indicate that the rate-controlling mechanism at low 
temperatures is the intersection of dislocations. However, Cottrell’s 
(1952) simple model must be modified to include an effect of stress on the | 
activation volume. This effect of stress appears to be due primarily to the 
fact that the force between two intersecting dislocations varies with their 
distance of separation in agreement with Basinksi (1959), rather than a 
change in the effective forest spacing as proposed by Friedel (1956). 
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Radiation-enhanced Helium Precipitation in Copper; 


By A. N. Gotanpt 
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[Received July 16, 1960] 


ABSTRACT 


The effect of neutron irradiation on the precipitation of helium atoms in 
copper has been studied at temperatures up to 425°c. At 160°c and 300°c 
no effect was observed. At 425°c the irradiated sample exhibited considerable 
precipitation around grain boundaries as compared to an unirradiated 
sample held at the same temperature. It is concluded that the additional 
vacancies introduced during the irradiation enhance the precipitation process. 
The resultant voids are equal in size to those obtained at higher temperatures 
in the absence of radiation. 


AccCORDING to the theory of radiation-enhanced diffusion (Dienes and 
Damask 1958) neutron irradiation produces a steady-state vacancy 
concentration in excess of that which would have existed at a given 
temperature in the absence of the radiation field. Consequently, the 
theory implies that processes whose rates are controlled by vacancy 
diffusion should proceed more rapidly during irradiation than under 
normal conditions. This has been clearly demonstrated for the case of 
partially ordered alpha brass (Dienes and Damask 1958). If possible, 
one would like to provide some visual evidence of enhanced diffusion, 
and the phenomenon of helium precipitation in copper (Barnes e¢ al. 
1958) seemed suitable. Barnes and his co-workers have demonstrated 
that helium will precipitate in copper and produce visually observable 
voids when the copper is heated to temperatures where self-diffusion is 
important. These workers believe that the precipitation results from the 
diffusion of vacancies into the vicinity of the helium atoms, and they have 
performed numerous experiments which tend to confirm this hypothesis. 
If one accepts their supposition, then it ought to be possible to induce 
precipitation at a temperature below the self-diffusion range _by 
irradiating copper with neutrons to increase the vacancy concentration, 
provided only that most of the vacancies live long enough to migrate to 


+ Communicated by the Author. Work performed under the auspices of the 


U.S. Atomic Energy Commission. 
t+ Guest Scientist, Ordnance Materials Research Office, Watertown, 
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the nucleation centres. The initial results of such an experiment are 
described in this note. 

Figure 1} shows the successive steps in the experiment. A sample of 
copper 0-040 in. thick with a nominal purity of 99-999% was bombarded 
with 40 mev alpha-particles from the 60 in. cyclotron until a total charge 
of 150 pamp hrs had been injected. This is equivalent to a concentration 
over the bombarded volume of approximately 0-2 atm. % helium. During 
the «-particle bombardment, the sample temperature never exceeded 
180°c, and it was probably considerably lower. After sufficient time had 
passed for the sample to reach a safe level of radioactivity, it was 
unsoldered from the target blank and cut into pieces approximately 
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described in the text. 


tin. wide by lin. high. Autoradiographs of these samples were taken 
to determine which ones had been uniformly bombarded with a-particles. 
This was a necessary task because some regions of the copper target had 
been masked by the grid which held it in the cyclotron beam. To check 
a eee eee 


} Figures 2-4 are shown as plates. 
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the experimental technique, some of the experiments of Barnes et al. 
(1958) were repeated with excellent results. 

After suitable samples had been chosen, one sample was irradiated at 
a selected temperature in the reactor, while another was held at the same 
temperature in a vacuum furnace. Each sample was maintained at 
temperature in its environment for the same period of time, approximately 
10 days. Most of the neutron irradiations were carried out in a furnace 
in hole E-11 of the Brookhaven graphite reactor. The epicadmium flux 
in this experimental facility was approximately 2 x 10" nv as determined 
by activation of bare and cadmium-covered aluminium—manganese foils. 
Following this treatment, the samples were stored until the in-pile sample 
could be safely handled once more. Finally, standard metallographic 
techniques were applied to the edges of the samples in order to expose 
voids which might have been introduced. The samples were mechanically 
polished first and then either electrolytically or chemically etched. 
Whenever voids were exposed by such polishing, they appeared in a 
band approximately 20 to 30 u wide and at a depth below the surface 
of about 0-22.mm. 

The experiments described above were performed at 300°c and 425°c 
in hole E-11, and in a higher flux hole at 160°c. No evidence of void 
formation was found in either the irradiated or contro! samples main- 
tained at 160°c or 300°c. However, the samples held at 425°c both 
revealed that the precipitation of helium and the formation of voids in 
the grain boundaries was proceeding during the experiment. Moreover, 
the sample which was irradiated during the ageing process showed 
considerably more precipitation than the sample which was merely aged 
in vacuo. The surfaces of these samples were polished and etched ten 
or twelve times to ensure the validity of the observations. A total of 
from 5 to 10 mils was removed from the surfaces during this procedure, 
and the original observations were confirmed each time. Figure 2, 
which was taken with the aid of a light microscope, shows this result 
clearly. The surfaces of all other samples including those which were 
untreated except for «-particle bombardment resembled the clear regions 
in the figure. Thus, the results furnish one of the first visual observations 
of enhanced diffusion. 

While providing superficial evidence that an interesting effect exists, 
the light microscope does not supply the detailed image which is really 
desired. For this purpose, one must resort to the electron microscope. 
Consequently, direct carbon replicas shadowed with tungsten dioxide 
were made of the two samples shown in fig. 2 so that the nature of the 
precipitation region could be investigated more fully. Since replication 
techniques generally are peculiar to the laboratory in which they are 
performed, a few words about the techniques employed in this work are 
in order. To begin with, the samples were mechanically polished and 
then etched with a solution consisting of one part H,O., one part H,0, 
and two parts NH,(OH),. This solution must be fresh when it is applied 
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to the sample or oxidation will result. Following a successful etch, the 
replication is carried out under vacuum. A commercial wetting agent, 
Victa-wet, is heated in a tungsten basket at a moderate temperature so 
that any water present is driven off. Then the Victa-wet is heated to 
incandescence for at least one minute, and deposited by evaporation 
onto the surface of the sample. The sample is shadowed, and subsequently 
covered by an evaporated layer of carbon. Finally, the shadowed 
carbon replica is floated off of the copper in distilled water. Replicas 
made in this way resulted in the electron micrographs shown in figs. 3 
and 4. 

Figure 3 (a) is an electron micrograph of the normal grain structure 
of the out-of-pile sample, and fig. 3 (6) is a similar micrograph of the 
in-pile sample, In each case the area shown is immediately adjacent to 
the area in which the helium band is prominent. Each micrograph 
reveals just the normal grain structure of a polycrystalline material, 
the background texture being a result of the etching process. Figures 4 (a), 
(b) and (c) are electron micrographs of the helium bands in the same two 
samples. By comparison with fig. 3 (a), fig. 4 (a) exhibits grain boundaries 
which are more sharply defined due to the precipitation of helium. 
However, this precipitation is on a relatively fine scale. Im marked 
contrast to fig. 3 (b), figs. 4 (b) and (c) reveal the presence of much larger 
voids in the grain boundaries of the reactor-irradiated sample as a 
consequence of helium precipitation during reactor irradiation. Figure 
4 (b) is an example of the phenomenon, already well established by Barnes, 
that certain grain boundaries are either sources or conductors of vacancies 
while others are neither. Following Barnes’ explanation (Barnes 1959) 
of the void formation, we can consider several explanations for the fact 
that the voids are in the grain boundaries rather than separated from 
them by a bubble-free zone of copper. One possibility is that irradiation 
has permitted enhanced localized diffusion of the helium to the grain 
boundaries. However, this seems unlikely because all grain boundaries 
should then be similarly affected. Figure 4(b) contradicts this. 
A second, more plausible, explanation is that the irradiation somehow 
causes some grain boundaries to become good conductors of vacancies 
by increasing their energy. In addition, the vacancy supply does not 
become exhausted readily because it is supplemented by the vacancies 
produced during the neutron bombardment. It is interesting to note 
that the sizes of the voids in figs. 4(b) and (c) are in the range of 
1:5-3x10-°cem. This is the same range as Barnes obtained in copper 
after two hours of heating a sample at 800°c. 

This work will be continued for different fluxes and temperatures in 
order to place the results on a quantitative basis. The theory of 
diffusion-limited precipitation (Ham 1958) predicts the rate of growth 
of voids in terms of the diffusion constant of the material. It is hoped 
that a measure of the enhanced diffusion may be obtained by applying 
this theory to the present technique. 
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ABSTRACT 


Metallographic studies on compressed bi-crystals of magnesium oxide 
have revealed cracks formed at the grain boundary by piled-up groups of 
edge dislocations. However, these were produced only when two slip bands, 
one from each grain, met at the boundary ; cracks associated with a single 
pile-up were not observed. Crack propagation was transgranular, and by 
crystallographic cleavage. 

A mechanism for crack formation at low strains is proposed in which a 
source, situated on the opposite side of the grain boundary and some two 
to eleven microns away from a piled-up group of edge dislocations, can be 
activated by the stress field of this pile-up and produce dislocations which, 
under certain conditions, may cooperate with the activating group to form a 
crack. 


§ 1. LyrRopUCTION 


RECENTLY it has been demonstrated that single crystals of magnesium 
oxide can exhibit a potentially useful measure of ductility (~5%) 
providing that such factors as environment and surface condition are 
carefully controlled (Gorum eft al. 1958, Stokes et al. 1959 a). However, 
in the polycrystalline form, this and other ceramic materials are commonly 
noted for their lack of ductility. This suggests that the presence of grain 
boundaries may be responsible in some way for the brittle nature of 
polycrystalline magnesium oxide. A possible explanation was provided 
by Zener (1948) several years ago: grain boundaries may act as stable 
barriers against which edge dislocations can pile up and coalesce, producing 
a crack nucleus. The theory of this type of crack formation has been 
developed by Stroh (1954, 1955, 1957) and others, and cracks caused by 
dislocation coalescence have been observed at the intersections of bands 
of edge dislocations in single crystals of magnesium oxide (Stokes ef al. 
1959 b) and lithium fluoride (Westwood 1960). Such cracks were 
responsible for the eventual failure of the single crystals. _ However, 
in polycrystalline samples there are several ways in which failure could 
occur: by cracks initiating inside individual grains and then propagating 
transgranularly ; by Zener-type cracks initiating at a grain boundary 
and propagating either intergranularly or transgranularly ; or by failure 
of the grain boundary itself due to impurity embrittlement or the presence 


of voids. 
+ Communicated by the Author. 
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In this work the fracture behaviour of bi-crystals of magnesium oxide 
in compression has been studied. The specific objectives were to elucidate 
the role of the grain boundary in fracture and to determine whether a 
Zener-type fracture mechanism is operative in this material. 


§ 2, EXPERIMENTAL TECHNIQUES 


Polycrystalline magnesium oxide blocks were obtained from the Norton 
Company and bi-crystals, approximately 25 mm x 5 mm x 4 mm in size, 
were cleaved from these blocks so that the boundary was close to the 
centre of the specimen. The change in lattice orientation across 
boundaries was typically about ten to twenty degrees and the rotation 
component was small (+ 5°). 

Specimens were mounted in lucite and parallel compression faces, one 
of which was a {100} cleavage plane, produced by grinding and mechanical 
polishing. After the lucite had been removed the bi-crystals were 
chemically polished for two to four minutes in 85% orthophosphoric 
acid at 100°c, washed in distilled water and carefully dried using ethyl 
alcohol and ether rinses. 

Specimens were then compressed in an Instron testing machine at 
strain rates of order 10-4 per sec, and a birefringence technique was 
utilized to study the strain pattern of the bi-crystal during testing. 
After an appropriate strain, the distribution of dislocations and cracks 
was revealed by the etching technique of Stokes et al. (1958). 


§ 3. EXPERIMENTAL OBSERVATIONS 


Bi-crystal boundaries were found to be extremely effective barriers to 
the movement of dislocations. In particular, the piling-up of groups of 
edge dislocations at the boundary was regularly observed. Often cracks 
were associated with such piled-up groups and these formed in a variety 
of ways, some of which are illustrated in figs. 1-4+. These typical 
micrographs suggest that a crack is formed in the grain boundary only 
when two slip bands, one from each grain, ‘meet’ at the boundary. 
Although a crack is not necessarily formed under such circumstances, no 
cracks associated with a single piled-up group of dislocations were 
observed. In fig. 1 for example, cracks have formed at A and B but 
not at C, D, For F. Similarly, in fig. 2 cracks have not formed at G or H 
where relatively broad bands of edge dislocations are arrested, nor at J 
where two slip bands meet. However, the concept of slip bands meeting 
at the boundary must be more closely defined ; examination of the 
cracks formed at B in fig. 1 and K,, K, and K, in fig. 2 will show that the 
slip bands responsible for these cracks were actually displaced some one 
to two microns from one another at the boundary. This displacement 
is thought to be a significant factor in the fracture mechanism for 
sometimes cracks are not formed when slip bands meet, apparently, 
head-on as in fig. 3 at R, and Rs. 


} Figures 1-6 are shown as plates. 
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Although cracks were initiated frequently at the boundary and only 
rarely in the interior of the invididual grains, these cracks did not. 
propagate along the boundary but always into the grains, usually on 
a {100} cleavage plane, fig. 2, L, or a {110} slip plane, fig. 3, R,. 
Metallographic observations indicate that cracks propagate away from 
the boundary as soon as they reach the edge of the promoting slip band 
as, for example, in fig. 4 at N. 

Slip bands sometimes originate from sources in the boundary as at Q 
in fig. 4. Figure 5 shows that a boundary can be a prolific source of 
dislocations, but this phenomenon was not commonly observed in these 
experiments—dislocations usually originated near the compression faces. 
However, the former may be typical of the mode of deformation of 
polycrystalline specimens. 

The cleavage crack at P in fig. 4 may have resulted from the interaction 
of dislocations in slip bands Q and T. Cracks believed to have been 
formed in this way have been observed previously by Washburn e¢ al. 
(1959). 

§ 4. Discussion 

It has been demonstrated that cracks can be formed at grain boundaries 
by the action of piled-up groups of edge dislocations but that two groups, 
and not one as predicted by Zener (1948). for metals, are necessary for 
this to occur. The meeting of bands of edge dislocations at a boundary 
is inevitable at © high’ strains but there are two ways in which this can 
arise so that cracks may be formed at low strains. First, two sources, 
both far away from the boundary, can emit dislocations which meet 
fortuitously at the boundary. Second, a source close to the boundary 
can be activated by the stress associated with a primary piled-up group 
of dislocations in the neighbouring grain. Stroh (1957) has shown that 
sources most likely to be activated will lie at a distance r=aGb/oy away 
from the primary pile-up, where « is a term of order, but greater than 
unity ; G@ is the rigidity modulus; b the Burgers vector and oy the 
‘upper’ yield stress. For a single crystal of magnesium oxide 

Gia10) =S11 — Siz + 2522 (Voigt 1928) 
and using values of the principal moduli determined by Durand (1936) 
a value of Gi449)= 1-22 x 10 dynes/cm? is obtained. Putting ca seas) fe 
b=3x 10-8 cm and oy = 3-3-5 x 108 dynes/cm?, we find that r lies between 
2 and 11 microns. Figure 6 is a possible example of this phenomenon. 
A half loop source at S, lying some 8 microns from U, was activated by 
the stress from piled-up dislocations in slip band V and the dislocations 
thus produced then cooperated with the primary dislocations to produce 
a crack at U. In this example the shear stress at 8, given by 
og=fo,(L/r)? (Stroh 1957) 

where f is a numerical factor approximately unity, o, the resolved shear 
stress on the primary slip plane and L is the length of the primary slip 
line, was some 20c, and thus was more than sufficient to activate the source. 
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Figure 7 illustrates schematically a situation similar to that in fig. 6 
and suggests one way in which grain boundary cracks may be formed 
at low strains. In the diagram, specimen ABCD is compressed on 
faces AD and BC. A suitably oriented source S, in grain I is activated 
by the resolved shear stress and emits negative edge-dislocations in the 
direction of the boundary which become piled-up at G. The stress-field 


Fig. 7 


associated with this primary pile-up then activates a second source, S,,, 
situated an appropriate distance away from the boundary in grain II 
(Stroh 1957), and oriented so that it emits dislocations of opposite sign 
(positive) to those in the primary pile-up in the direction of the boundary. 
een ener es ae 

} It is considered unlikely that a source would be activated which would 


emit dislocations in the direction of the boundary of th : 
in the primary pile-up. y e same sign as those 
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Because the activation stress on §,, is expected to be large, possibly 
20c,, this source is likely to produce an avalanche of dislocations (Johnston 
and Gilman 1960). These will also pile up at the boundary forming a 
secondary pile-up displaced some distance away from the primary pile-up, 
for example at H. 

Now if H lies between G and , fig. 7 and, from the present metallo- 
graphic observations, if GH is not more than one or two microns in 
length, then a crack is likely to be formed in the boundary between the 
two promoting slip bands. This will be a consequence of the cooperation 
of the two piled-up groups of dislocations in producing a high tensile 
stress over that small area of boundary, GH, between them. Furthermore, 
coalescence of leading dislocations in the two pile-ups also favours crack 
formation in this area. 

However, if H lies between G and F, then it is less likely that a crack 
will be formed for the pile-ups now cooperate to put the boundary 
between them into compression. Moreover, any cracks formed by 
coalescence of piled-up dislocations will be separated by the distance 
GH instead of being located within GH. 

Stroh (1957) has suggested that if a crack is to form near a group of 
piled-up dislocations, then it is necessary for the material round the group 
to remain unyielding ; for plastic flow will tend to relieve the stresses 
associated with the pile-up and in this way prevent the initiation of a 
crack. The present observations suggest that the converse is sometimes 
true in magnesium oxide for, under certain conditions, edge dislocations 
on each side of the boundary may ‘ cooperate ’ to produce a crack. 

In fig. 1 at A and in fig. 2 at K,, K, and Kj, cracks have formed between 
slip bands believed to contain edge dislocations of opposite sign and which 
were suitably displaced from each other at the boundary. Taking into 
account the direction of the applied compressive stress, the nature and 
position of these cracks appears consistent with this view. However, 
other experiments are planned in which an x-ray technique will be used 
to examine the veracity of this opinion. 

The investigation has demonstrated that grain boundaries in magnesium | 
oxide are not inherently weak in compression but that cracks can be 
formed at the boundary by the action of piled-up groups of edge 
dislocations. ; 
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ABSTRACT 


In an attempt to extend the Edwards—Matthews (1957 a) theory of 7—N 
scattering to higher relativistic energy it was discovered that their evaluation 
of the fourth-order Feynman graphs required was incorrect, involving an 
expansion which fails to converge at the threshold of physical scattering 
energy. A new evaluation of these graphs is then carried through and the 
result substituted into their scattering formula. The new values do not 
show the agreement with experiment found in Edwards and Matthews’ 
original results. A study of the Feynman graphs and the method of their 
evaluation are briefly described in § 2. 


§ 1. Backerounpn, Discussion AND SUMMARY OF RESULTS 


THE most prominent feature of 7—N scattering at energies of a few hundred 
Mev is the sharp resonance-like behaviour at around 190Mev. The 
resonance has been assigned with great certainty to the 7’'= 3, J =3+state. 
Another surprising feature is the smallness of the s-wave scattering at 
low energy. A straightforward weak coupling perturbative approach 
failed completely to reproduce either of these features and recently several 
theories have been proposed to try to predict them from accepted field 
theoretical principles. The Edwards—Matthews (1957 a) approach is one 
of these. 

In an earlier paper, Edwards and Matthews (1957b) had proposed a 
very simple model to the problem and they obtained the Chew—Low 
(1956) plot. The work under discussion was an attempt to extend the 
same sort of treatment to a relativistic system. Starting from the 
Schwinger (1951) Functional Differential Equation for the single-nucleon 
propagator, they found an approximate soiution for the scattering ampui- 
tude in terms of the first and second Born approximations of the type that 
Cini and Fubini (1954) had proposed for a non-relativistic system, namely 

os” 
where S,“, S,® are the first and second Born approximations to the 
S-matrix projected into the various partial-wave states a. Readers are 
referred to the paper of Edwards and Matthews (1957a) for details of 
their derivation of this formula. 
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Edwards and Matthews did not mean the derivation to be rigorous 
but merely as a guide, and several rough approximations had been made 
which were not really justified. But the agreement of their final formula 
with experiment was extremely good. Not only did they predict the 
33-resonance at the right energy and of the right shape but they also gave 
the correct order of magnitude of the s-wave and other p-wave phases. 
This is remarkable since they started with the fully relativistic and re- 
normalizable y,-coupling. There being no cut-off parameter, they had 
only one free parameter to fit the experimental results and this came out 
with the value g?/47 216 in agreement with the photo-meson production 
data (Kroll and Rudermann (1954)) and dispersion theory (Cziffra et al. 
(1959)). This compares favourably with the Chew—Low (1956) theory, 
which requires two parameters, g and a cut-off parameter. Also Chew 
and Low started with a o. V interaction and their theory thus gives no 
s-wave scattering whatever, whereas experiments show some, though 
small, s-wave scattering. 

The application of the Edwards—Matthews scattering formula involves 
the calculation of the renormalized Feynman graphs of fourth-order in 
the coupling constant. These Edwards and Matthews calculated by an 
expansion method which will be briefly described and discussed in § 2. 


Fig. 1 


The present investigation was started with the intention of extending 
the Edwards—Matthews calculation to higher relativistic energy. The 
interest was mainly in the application rather than the justification of the 
formula. In the course of the investigation, it was found that the 
Edwards and Matthews evaluation of the fourth-order Feynman graphs 
was incorrect and involved an expansion which would fail to converge at 
the threshold of physical scattering energy. It might also be recorded 
here that there exists another calculation of these same graphs by Wyld 
(1954) which is also incorrect, for the result does not have the crossing 
symmetry. 

We are thus led to a detailed study of the properties of these graphs and 
a method of their evaluation. Apart from the object of testing the 
Edwards—Matthews theory, it seems that a correct evaluation of these 
graphs is valuable though they have no direct physical significance. The 
Investigation is unfortunately rather lengthy, but methods were found 
whereby Self-Energy and Vertex-Diagrams can be evaluated exactly 
except for an almost trivial approximation in the vertex part. The 
remaining diagrams (fig. 1) are more difficult but an expansion method is 
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found which is expected to converge for energies for which k?/m?<1. The 
calculation was carried through and the result substituted into Edwards 
and Matthews’ formula. The result, unfortunately, differs widely from 
their original values. We still have a resonance in the 33-state the 
position of which may be fitted to the experimental value by taking 
g?/47~10. But for this value of the coupling constant g, the 11 and 13 
states both show resonances at laboratory scattering energies of 270 Mev 
and #700Mev respectively. The 11 resonance at 270Mev is clearly 
inadmissible, and though the z—N scattering curve does show a T=4 
resonance at *615Mev this has been assigned to J=3— from angular 
distribution data. We claim therefore, that Edwards and Matthews’ 
agreement with experiment was accidental. 

It is believed that the method of evaluation of the Feynman diagrams 
developed in this work and their properties may be useful to some readers 
and the main steps are briefly described in the following section. 


§ 2. RENORMALIZED FEYNMAN GRAPHS OF THE FouRTH ORDER 
IN THE COUPLING CONSTANT 
The Lagrangian on which the Edwards—Matthews theory was based 
was the pseudo-scalar Lagrangian 


L=L,—ghyystby -D., . . ° . ° (2.1) 
where vector-notation was used for isotopic-spin space. The diagrams 
under discussion are those given in fig. 2 (a), (b), (c), and the corresponding 
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crossed ones (a’), (b’), (c’). It is often convenient in what follows to 
distinguish between the masses of external particles im for nucleons, j a: 
mesons) and internal particles (m’ for nucleons and ,’ for mesons). The 
aphs are obtained by putting m’=™m, p’ =p. 
Sea note that the Aaa (a) and (b) are functions of one scalar 
variable only which may be taken to be s= — (p+k)? for the uncrossed 
diagrams and §= —(p—k’) for the crossed ones. The diagram (c); ig08, 
function of two variables, both s and §. Because of the crossing sym- 
metry, the crossed and uncrossed diagrams are actu ally the same eae 
of the scalar variables s and § but with their roles interchange - Thus, 
e.g. (a’), the crossed Self-Energy diagram is, apart ae isotopic spin 
dependence, the same function of § as (a) the uncrossed vedere Pe i 
The physical values are to be taken for the appropriate values of the 
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variables. s being the total energy, has physical values > (m+)? and 8, 
the four-momentum-transfer, has physical values <(m—,p)?. 

Next we note that the uncrossed diagrams (a) and (b) contribute only 
to the partial wave 7 =} and J=}. This is because we have as an inter- 
mediate state a single nucleon. The diagram (c) and all the cross diagrams 
contribute to all partial wave states. 

The matrix elements represented by the diagrams are given as integrals 
over momentum space. These integrals can be converted into integrals 
over auxiliary variables by the Feynman method. Also the diagrams («) 
and (b) are divergent and these have to be renormalized following the pre- 
scriptions of Dyson. The method used in this part of the calculation is 
familiar and will not be presented. 

The resulting integrals are of the kind typical of perturbation theory, and 
are of the form 


O(¥y 22 =) : 
a. < ee eee 
| late oi 
where ¢(z,....) is some polynomial while a(x,....) is a quadratic of the 
variables w,.... The number of auxiliary variables necessary is 1, 2, and 


3 respectively for the graphs (a), (6), (c). The analytic properties of these 
graphs are readily studied by studying the zeros of the quadratic denomi- 
nators a(z,...). Techniques for this are mainly algebraic-geometrical 
and are now highly developed by several authors (e.g. Landau (1959)). 
Only the most primitive techniques are necessary in our very simple 
cases here. The zeros of a(x,....) lie on a quadric surface in general, the 
position and shape of which depend on the variables s and §. The values of 
s (or 8) at which the quadric surface a=0 enters the region of integration 
of the auxiliary variables x,... usually correspond to branch-points. 
In the graphs under consideration these branch-points are square-root 
branch-points. We shail not give details of this investigation but shall 
only give the result. 

(a) and (b) are analytic functions in s on the cut s-plane, cut from 
s=(m+p)? to oo, with the square-root branch-point at s=(m+p)*. 
The crossed diagrams (a’) and (b’) have the same analytic properties in 
the variable 5. Note that there is no pole term in spite.of the single-nucleon 
intermediate state. This has been absorbed into the coupling constant 
in the subtraction procedure of renormalization. 

The analytic properties of (c) are best stated in the two variables s, and 
t= —(k—k’)P=2m? + 2u?—(s+8). For fixed ¢, (c) is analytic in s on the 
cut s-plane cut from (m+ )® to oo. For fixed s, (c) is analytic on the 
cut t-plane cut from t=4m2 to co. The crossed diagram (c’) has the 
same analytic properties in the variables § and ¢. 

_These analytic properties are well known now since the intense work on 
dispersion relations and the subsequent work of Mandelstam (1958) and 
others. They are just the usual threshold branch-points. It is however 
exactly because of these branch-points that the Edwards—Matthews 
expansion method breaks down. They have expanded the diagrams 
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about w=0 in the unphysical region, (w*=k*?+ 7). Because of the 
branch-point at w=, the expansion fails to converge at just the physical 
energy of scattering. We shall give a more specific example. In the Self- 


Energy diagram, we have a typical integral of the form 


» 


if a2 —te 
I= | dx log — eee ees .. (2:3) 
0 Gm 
where 
A? = sx? + (m?—p?—s)u+p, 
eR sy ty? — eae ae 
An = ML — we 


Edwards and Matthews expanded the integrand into a power series in w 
and integrated term by term, keeping the first non-vanishing term in a 
sum of such integrals. Such an expansion is justified only if the series 
is uniformly convergent in the region of integration. It can be seen, 
however, that uniform convergence fails just at threshold. The expansion 
method is thus invalid. In fact it gives a real value for the diagram 
whereas above physical threshold it should necessarily be complex. 
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Finally it might be pointed out that from the considerations above, at 
can readily be seen that the physical values of ee crossed diagrams (a’), 
(b') are real, while the uncrossed ones are compiex. A convenient way to 
i is is the following. 

Baise the uncrossed Ect diagram gives a diagram which can ede 
for physical energies, hence the existence of a threshold ee ae an 
complex values. Cutting the crossed vertex diagram gives a diagram 
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which requires a free nucleon to decay into a nucleon and two pions which 
is kinematically impossible in physical region. Hence we have real 
physical values for (b’). These considerations may be applied to more 
complicated diagrams. 

We next come to the task of actually evaluating these diagrams. We 
shall not give any details of the lengthy calculation but shall only briefly 
describe the techniques evolved. Two general points to note are the 
crossing symmetry and the infra-red divergence. Because of the crossing 
symmetry, the crossed and uncrossed diagrams are the same analytic 
functions of s and § as pointed out above. Hence any analytic evaluation 
of one also holds for the other. It was found convenient always to keep 
to these variables so as to make the crossing relations explicit. Such 
seemingly simplifying procedures as choosing on the onset a special co- 
ordinate system tend to make the evaluation more complicated. For 
the same reason, it is more convenient to leave partial wave projection, if 
required, until the end of the evaluation. Secondly, we shall in what 
follows neglect the external meson mass yp. This gives errors of order 
u2/m? which for the pion—nucleon system is only about 2%. We must, 
however, not neglect the terms involving »’, the internal meson mass, for 
this will give the well-known infra-red divergence for zero-mass mesons. 
Even when the infra-red divergence is not apparent, the neglect of .’ 
often leads to errors of order log ,.2/m? or higher. 

The Self-Energy diagrams can be reduced to integrals of the typical form 
given in eqn. (2.3). All these integrals can be integrated exactly. We 
shall give the integrals J as example: 


L=Q(s)—Q(m*), 


Wit —— py Les 
= 
2s 


m2 
log + R(s), 
fad 


l Pie ea A te 
Riya uh lone L a pe+s 2 
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A= [(m? — py? —s8)? — 4sp?]12 si) 80, ny Rate eyelet tated 


where the value of the logarithm is to be taken on the principal sheet for 
real values of s<(m+)?. This is readily seen if we follow carefully the 
Feynman prescription (given by ie). The integral has a square-root 
branch-point at s=(m+ )?. The seeming branch-point at S=(m—p)* 
cancels off on the ‘physical’ sheet. Related integrals have similar pro- 
perties and will not be separately discussed. 

A typical integral in the vertex diagrams is of the form 


1 lee 1 
J=| da| yee 
0 OC aave 


ay? =a,7+ K(1—a)y, ) 
Rete ee, 


where 
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Terms of order of the external meson mass have been neglected. After 
integrating over y, we have 
a,” —te 

eee ced es (2.8) 


m 


] 
lo °8 a 


This integral cannot be done in terms of elementary functions but can be 
reduced to the Spence Function defined by 


fle) = [dy SC Seve ae ers. (2:77) 


This function has been studied and tabulated by Mitchell (1949), to whose 
work we refer the reader for more information. With the knowledge of 
the Spence Function, the whole vertex diagram can be evaluated. 

The diagrams (c) are by far the most difficult to evaluate, not only 
because of mere complication but also because the integrals are functions 
of two variables s and §. A typical integral is of form 


Pl x y 1 
L=| dar| dy | de Ste eros Ea 
0 0 O= (477346 
where 


a2 = mx? + w2(1—x)+ K(1—ax)y+{(K+K')y—K’x— Kz, 


K=m?—s, K’'=m?—8. 


Terms of order of the external meson mass have again been neglected. 
Integrals of this type were calculated in this work by an expansion method. 
For K’+K=0, the integral Z and related integrals can all be reduced 
algebraically to integrals already encountered in evaluating (a) and (d). 
Now Lis analyticin K’ at K'+ K=0. We may thus expand in K’ about this 
point. The nearest singularity on the K’-plane is at t= — (K+ K’)=4m?. 
Since t=2k?(1—cos6@), where k is the centre-of-mass momentum, 
we expect the expansion to converge uniformly for all physical values of 
cos@ for k2<m?. This criterion for convergence is readily satisfied for 
laboratory energies of up to several Bev. Thus 
)=L(K 2(c0s 8-1) 2 eh ae Be 
L(K, K')=L(K, — K) + 2k?(cos 6 ae 1 een Fe (2.9) 
Similarly we may expand the crossed diagram in K about K + K =O The 
expansion coefficients can then be reduced after a rather tedious calcu- 
lation to known integrals. In the actual calculation the expansion has 
been carried to the second term so that p-wave contributions are taken 
consistently to the leading order. 

It may be remarked here that in all the calculation above only terms 
of order p2/m? have been neglected, the internal meson mass has been 
retained throughout. The same method of evaluation can thus be applied 
without alteration to diagrams where the virtual particles are different 
from the incoming ones; e.g. strange particles, whose masses are too large 


to be neglected. 
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The projection into partial-wave states follows standard procedure 
(see, e.g. Jacob-Wick (1959)). Uncrossed diagrams are readily handled 
since the integrals depend only on the energy and not the momentum- 
transfer. For the crossed diagrams since the integrals are functions of 
cos 8, the projection into angular momentum states had to be done numeri- 
cally. This is, however, not difficult since the functions are regular 
everywhere for physical cos 6. 

The implications of the result of the calculation on the Edwards— 
Matthews method have already been discussed in §1. A note containing 
the complete result of the evaluation and some intermediate steps is 
deposited in the Department of Mathematical Physics, University of 
Birmingham. Duplicates of this are available on application to the 
department. 
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ABSTRAOT 


The collapse of narrow V-shaped notches formed in glide dislocations at 
points of close approach to dislocations threading the glide plane is shown 
to facilitate intersection without thermal activation. The force effecting 
intersection, which is a component of the line tension of the notches, sub- 
sequently assists the glide dislocation to overcome the drag of jogs formed in 
the process. A consequence of this dual function of the force is the Cottrell— 
Stokes law, 1.e. the strain invariance, at any given temperature and strain rate, 
of the ratio 7p/7g, where 7p and 7g are the temperature dependent and 
independent parts of the flow stress respectively. Linear and parabolic 
work-hardening laws are deduced from the model. On equating the slip 
distance to grain size in the parabolic relation the flow stress is found to be 
linear in (grain size)—1/?, which is known to apply for example to zinc. The 
activation energy of the jog migration associated with 77, estimated from low 
temperature flow stress data, is equal to about 4 ev in copper. 


§ 1. [yTRODUCTION 


In his theoretical treatment of work-hardening Seeger (1954) suggests 
that the flow stress of a pure metal, 7(7Z'), should be resolvable into two 
independent parts at low temperatures : 

(T)=TetT7, Cys, ae AS eee CL) 
where 7, depends upon temperature only as the shear modulus G', while 
Tp increases approximately linearly with decreasing temperature, and 
depends logarithmically on the strain rate. Cottrell and Stokes (1955), 
and Cottrell and Adams (1955), who investigated the flow stress of alu- 
minium and copper crystals, found that the ratio of the flow stresses at 
two temperatures 7(7',)/7(Z',) was the same at any given strain, implying 
a proportionality between 7, and t7. The relation was subsequently 
confirmed to hold also in polycrystaline metals (Basinski 1959), as well 
as in precipitation hardened alloys (Kelly e¢ al. 1959). The object of the 
present work was to study this proportionality, and to attempt to elucidate 


its significance. 


§ 2. THe Distocation MopEL 


We shall relate our discussion to the dislocation configuration shown in 
the figure (a). The dislocation, in the form of circular ares, is held up 
at lattice defects such as J; the precise nature of the defects will be con- 
sidered later. For simplicity the dislocation will be assumed to be un- 


dissociated and of the edge type; the generality of the argument is not 


+ Communicated by the Author. 
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thereby affected. In order to migrate forward the dislocation must cut 
through the obstacle or drag it along. An activation energy, Q(T), 
which may be a function of temperature, would therefore be required. 


Early (a) and (5) late stages in the process leading to the intersection of the 
forest dislocation J; (c) idealized form of the cusp prior to collapse. 


Now, if U is the line tension of the dislocation, and @ the angle between 


the tangent at J (figure (a)) and the direction of the Burgers vector (IF), 
the force on the defect is 


F=2U cos@. 


Expressed in terms of the curvature of the are 


U =pb(r—7T4), a ee Re ane) 
where, as in Seeger’s model, rg is. the temperature independent elastic 


stress opposing the applied stress, and 6 is the Burgers vector. As 


cos =1/2p, with / equal to the mean spacing between pinning points, we 
have 


F=(r-7,)bl Sets ee a ES 
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and hence 
Q(0) = (7 —74)07. 
If flow is to take place at an observable rate then 
OTL)=(t#—Te)blamkT, . 2. s wae (4) 


where m is a constant approximately equal to 30; the actual value chosen 
will determine the flow rate at which the deformation takes place. 
On writing 


WT) =OL) nk, * 0 ele renee (5) 
the flow criterion becomes 
bir 7)= WT) 2. a een (6) 
which, on putting 
Ta=aGb/l, oer) eee ak ER 
yields 
TWAT) [O° eeG Diba ion oe tA eee nen (S) 


For the present we take « to be a constant of order 0-1-1; a more precise 
estimate will be obtained below. Fluctuations of internal stresses will 
be small in a work-hardened crystal, for at spots where « differed 
appreciably from the most probable value, e.g. where it is smaller, - 
work-hardening would tend to occur preferentially, largely removing the 


discrepancyy. 
On comparing eqns. (1) and (8) one finds 
tp = W(T)/b*4, si Seaeticat eon (Oo) 
so that the ratio 
T q(T = ab? /W(T) i cee cam heer etabl 0, 


is independent of /, and hence of the state of work-hardening of the crystal. 
The ratio of the flow stresses at two temperatures, obtained with the same 
specimen without intermediate recovery or deformation, i.e. with / constant, 
is also independent of the strain, as is clear from eqn. (8). The 
Cottrell-Stokes law is therefore obeyed by the modelt. 


§ 3. WorK-HAaRDENING Laws 
Specific relations between / and the strain y have to be assumed in 
deducing the work-hardening law of the metal. To do this we shall take 
as model a crystal cube of side L, with one set of parallel active slip planes 


+ This argument would not necessarily apply under conditions of * easy- 
glide ’, when only a few widely spaced slip planes are active. In fact Seeger 
et al. (1957) and Basinski (1959) find deviations from the CottrellStokes 
relation at low strains. ere 

+ The Cottrell-Stokes law is in fact implicit in power laws of the type 
r=K(T, )(y—yo)™ irrespective of the magnitude of the exponent m, for the 
ratio 7(7',)/7(T'2) Af(y) at y= const. 

02 
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d apart, each plane containing a group of equidistant dislocations of equal 
sign, with a mean spacing J. The total shear strain due to the dislocations 


is then 


y= (L?/Id)(4b/L) =bL/2ld, enh Sommirimited tei 2 Gt 
so that 
L=bL/2yd, = Bl Reh eo lae ee ee 
Or, t=a: 
T= (bL/2y)!. ol he ae ORD ee eth 


Equation (12), with d>/, should be applicable to single crystals or well- 
annealed polycrystals in the early stages of deformation, while eqn. (13) 
would be expected to represent the relation between the strain and the 
dimensions of the sub-structure more correctly at strains at which more 
than one slip system have become fully operative. On eliminating / 
between eqns. (8) and (12), and between (8) and (13) respectively, one 
obtains the linear and parabolic work-hardening laws 


tT=X(T)y, at eer eee ee PN 
with 

A=[W(T) + «Gb? |(2d/ Lb), . lap ek. Le eee 
and 

= l(T) yt wi ke are eee) 
where 

C=(W(T) + aGb?](4b°) 7242, ee i 


If the slip distance is equated to the grain size, then the flow stress given by 
eqn. (16) is linear in (grain size)—"*, which is known to be the case in zine 
(Zein 1955) and molybdenum (Johnson 1959). Subgrain dimensions 
(Ball 1957), or even smaller elements of substructure (Feltham and 
Meakin 1957) may determine L in face-centred cubic metals. In eqn. 
(15) L and d, but not L/d may vary with the strain. 


§ 4. THE Activation Enercy Q 
At 42° the ratio 74/7 in copper crystals is about 10, (Basinski 1959), 
so that eqns. (5) and (10) yield, on taking 7'=0, 
Q=0-laGb? (ab 42K), 2. . . . . . (18) 
As Gb? ~ 5ev in copper we have 
OS 0GKeT! 6). sn ieedeenad teat 


Now, a further estimate of @ can be made by considering that the pro- 
nounced temperature dependence of the flow stress in copper on lowering 
the temperature begins to be observed at about 150°K; at still lower 
temperatures W(7’) in eqn. (16) must therefore be positive, i.e. non-zero. 
‘Thus, taking W(7')=0 at 150°K, and m=30, eqn. (5) yields 


Q=0-37ev (abi160°K).. . g au ee (20) 
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If the temperature dependence of Q is assumed inappreciable between 4 
and 150°K, comparison of eqns. (19) and (20) shows that 


OER he lo 20) en Ce 


suggesting that the applied stress at low temperatures is close to the 
‘Frank—Read stress’ Gb/l. This result indicates that the obstacle (I in 
figure 1 (@)) is not intersected by a cutting mechanism of the type considered, 
- for example, by Thornton and Hirsch (1958), which would require an 
activation energy considerably higher than the value of Q indicated above 
(Friedel 1956, Seeger 1955 a, Seeger ef al. 1959). We believe that the 
mechanism by which the glide dislocation passes through the obstacle /, 
which we shall now take to be a dislocation on an intersecting glide 
plane, probably occurs as follows. 

Considering the edge dislocation in the figure (a), we see that on in- 
creasing the stress the arms of the dislocation cusp with node at J will 
reach a critical configuration (figure (6)) at which the segments JP and JQ, 
which are virtually pure screw dislocations of opposite signs, will cancel, 
leaving a small jogged loop around J. If J is a screw dislocation a jog 
will also remain on the glide dislocation. As JP and IQ lie on adjacent 
slip planes, a line of vacancies or interstitial atoms will remain over the 
length over which the screw dislocation segments annihilated. The 
jogged loop around J would shrink under its own line tension, and dis- 
appear, leaving a jog in the forest dislocation J. This process of dis- 
location intersection would not be significantly temperature dependent, 
but it would not lead to flow unless the glide dislocation can migrate from 
points such as J after intersection. Ifa jog is left in the glide dislocation 
after intersection then the subsequent migration will be assisted by thermal 
activation of the movement of the jog, the latter being subject initially to 
the same force as the forest dislocation J at the instant of ‘intersection’. 
The rate determining process of the flow can therefore be regarded as 
consisting of two stages; the first, which is not thermally activated, 
results in the formation of jogs in (glide) dislocations, while the second 
involves slip of the glide dislocation. The second should be temperature 
dependent, with an activation energy equal to that of jog migration, NES 
numerically equal to that of vacancy formation for screw/screw jogs 
(Seeger 1955b) or somewhat less for jogs in edge dislocations (Maddin 
and Cottrell 1955). Activation energies determined from flow stress and 
logarithmic creep data by Thornton and Hirsch (1958) are 0-7 + 0-35 ev 
for copper at 80°K; Feltham (1961) associates about ev with the 
migration of jogs in edge dislocations in copper, with somewhat higher 
values for alpha-brasses. These experimentally determined activation 
energies are therefore comparable with the Q's deduced above (eqn. (20)) 

A further inference of the anlysis relates to the dislocation configuration 
near I (figure (b)) just prior to collapse of the cusp; eqn. (21) indicates 
that a narrow cusp is to be expected, with PIQ only a few degrees. Now, 
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the approximate shape of the cusp as the flow stress is being attained can 
be established from the following considerations. 

We assume PI and QI to be coplanar screw dislocations of opposite 
signs, and take for the stress between them 


7(8)=G@b[27(280+w), . . - . = = =) (22) 


where s is the distance of a point along JQ (or JP) measured from J along 
the curve, 6 is now taken to be the angle of a notch with straight sides 
obtained by joining JP and IQ as shown in the figure (c), and w is a length — 
of the order of the diameter of the cusp at the node. Taking the applied 
stress to be aGb/l, the radius of curvature is 


p(s) =(U/b)/[aGb/l + Gb/2a(2s0+w)], . . . . (28) 
or, with U=46G0?, 
L/p=2a/l + Ifw(2s0@+w). : . . . . . (24) 


The curvature p~ is thus high for small s+. The interaction will therefore 
be negligible beyond a small region of the cusp close to J, and we shall 
confine our attention to a small zone P/Q, approximating to the arms of 
the cusps by straight segments bent sharply at P and Q figure (c). Again 
measuring s along JQ, the work done on the field between the dislocation, 
in moving them from PIQ, to P’IQ’ as shown in the figure (c), is 


2 
Weed tore [S“ Sw (29) 
which, omitting w in the denominator of the integrand, gives 
dH (0)= —(8;gGb?/270)d@. 2. 2. 1 1... (26) 
The associated increase in the line energy of the cusp is 
dE(U)=2Us;9d0=Gb*s,,d0. . . . . . . (27) 
The criterion of instability 
a, B(0)+H(G)/d@=0 ~~ ie « © (28) 
then yields 
D cay eel | Derk lsat . en 


The ‘clipping-off’ process of intersection will therefore occur when the 
angle P/Q has diminished to about 18°. This estimate agrees with the 
implications of eqn. (21). 


ee ee eee eee 
} On using the intrinsic equation dy/ds=1/p, eqn. (24) is easily shown to 
represent the spiral 


y= 2as/1+ (1/6) In (14+2s0w-1), 
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ABSTRACT 


Several kinds of etch pit are formed when NaCl is etched with a mixture of 
acetic acid and methyl alcohol. In particular, slip bands contain shallow . 
pointed pits in addition to the normal pits formed at dislocations. The 
density of these shallow pits is reduced by an anneal at around 300°c before 
etching, and it is believed that they originate at ‘ debris ’ left behind moving 
dislocations. Debris is also produced by the to-and-fro motion of edge and 
screw dislocations in crystals vibrated at high amplitudes at 90 ke/s, and this 
effect is related to a decrease in the internal friction. Experiments on LiF 
show that three etches which reveal dislocations differ in the extent to which 
they also reveal debris. 


§ 1. [INTRODUCTION 


Tue technique of etching ionic crystals to reveal dislocations is now well 
established (see for example Gilman and Johnston 1957); and the present 
paper describes the results of a series of etch-pit observations on crystals 
of NaCl. The main feature of these experiments is that the etch produces 
pits of various types, only some of which can be attributed to the presence 
of normal dislocation lines. Slip bands in crystals etched after static 
deformation or vibration at 90 kc/s have often been observed to contain 
pits which do not appear to be formed at dislocations ; the presence of these 
additional pits suggests that under certain conditions ‘ debris’ is left 
behind moving edge or screw dislocations. 

In the case of LiF somewhat similar results have been reported by 
Johnston and Gilman (1960), who observed indistinct pits on the slip 
planes of expanded dislocation loops. They attribute these pits to debris 
produced by moving dislocations, and Keith and Gilman (1959) suggest 
that this debris consists of point-defect clusters. The extent to which pits 
are formed at debris depends on the etch used, and in this respect the 
etches of Gilman and Johnston are not very effective. 


§ 2. ExpERIMENTAL TECHNIQUES 


The experiments were performed on specimens cleaved in {100} orienta- 
tion from single crystals grown from Analar grade N aCl by a modified 
Kyropoulos method. These specimens were annealed at 650°C to remove 
ie EE EE ees 
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the internal stresses produced by cleavage, and then cooled to room tem- 
perature over a period of several hours. 

The specimens were etched in a mixture of 1 part of methyl alcohol 
with 2 to 4 parts of glacial acetic acid for 10 to 30 sec, washed in anhydrous 
ether and dried in a stream of cold air. This method was found to be more 
satisfactory than those described by Amelinckx (1954) and Moran (1958), 
which required very short etching times. 

A chemical polish was used to remove the surfaces of specimens so that 
internal surfaces could be etched. It consisted of a mixture of 1 part of 
methyl alcohol with 2 parts of ethyl alcohol, and removed about 0-7 per 
min leaving the surface covered with extremely shallow overlapping 
pits. 

The etch pits were photographed in oblique incident light, except in a 
few cases where incident phase-contrast microscopy was used. 


§ 3. OBSERVATIONS 
3.1. Hitch Pits in Undeformed Crystals 


Figure 1} shows part of the etched as-cleaved surface of an undeformed 
crystal. The pyramidal pits have their edges in (100) directions, and the 
faces of most of these pits make an angle of about 10° to the crystal 
surface. However, some of them are shallower than this, and a similar 
effect has been reported by Amelinckx (1954). It is believed that all 
these pyramidal pits are formed at dislocations, and this is supported by the 
observation that there is a one-to-one correspondence between such pits 
on the matched cleaved faces of undeformed crystals. The etch often 
produces flat-bottomed pits along cleavage steps, as shown in fig. 1, and 
these pits may originate at very small dislocation loops formed at the step. 
Similar flat-bottomed pits also form at random sites, particularly on 
annealed surfaces. 

In addition to pyramidal and flat-bottomed pits, a high density of very 
shallow poorly resolved pits was observed in a specimen that was 
x-irradiated to produce some permanent coloration and then polished and 
etched. Pits of this type have been observed by Aerts et al. (1959). who 
attributed them to point defect clusters produced during the irradiation. 
Similar unresolved pits are formed in crystals containing 0-02% CdCl,. 


3.2. Statically Deformed Crystals 
3.2.1. Dislocations produced at indentations 
A convenient way of introducing a controlled number of dislocations 
into a specimen is to deform it with a microhardness indenter; subsequent 
etching produces a ‘rosette’ of pits, and fig. 2 shows part of such a 
rosette. Practically all the pits formed in rosettes are identical, both on 
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external surfaces and on surfaces polished after indentation; these pits 
are similar to the majority of those attributed to dislocations in undeformed 
crystals. When indented and etched specimens were strained slightly to 
move the dislocations and re-etched, many flat-bottomed pits were 
observed together with similar numbers of smaller pointed pits. These 
observations indicate that the pointed pits originate at dislocations, and 
they will hereafter be called normal pits; it seems reasonable to assume 
that all fresh dislocations should produce this type of pit. 


Many experiments have been performed to study the pits formed by 
etching crystals which had been deformed until a few slip bands were 
produced. The surface of the crystal present at the time of deformation 
will be called an external surface, and one subsequently exposed by cleav- 
age or polishing an internal surface. A slip band in a (110) direction 
where pits are formed at edge dislocations will be called an edge slip 
band, while one in a (100) direction with pits at screw dislocations will be 
called a screw slip band. 

Figure 3 shows etch pits along part of an edge slip band on an internal 
surface, which was exposed by cleavage and then etched. It is clear 
that there are at least two types of pit in this band: normal pits, and 
others which are shallower and somewhat smaller. The proportion of 
normal to shallow pits varies considerably, there being rather more shallow 
pits than average in fig. 3. There are, however, almost always some shallow 
pits in slip bands on internal surfaces, unless the density of normal pits is 
so large that the others are masked. The density of shallow pits on 
external surfaces is usually much smaller than that on internal surfaces. 
Both normal and shallow pits sometimes become flat-bottomed during a 
single etch, and in some cases they start to grow again at a point slightly 
displaced from the centre of the original pit. This intermittent etching 
is similar to that described by Johnston and Gilman (1960), and some 
examples are shown in fig. 4. 

Relatively few shallow pits form on slip bands in crystals that have been 
annealed at 200-300°c subsequent to deformation. To demonstrate this 
effect a specimen was lightly deformed, cleaved in half to expose an internal 
surface, polished for 1 min and etched. It was then annealed for 1 hr at 
250°c, furnace cooled, polished for 1 min to remove all the etch pits, and 
re-etched. Figure 5 (a) shows a region of the cleaved surface after the 
first etch, in which several edge slip bands are visible with normal and 
shallow pits along them. The same region photographed after the second 
etch is shown in fig. 5 (b), where the shallow pits are considerably fewer 
in number and less distinct than in fig. 5 (a). . 

Similar investigations have shown that the number of shallow pits of 
this type formed on a given surface can also be reduced by irradiating the 
surface with unfiltered 50 kv x-rays until some permanent coloration is 


produced. 
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The different appearance of normal and shallow pits and their different 
behaviour in the annealing and irradiation experiments indicate that they 
do not originate at the same type of defect. The normal pits can fairly 
reliably be attributed to ordinary dislocation lines or to pairs of edge 
dislocations of opposite sign as proposed by Johnston and Gilman (1960). 
The origin of the shallow pits, however, must differ from that of the 
normal pits, and they are probably formed at some kind of ‘ debris ’ 
left on the slip planes by moving dislocations. The effect of annealing or 
irradiation must be to disperse this debris, or to modify it so that it is no 
longer attacked by the etch. 


3.2.3. Screw slip bands 


Pits other than normal etch pits are formed on screw as well as edge 
slip bands, but in this case the situation is more complicated. When the 
external surface of a lightly deformed crystal is etched, screw slip bands 
are found to contain pits of many depths and sizes ranging from normal 
pits down to ‘streaks’ comprised of very shallow pits lying in (100) 
directions. Part of such a slip band photographed in oblique incident. 
light is shown in fig. 6 (a), and fig. 6 (b) shows the same area under phase- 
contrast illumination, which reveals the streaks more clearly. After 
polishing or cleaving through an existing slip band, the streaks formed on 
etching are fewer and much less distinct, although both normal and shallow 
pointed pits are still formed; figs. 6 (c¢ and d) show the same areas as 
figs. 6 (a and b) after polishing for 1 min and re-etching. The dislocations 
revealed as normal pits have moved during handling between the two 
etches, and there are indications that some dislocations have also moved 
during etching. 

In contrast to the streaks, the shallow pits observed on screw slip bands 
occur both on external and internal surfaces; they will be attributed to 
debris, and this may be related to that giving rise to shallow pits on edge 
slip bands. The density of shallow pits is considerably reduced in 
specimens etched after annealing at 250-350°c subsequent to deformation, 
this temperature range being slightly higher than that quoted for edge 
slip bands. 

The streaks of very small pits on screw slip bands are observed on 
surfaces across which screw dislocations have moved, whether the surfaces 
be annealed, as cleaved or polished, and it is interesting to note that slip 
steps should occur under these conditions. Streaks have not been seen, 
however, in the vicinity of screw dislocations in rosettes. It is not clear 
whether streaks on internal surfaces, like those in fig. 6 (d), arise from de- 
fects present in the crystal, or whether they are associated with movement 
of dislocations across the newly formed surface during handling. On 
several occasions similar streaks have been seen on freshly cleaved and 
etched surfaces in regions where there are no dislocations other than those 
in sub-grain boundaries. An example is shown in fig. 7, where the streaks 


are seen to consist of small pits in (100) rows which terminate at a sub- 
grain boundary and a cleavage step. 
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3.3. Crystals Vibrated at High Amplitudes 

In a previous paper Whitworth (1960) described some changes in the 
internal friction of NaCl crystals vibrated at high amplitudes at a fre- 
quency of 90 ke/s. In particular, these experiments showed that the ampli- 
tude-dependent damping of lightly deformed specimens is reduced by 
vibration at strain amplitudes above about 3x 10-5, Specimens given 
this vibration treatment have now been studied by etching, and the 
observations suggest that debris is produced by oscillating edge and screw 
dislocations; this result probably accounts for the reduction in damping 
by a mechanism involving the pinning of dislocations. Somewhat 
similar observations of the production of debris by dislocations moving in 
LiF under a cyclic stress have been made by Keith and Gilman (1959). 

The results of § 3.2.1 indicate that very little debris is formed around 
indentations, and the specimens used in the present work were therefore 
deformed in this manner. In a typial experiment a specimen was 
polished, indented, and vibrated at a maximum strain amplitude of 
2-6 x 10-4, where the maximum reduction in damping was obtained. This 
crystal was then etched, and a portion of one of the rosettes so formed is 
shown in fig. 8 (a). Figure 8 (6) shows the same area after polishing for 
3 min and re-etching. These figures show that shallow pits are produced 
on the slip planes of edge and screw dislocations, both on external and 
internal surfaces. Streaks are also produced by individual screw dis- 
locations, but only on the external surface. Vibration of crystals at the 
amplitude reached in this case creates new dislocations, as was shown by 
Whitworth (1960), and it is possible that this effect has caused some 
multiplication of the dislocations in the rosette, giving rise to more 
normal pits than are usually observed. 

Figure 9 shows the shallow pits associated with some edge dislocations 
in a rosette in another crystal vibrated at a smaller strain amplitude 
(2:0 x 10-4), where nucleation of new dislocations does not occur; the 
pits shown are on an internal surface, the crystal having been polished for 
3 min between vibration and etching. Shallow pits of this type have never 
been observed at rosettes in any crystal which had not been given a high 
amplitude vibration. | When the specimen of fig. 9 was annealed for 1 hr 
at 300°c, polished and re-etched, very few shallow pits were formed. 

From these observations it appears that a single dislocation, moving 
backwards and forwards on its slip plane under the oscillating stress, 
produces some form of debris in its vicinity. The length of each group of 
pits is of the order of 10-* cm, which is quite comparable with estimates 
from internal friction measurements of the distances through which the 


dislocations should move. 


3.4. Experiments on Lithium Fluoride 


A few experiments on LiF have been performed, to look for evidence of 
debris which could be revealed by etching. The LiF used was supplied 
by Taylor, Taylor and Hobson Ltd. and was harder go probably more 
impure than the crystals from the Harshaw Chemical Co. used by Gilman , 
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and Johnston. Before use, the material was annealed for 12 hr at 800°c 
and cooled over 12 hr. 

Three reagents were used to etch the crystals: the etches A and W of 
Gilman and Johnston (1957), and another etch, ‘H ’, consisting of a 6% 
aqueous solution of hydrogen peroxide, similar to that used by Urusovskaya 
(1958). It was found that only etch H formed an appreciable number of 
shallow pits on slip planes in lightly deformed crystals. To illustrate this 
observation, fig. 10 shows the matched faces of a crystal which was cleaved 
in half after deformation; one half (a) was etched with W, and the other 
(b) with H. It is seen that on edge slip bands H forms poorly defined 
pits which are additional to the normal pits formed by both H and W. 
Experiments to study the effect of annealing on the debris are unsatis- 
factory, because etch H, like etch W, produces poorly defined pits at 
dislocations in crystals annealed at temperatures above 200°C. 

Although etch H is more successful than A or W in revealing debris on 
active slip planes in LiF, it compares unfavourably with the present 
etch for NaCl. Nevertheless, the results provide evidence that debris is 
formed by the deformation of both LiF and NaCl. 


§ 4. DIscUSSION 


The preceding paragraphs describe the production of etch pits in variously 
treated crystals of NaCl. Three main types of etch pit are formed as a 
consequence of plastic deformation: normal pits due to dislocations, 
shallow pits, and streaks. The streaks appear to be mainly a surface 
phenomenon, and as their origin is obscure, they will not be discussed 
further. Shallow pits, however, occur on internal as well as external 
surfaces, and thus show that some defects other than normal dislocation 
lines are left on slip planes during deformation; these defects are referred 
to as debris. It should be noted that very few shallow pits occur in 
rosettes, and the nature of the deformation must therefore be important 
in determining the amount and distribution of the debris produced. 

Some indication as to the nature of the debris can be obtained from the 
shape of the shallow pits, which are pointed and fairly well defined. In 
this respect they differ from the very shallow poorly resolved pits, which 
are formed in irradiated and impure crystals, and believed to be due to 
clusters of point defects. A piece of debris therefore differs from a simple 
cluster of defects, probably because it extends into the crystal for a 
distance at least equal to the depth of its etch pit (~ 0-1 y in fig. 3). 

One possible mechanism for the production of these linear defects is the 
process of dislocation multiplication by the cross-glide of screw disloca- 
tions discussed by Johnston and Gilman (1960). When the distance of 
cross-glide is large, {100} dislocation loops are produced, and these should 
etch-up on screw slip bands. However, small amounts of cross-glide lead 
to the formation of pairs of edge dislocations of opposite sign lying on 
nearby slip planes, and these pairs, or dipoles, probably give rise to many 
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of the normal etch pits in edge slip bands in LiF and NaCl. This belief is 
supported by observations of intermittent etching, although in NaCl it is 
hard to distinguish this effect from the movement of dislocations during 
the etch. When the dislocations forming a dipole are only a few atomic 
spacings apart, their elastic fields largely cancel out, and in the limiting 
case the dipole becomes a column of vacancies or interstitial ions. Such 
close dipoles may etch up differently from dislocations or wider dipoles, 
and could be the origin of the shallow pits in edge slip bands; the depth of 
the pit might then depend on the precise nature of the dipole. 

The motion of jogged dislocations can produce similar trails of defects 
lying in any direction on the slip plane other than that parallel to their 
Burgers vector. These trails should, however, be most stable when the 
density of defects along them is high, and this would favour the (100) 
orientation. Shallow pits may form at such trails, and could thus occur 
in both edge and screw slip bands, although the former will be favoured if 
the stability criterion is important. 

The above suggestions about the nature of the debris are only a few of 
the many which could, no doubt, be proposed. They do suggest, however, 
that the pits in edge and in screw bands may be formed at different kinds 
of debris. Almost any type of debris, such as close dipoles, columns of 
defects or very small dislocation loops, would be expected to disperse on 
annealing at some temperature lower than that required for large scale 
dislocation climb, and this is consistent with the observed annealing 
characteristics of the shallow pits. The results indicate that x-irradiation 
also breaks up debris, either directly or by the motion of the point defects 
which it produces. There does not seem to be sufficient information 
available to make deductions about the precise nature of the debris from 
the annealing and irradiation experiments. . 

No process has yet been proposed which would give rise to debris by the 
backward and forward motion of a pure edge dislocation, and yet debris 
is observed around edge dislocations in vibrated crystals. A possible 
explanation is that the motion of edge dislocations may be rather complex 
in this case, involving the production of kinks in screw orientation which 
can then cross-glide. The formation of large amounts of debris around 
dislocations should interfere with their motion, and the production of 
debris may therefore account for the decrease in damping of vibrated 
crystals that was observed by Whitworth (1960). 


ACKNOWLEDGMENTS 


This work was supported by a D.S.I-R. Grant for Special Research, one 
of us (R.W.W.) holding a research fellowship on this grant, and the other 
(R.W.D.) a D.S.LR. research studentship. We wish to thank Professor 
J. G. Ball for the provision of research facilities in the Department of 
Metallurgy at Imperial College, and the leader of the Research Group, 
Dr. P. L. Pratt, for his encouragement and his criticism of the manuscript. 


224 On Dislocation Debris in Sodium Chloride 


REFERENCES 


Asrts, E., AMELINCKX, 8., and DEKEYSER, W., 1959, Acta Met., 7, 29. 

AMELINCEX, S., 1954, Acta Met., 2, 848. 

Gitman, J. J., and Jounston, W. G., 1957, Dislocations and Mechanical 
Properties of Crystals, ed. Fisher et al. (New York : Wiley), p. 116. 

JOHNSTON, W. G., and GitmaAN, J. J., 1960, J. appl. Phys., 31, 632. 

Kerry, R. E., and Gunman, J. J., 1959, A.S.T.M. Symposium on Basic 
Mechanisms of Fatigue. Special Technical Publication No. 237, p. 3. 

Moray, P. R., 1958, J. appl. Phys., 29, 1768. 

Urusovskaya, A. A., 1958, Kristallografia, 3, 726. 

WuitwortH, R. W., 1960, Phil. Mag., 5, 425. 


Demagnetization of Igneous Rocks by Alternating Magnetic Fields} 


By E. Irvine, P. M. Storr and M. A. Warp 


Department of Geophysics, Australian National University, 
Canberra, A.C.T. 


[Received May 31, 1960] 


ABSTRAOT 


The presence of secondary components of magnetization in basalts often 
limits the usefulness of their directions of natural remanent magnetization 
for palaeomagnetic work. These secondary components may be removed 
by alternating magnetic fields. The procedures for doing this are 
described and tests based on the internal consistency of the results are 
devised to judge their reliability. These are applied successfully to a series 
of results from specimens collected from the Tertiary basalts of New South 
Wales, Australia, and a determination of the geomagnetic field in 8.E. 
Australia during the Lower Tertiary is obtained. A study of the behaviour 
of these basalts in high alternating fields suggests that most of the iron 
mineral grains have a multi-domain character, and contain a hard and 
soft component in the ratio of about 1: 10. 


§ 1. [yTRODUCTION 


In order to obtain the direction of the ancient geomagnetic field from a 
study of rock magnetism it is necessary to measure the direction of the 
primary component of magnetization acquired at the time the rock is 
formed. This primary magnetization is often seriously affected by 
secondary components, which are added later, and which render the 
direction of the natural remanent magnetization (NRM) of little use for 
precise palaeomagnetic work. 

Lavas, on cooling through their Curie point, acquire a primary thermo- 
remanent magnetization (TRM) in the direction of the earth’s field. If 
magnetic particles with relaxation times less than the geological age of 
the rock are present, a secondary component is acquired different in 
direction from the primary one if the direction of the geomagnetic field 
has changed since the latter was acquired. ‘This is viscous magnetization 
and such components have the form of an isothermal remanent magnetization 
(IRM). Rock specimens containing grains with relaxation times in the 
earth’s field of the order of 10? or 10% years will contain components of 
magnetization along the present dipole field. Relaxation times of 1 to 
100 years will produce magnetizations close to the present field, whereas 
shorter relaxation times produce further viscous components dependent 
on the conditions of storage. This latter type of secondary components 
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has been termed temporary components by Creer (1957). Secondary 
components are the most substantial source of error in palaeomagnetic 
studies and may be removed by alternating magnetic fields. The problem 
is to find the alternating field conditions necessary for their removal and 
yet ensure that a measurable component of primary origin remains. 
Because of the variation in magnetic minerals the field strength will not 
necessarily be the same for all lava types even within the same formation. 

In general the removal of unwanted secondary components from 
specimens collected from a lava sequence will decrease the scatter in 
the directions of magnetization. We suggest that certain consistency 
conditions should be fulfilled before the secondary components may 
be judged to have been removed. Compared with the directions of 
NRM (M,) the directions after removal of the secondary components 
should have these characteristics. 


(1) The scatter between specimens taken from the same lava flow 
(within-site dispersion) should either remain constant or should decrease. 


(2) The agreement between the mean directions of different lava flows 
(between-site dispersion) should be much improved. 


(3) Assuming the average dipole character of the earth’s field the mean | 
direction of magnetization for all flows should agree with similar results 
from rock formations of comparable age in the same region after correction 
for differences in geographical location. 


These conditions apply to the great majority of lava sequences which 
were laid down in a period of time during which the geomagnetic field 
had an approximately constant direction (irrespective of sign) different 
from the present field. The first will in most cases be necessary, but it is 
possible to envisage situations in which the within-site scatter increases 
slightly after removal of the secondary component. The second is a 
necessary and the most important condition. The third is also necessary 
given the stated assumption. As a test of these conditions basalt speci- 
mens from several lava flows have been treated in alternating magnetic 
fields (§§ 2, 3, 4 and 5). 

Further experiments on the same specimens were carried out in alternat- 
ing magnetic fields to study the behaviour after removal of most of the 
primary magnetization and it was not found possible to reduce the moment 
beyond a certain limit, called the minimum intensity. This property, 
which may be of some interest in the theory of NRM of basalts, is discussed 
in §7. In §8 an attempt is made to relate the primary and secondary 
components to the ferromagnetic grains present. 


2. DEMAGNETIZATION BY ALTERNATING MAGNETIC Freups 


The procedure is to subject a rock specimen (in the absence of a steady 
field) to an alternating field which isreduced smoothly. The demagnetizing 
instrument used is similar to that recently described by Creer (1959). 
Specimens are placed in a holder which spins inside a set of coils which 
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annul the earth’s field. The alternating field is provided by a solenoid 
mounted on a trolley running on an aluminium tramway. The field at the 
specimen is reduced by withdrawing the trolley to a distance of 3m at 
which point the field is reduced by a factor of 3 x 105. 

Two pairs of square coils (Parry 1957) side length 180-2 and 165-9 cm 
and separations 98-2 and 90-4cm and each with 100 turns annul the 
horizontal and vertical components of the earth’s field respectively. 
Chokes are included in both circuits to reduce the effect of any mutual 
inductance with the solenoid, which in any case should be zero since the 
solenoid axis and tramway are set perpendicular to the magnetic meridian. 
and to the axes of the Parry coils. The demagnetizing solenoid (i.d. 
14-7cm, o.d. 20-8em, length 12-5em)-contains 1920 turns in 20 layers of 
18s.w.g. enamelled copper wire and is similar to that described by Parry 
(1957). The resistance is 15-6Q and inductance 0-5H. The coil has 
been calibrated by a.c. and d.c. methods and both methods agree to 1% 
and give a calibration of 105 oersteds/amp. Power is provided by public 
supply mains (240v; 50c/s) through a bank of condensers (20 nF), which, 
with the coil, form a circuit tuned to 50 cycles. This reduces undesirable 
harmonics and as an extra precaution a series-resonant circuit tuned to 
100 cycles has been connected in parallel. The required voltage is obtained 
from a variable transformer. The rotations about the vertical and hori- 
zontal axes are in the ratio 15:16. Gears provide primary rotation rates 
about the vertical axis of 140, 350, 525 and 700rev/min, which do not 
resonate with any simple fraction of the mains frequency. It was found 
experimentally that the results are unaffected by changes within this range 
and the lowest rate has been used in these experiments. 

The direction and intensity of magnetization of each specimen is first 
measured using an astatic magnetometer. Directions are determined to 
3° and intensity to 5%. The specimen is treated in a small alternating 
field and then remeasured. Treatment in a higher field is then given and 
further measurements taken, and so on. The peak values (H,,~) of the 
alternating fields used are always quoted. The field strengths generally 
used are 75, 150, 300, 450, 600 and 750 oersteds, although in many cases 
intermediate steps are included. Provided the field used is less than that 
required to reduce the specimen to its minimum intensity (a situation which _ 
is explained later) repeat treatments give the same results. 


§ 3. Rock SPECIMENS STUDIED 


The rock specimens chosen for study were collected from eight sites 
(fig. 1) in the Tertiary Basalts of New South Wales (see general account of 
geology in David (1950)). The specimens were chosen because they showed 
different degrees of magnetic stability. The NRM consists of a primary 
component together with secondary components of varying magnitude 
directed approximately along the present geomagnetic field or along the 
field during storage in the laboratory. 

Natural exposures in these basalts are always extensively weathered 
and the collections have been made in fresh exposures in quarries and 
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Fig. 1 


SYDNEY 


f 
«! 


MELBOURNE TASMAN 


SEA 


Older Votcanics, 
Victoria 


Locality sketch-map. The basalt outcrop is in solid black. The collecting 


localities are numbered | to 8 as in the text and tables. The Older 
Volcanics of Victoria which has been studied previously are also marked. 
I. Road cutting 2km N of Nimmitabel on Monaro Highway. 


2. Road cutting on Hume Highway, 5 km NE of Berrima (Site B, 
Green and Irving 1958), (Map reference : Mittagong 370437) 
3. Three sites, one in Robertson village (Kiama 627294), one on 


Robertson-Moss Vale Road, 3km W of Robertson (Kiama 560284) 
and the third on Belmore Falls Road (Kiama 594275). 


4. Road cutting, Nundle-Timor Road, 18 km SW of Nundle. 
5. Road cutting on Oxley Highway, 53 km E of Walcha. 
6. Road cutting on the Guyra Road, 5 km S of Ebor, 


7. Two quarries, one on site of University of New England and the 
second 24 km NW of Armidale. 


8. Road cutting 50m W of bridge over Myall Creek on Delungra- 
Bingara Road, 24 km E of Bingara 
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road-cuttings spread over a meridional distance of 650km. The basalts 
are flat-lying and no geological correction to the directions of magnetization 
is necessary. In six cases, each site is a quarry or road-cutting in which 
the samples cover 5-10m in a single flow. At Armidale (Site 7) samples 
were obtained from two quarries about 1km apart. At 3, three road- 
cuttings 2-3km apart west of the village of Robertson were sampled. 
At both these sites more than one flow may be represented, although if 
this is so they cannot be much separated in time. Several separately 
orientated samples have been taken at each site and from each sample 
one or more cylindrical specimens (3-5cm diameter, 3-5cm height) have 
beencut. Thedetailsaregivenintable1. Altogether there are 37 samples 
and 66 specimens. 


§ 4, RESULTS 


The demagnetization procedure has been carried out on all specimens, 
about 500 determinations of the direction and intensity of magnetization 
being necessary. The effect on the directions of magnetization in specimens . 
from the same site is shown for two sites in figs. 2 and 3. The variation in 
the within-site dispersion may be illustrated by plotting £/ksig against the 


Fig. 2 
TN 
Effect of alternating fields on the directions of magnetization at site 2. 
(a) NRM, and after treatment in (6) 150, (c) 300, and (d) 450 oersteds. 
The projections used in this and figs. 3, 7 and 8 are polar stereographic. 
North-seeking directions are plotted as dots on the lower hemisphere 


i en on the upper hemisphere. The inclination scale is given in 
a). 
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Fig. 3 


(b) 


(c) 


Effect of alternating fields on the directions of magnetization at site 5. 
(a) NRM, and after treatment in (b) 75, (c) 150 and (d) 750 oersteds. In 
(a) P is the dipole field direction and F the present field both on the 
upper hemisphere. The great circle (continuous on the lower hemisphere 
and dotted on the upper) through P and the mean direction after 
treatment in 75 oersteds is given in (a) and shows the ‘ stringing’ of 


Mp. 


Fig. 4 


400 
Oersted 


Effect of alternating magnetic fields on the within-site dispersion. This is 
illustrated by plotting the parameter k/hg (see text), which is 
proportional to the inverse of the dispersion, against alternating field. 
The horizontal axis’ is drawn at k/k.g=1 and values below this indicate 


a random within-site dispersion. 


E. Irving et al. on 


bo 
ww) 
bo 


Fig. 5 


Oo 
So 


ial ee Lae ae 
400 600 800 
Oersted Hp~ 


=a 
200 
Demagnetization curves for specimens from sites 2, 3, 5 and 8. The value 


M,, is plotted against alternating field H,~ and M,, is the intensity 
after treatment, normalized to M,.=1 at H)~ =0. 


Fig. 6 


1-0 


“ac O.5 


200 400 600 800 He 


Oersted 


Demagnetization of IRM and TRM. The value of J/,, is plotted as in fig. 5. 
The initial intensities of TRM in 0-6 oersted (Mto.g)) are 0:54 x 10-3 
e.m.u./cem’ for a specimen from 8 and 22:5 x 10-3 for a specimen from 3. 
The IRM (Migoo)) acquired by another specimen from 2 in a field 
of 300 oersteds is 108 x 10-°. 
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strength of applied alternating field. This has been done for all sites and 
representative curves for four sites are given in fig. 4. The parameter k 
is an estimate of « the dispersion parameter occurring in the direction 
distribution exp («cos@) (Fisher 1953), and ksig=(N — 1)/(N —R,) where 
f, is the value, tabulated by Watson (1956), of the resultant of V random 
vectors which will be exceeded in only 5% of cases. When k/ksig< 1 
the distribution may be considered random. 

_ The associated variations in intensity are illustrated by results from 
specimens from the same four sites (fig. 5). For comparison the demagne- 
tization of IRM and TRM are given in fig. 6. The specimens were first 
demagnetized in 900oersteds alternating field and then magnetized in 
300 oersteds steady field to give M,/..)), or heated in nitrogen to 625°C 


and cooled in the earth’s field (0-6 oersted) to induce M,,o.¢). 


Table 2. Some magnetic properties. Wy is the intensity of the natural 
remanent magnetization Mmin is the minimum intensity (see 
§ 7), Wisat) is the saturation isothermal remanence, Haz is the field 
required to saturate and Hey is the coercivity of remanence of 
Misat) obtained by the methods described by Creer et al. (1959). 
Intensities are given in units of 10-°e.m.u./em® and field strengths 
in oersteds 


Raah ae M, Min Meat oe ie 

| 10-48 0-06 oe 101-8 900 200 

2 5:95 0:82 188-0 1700 225 

3 32°25, 1:34 224-5 1000 205 

BS 4 13-74 0-26 354-0 1500 190 
5 2°52, 0:24 er 354-0 900 180 
fea 6 1-25 0-11 405-0 1200 180 
2 7 7-13 0:86 236-5 900 160 
8 0-014 0:005 0:32 4000 600. 


Certain other magnetic properties are listed in table 2; amongst these 
it may be noted that the ‘back field’ (the coercivity of remanence, A es) 
required to remove the saturation isothermal remanent magnetization 
Miat)) is usually about 200 oersteds. ; 
These results are now discussed dealing first with the effects in low 


alternating fields. 
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§ 5. Errects in Low ALTERNATING Frevps (H, ~ < H,,) 


The results within each site which are relevant to the first condition in 
§ 1 are considered first. At site 8 the directions of NRM (Mn) have large 
positive inclinations. The directions are well grouped and this is not 
improved by treatment in alternating fields. In fields up to 400 oersteds 
the demagnetization curve of My is convex upwards very similar to 
Mo.) induced in the laboratory. The small initial rise is presumably 
due to the removal of a soft component, along the earth’s field and which is 
approximately anti-parallel to the hard component. The coercivity of 
remanence (Her) is high (table 2). It may be concluded that the NRM is 
predominantly stable and of TRM type. 

At site 3 the directions My are random (see fig. 5 (b), Boesen et al. 1961), 
but these become grouped together with steep positive inclinations after 
treatment in alternating magnetic fields of 300oersteds (table 1). The 
demagnetization curve is convex downwards similar to that for IRM 
induced in the laboratory although the initial decrease is not so rapid. 
The directions in specimens after storage in the laboratory for three 
months can change by as much as 180°. The value of Her is low (table 2). 
We conclude that the NRM of this basalt consists of a highly unstable 
magnetization and a stable component of much smaller magnitude with 
steep positive inclination. 

At site 2 the directions Mn are well grouped with steep positive 
inclinations (fig. 2(a)). At 5 there is a group of specimens with similar 
directions but also there are others with low inclinations and north-westerly 
declinations strung out close to a great circle through the dipole field (fig. 
3(a)). At both sites after treatment in small alternating fields the disper- 
ion decreases and the directions all become southerly with large positive 
inclinations. The minimum dispersion occurs after treatment in 100—200 
oersteds. The demagnetization curves lie mainly between those for 8 
and 3 and between M,,,.4) and M/s) induced in the laboratory. The values 
of Heyarelow. The conclusion is that the NRM consists of a comparatively 
stable component of TRM-type and an unstable IRM-type component of 
similar or lesser magnitude. The stable component is reversed and is 
appreciably affected by low demagnetization fields since the demagnetiza- 
tion curves do not show the initial rise as at 8. The behaviour of specimens 
at the other sites (1, 4 and 7) is similar to that at 2 or 5. 

Turning now to the comparison of results between-sites (condition (2) 
in §1) it is immediately clear that the effect of treatment in low fields 
(300 oersteds for site 2 and 150 oersteds for the remaining sites) is to bring 
the widely scattered mean directions of NRM at each site into close agree- 
ment (table 1 and fig. 7(a) and (b)). The divergences remaining between 
the site mean directions are due mainly to experimental error and any 
residue may be satisfactorily explained as due to the secular variation in 
the geomagnetic field during the time spanned by the lavas. 

The third condition may now be tested by comparing the directions with 
those previously obtained from Victoria (see fig. 1). The directions 
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irrespective of sign are very similar to those observed in the Older Voleanics 
of Victoria (compare fig. 7 (b) with fig. 6 of Irving and Green (1957 )). The 
poles for each site are given in table 1 and the mean of these (63S, 137 E) 


Rigs 7 
Site mean directions. (a) NRM, (6) after treatment in low alternating fields. 


The initial within-site ee at 3 are random and there is no 
significant mean, but after treatment in 300 oersteds alternating field 
the directions become grouped together and the mean of these is included 


in (6). 
Fig. 8 
Oo 
Older Volcanics 
Victoria 
9OW Mees 9OE 


Tertiary Basal 
N.S.W. 


Pole position for the N.S.W. basalts. The pole for the Older Volcanics of 
Victoria is added for comparison. The ovals around the poles mark 
out the zones within which the pole lay at a probability of 95%. The 
south geographic pole is marked with a cross. 
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agrees well with the pole for the Older Volcanics of Victoria (67 S, 123 E) 
(see fig. 8). The spread of time of the New South Wales basalts is probably 
considerable, but does not necessarily overlap that of the Older Basalts of 
Victoria although both are Lower Tertiary. The Lower Tertiary sequence 
of Victoria has in some places been gently folded but the lavas have been 
sampled only at those places where the dip of adjacent sediments is 5° 
or less (Irving and Green 1957). Thus here is an example of agreement 
between pole positions obtained from rocks from the same land unit which 
have suffered no appreciable tectonic movement and which are of similar 
geological age. 

It is evident that after treatment in low alternating fields the three 
conditions listed in § 1 are satisfied and that the effects of the secondary 
magnetizations imposed by the present earth’s field have been removed. 
Since the stable magnetization, for instance, at site 8 where the basalt 
has not been appreciably affected by a secondary magnetization, is of 
TRM type, and in the absence of any other known geophysical force 
adequate to produce consistent directions of magnetization over a meri- 
dional distance of a 1000 km we may assume that the directions of magne- 
tization remaining after this treatment are the directions of the primary 
TRM and thus of the geomagnetic field at the time the lavas cooled. 

One further point. which for the sake of completeness needs to be 
considered here is the effect of burial stress, which could conceivably have 
had a systematic effect on the inclination. Except for laterite and soil 
cover, no deposits overlie the basalts. The thickness of the lavas is not 
usually more than 100m. The maximum is probably to be found in the 
Liverpool Ranges (fig. 1) where 500m of lavas occur which may be taken 
to be not more than a factor of two less than the maximum burial depth. 
To test this point the magnetization of pairs of dises cut from two samples 
at each of sites 5 and 7 was further examined in the laboratory. In one 
of each pair, IRM in a field of 100 oersteds, M;,,5)) and TRM in the earth’s 
field, M,,o.¢) was induced whilst the specimens were subjected to an axial 
load of 500kg/cm?, and in the other, the same magnetizations were pro- 
duced without load. The apparatus and method have been described by 
Stott and Stacey (1960). There was no measurable deflection of the 
magnetization due to the application of stress. The directions of magne- 
tization in each pair agreed to within the measurement error. The applied 
stress corresponds to a depth of burial of about 1500m and it is unlikely 
therefore that any systematic effect due to stress during burial is present. 
No self-reversal of TRM was found, with or without stress. 


§6. PROCEDURE FOR REMOVAL OF SECONDARY MAGNETIZATION 


We suggest that the following procedure may be applied to volcanic 
rocks. ; 


1. Ata particular site measure Hy for one specimen; there is not usually 
much variation between specimens from the same site. 
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2. Measure the directions of NRM of 5 (say) specimens preferably from 
different samples and the directions after treating in alternating fields 
(H, ~) of 4H or, 4Her, Her and 2H oy. Find the field at which the points 


are least scattered, either by plotting & against H,~ or by plotting on a 
stereogram. 


3. Treat all specimens in this alternating field. In those cases where 
the primary and secondary components are of the same magnitude (that 
is, when the directions Mn are distributed along a great circle between 
the present field and the directions of the primary magnetization) the 
strength of alternating field required is usually about 2H¢r. When the 
primary component is small compared with the secondary ones alternating 
fields as high as 14H cy may be necessary. 


§ 7. TREATMENT IN Hicu ALTERNATING Frexps (H, ~ > H,,) 


After treatment in alternating fields of 600 oersteds, which in most 
cases is several times greater than Her, the directions at each site 
approach a random distribution. The intensity does not diminish 
regularly with increasing alternating field, but levels off to an almost 
constant average value at about that field at which the dispersion begins 
to increase. There are however, fluctuations about this average. This 
value which is called the minimum intensity (Mmin), is now discussed. 

The effects may be due to instrumental imperfections or to inherent 
properties of the rocks. As and Zijderveld (1958) found that ‘ during 
the experiments of demagnetization by powerful alternating fields a 
disturbing magnetization appears to be induced in the samples”. Creer 
(1959), however, finds that spinning “prevents the specimens acquiring 
an IRM due to DC components’, but even so he still finds “stray 
magnetizations acquired during ‘cleaning’’’. We suggest that while 
the effects due to instrumental imperfections can be important if pre- 
cautions are not taken, there will always be some residual, fluctuating 
magnetization, even when the demagnetization takes place under perfect 
conditions. The reason for this is that spontaneous magnetization of 
domains is lost only at temperatures above the Curie point, so that the 
effect of alternating magnetic fields is one of rearrangement of domains, 
and because of statistical fluctuations the resultant moment of a specimen 
is never likely to be zero. The observed minimum intensity, and the 
scattered, almost random directions in alternating fields whose strength 
is much greater than the coercivity of remanence may therefore be a 
result of such fluctuations and will depend on the number of domains in 
a specimen. An estimate of this number may be made as follows. 

Assume firstly that there are n, ferromagnetic grains per cin® of equal 
size and moment M, whose directions are random in an isotropic 
demagnetized rock. Rayleigh (1919) has given the probabibty dis- 
tribution of the magnitude of the resultant 7) of n randomly oriented 
unit vectors for large values of n. The value of ry is given by ry>=Cn'? 
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where C is a constant having the value 1-61 for P(r >79) =0-05 and 0-89 for 
P(r>r,)=0-:50. The latter will give an average value for the mean 
random magnetization so that the random moment (Mr) per cm* is given 
by My=M,, .0-89.n,1?. Assume further that in each ferromagnetic grain 
there are n, domains of equal moment m so that the total number of 
domains per cc is n,.”,. We may say that domains have randomly 
arranged directions in a single demagnetized grain and the average 
moment of each grain will thus be M,=m.0:89.n2"?. However it is 
more likely that the domains behave as an assembly of anti-parallel 
vectors having the values +1 and —1 with equal probability. Rayleigh 
(1880) has given the probability for the length of the resultant in this 
situation also and we obtain r)= An? where A =0-83 for P(r >19) =0-50. 
The average moment per demagnetized grain in this case will be 
M,=m.0-83.n,'. The correct distribution probably lies somewhere 
between the random and anti-parallel situations but the difference for 
these purposes is negligible. Taking an average for the constant as 0-86 
the random moment per cc becomes: 


HM >in 016 nye), ee 
Now the saturation moment Ms=n,.n,..m and the ratio of this to the 
random moment will have the average value 

Mel Maxl:31 (njstar 
Taking the observed values of Migat) and Myin as estimates of MW, and 
M, respectively we find for a specimen from site 5 that Ms/M,;=1-5 x 10 
which give 1:3x10° domains per cm*. But Migat)<$Ms and also 
Mwmin> My, since the directions after the minimum intensity have been 
reached are not completely random due to the persistence of a very hard 
component even after treatment in high alternating fields. This will, 
therefore, be a minimum estimate. A better estimate may be obtained 
by putting Ms=2Migat) and M;y=4M in) which gives M,/M,=5-9 x 108 
and 20-4 x 10° domains per cm3. 

Examination of a thin section of this specimen shows that although 
there is wide range in the size of visible magnetic grains, about 80% of 
the magnetic material is contained in grains of size 25-150 and over 
99% in the range 25-175. A modal analysis of this section shows 
that the magnetic minerals oceupy more than 10% by volume of the 
rock, Assuming that the material of grain size > 25. is responsible 
for the observed minimum intensity and that these grains are spheres of 
average diameter 100,, there will be 1:9x105perem?. From the 
estimates of the total number of domains, this gives the order of 10 or 
100 domains per particle and suggests that the minimum intensity is 
due to the presence of multi-domain particles. For the reasons just 
mentioned the estimate of ten domains per particle will be a minimum. 
The estimate of the order of 100 is likely to be more correct and there is 
no reason to regard this as a maximum estimate. 

Since this paper was written 30 specimens from these basalts have been 
thermally ‘demagnetized ’, by heating to 610°C (more than 50° higher 
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than the Curie Points of the iron minerals they contain) in nitrogen in a 
non-magnetic oven followed by cooling in zero field. After this treatment 
the specimens still retained an appreciable moment. The values, which 
are to be compared with the minimum intensities obtained by demagnetiza- 
tion in alternating magnetic fields (table 2, column 3) are these: site 
1, 0-11; 2, 0-64; 3, 1-67; 4, 0-22; 5, 0-23; 6, 0-16; 7, 0-43; 8, 0-01. The 
values are an average of about four specimens per site and are given in 
units of 10e.m.u./em®. As with demagnetization by alternating fields 
these small residual moments could result from instrumental imperfections 
due, for instance, to incomplete cancellation of the earth’s field. The 
de-gaussing coils have a central space of about 8000 cm’ which is uniform 
to about I part in 10° and the oven, which has a working space rather 
less than this, holds about 20 specimens spaced sufficiently far apart to 
avoid mutual interference from their moments. The directions in each 
batch after treatment are random so that it appears that the residual 
field in the oven has negligible aligning force on the specimens. The 
important point in the present context is that the magnitudes of the 
moments, /,,,,, are the same whether they are obtained by thermal or 
alternating field demagnetization ; it would be a remarkable coincidence 
if instrumental imperfections in two different methods produce such 
closely comparable effects. These results therefore provide strong 
evidence confirming the view that the minimum intensity is an inherent 
property of these basalts and results from the statistical nature of the 
remanent magnetization. 

The accuracy of the procedure of removal of the unwanted secondary 
magnetizations by alternating magnetic fields will be limited by the 
presence of the minimum intensity, in those cases where the magnitude 
of the soft component is large, since the minimum intensity of this may 
exceed the magnitude of the hard component. However, since the 
minimum intensity directions are random, the average of several speci- 
mens will give an estimate of the direction of the primary component 
although the individual specimen directions will be widely scattered. ‘This 
is probably the case at site 3 where even after removal of the secondary 
component the directions are widely scattered but the mean is significant 
and agrees with the direction of the primary component obtained at the 
other sites. In rocks containing multi-domain particles which have become 
demagnetized over long periods of time a minimum intensity will always 
be present and may be the cause of the random scatter in the directions 
of NRM observed in specimens taken from the same sample of lava and 
which in our experience is usually greater than the experimental error. 
The presence of small temporary components could also explain this. 


§ 8. FERROMAGNETIC COMPONENTS 


In addition to the grain sizes notes above there are a large number of 
very small grains of the order of a micron in diameter and a small portion. 
‘of the stable NRM may be due tio these. The directions of magnetization. 
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do not become completely random (figs. 2 (d) and 3 (d)) even after treatment 
in high alternating magnetic fields suggesting that there is an exceedingly 
hard component present and these grains, acting as single domains, 
may be responsible for this. They could not be responsible for all the 
stable NRM since they represent less than 0-1% (by volume) of the total 
ferromagnetic content. It seems therefore that both the primary and 
secondary magnetizations reside in grains which, because of their size 
and the considerations in the previous sections, are multi-domain. The 
evidence presented in this paper is consistent with the view that in the 
natural rock these grains are for the most part either demagnetized 
and contribute only their minimum intensity, or possess a secondary 
magnetization, but that they contain within themselves regions of much 
higher magnetic stability responsible for the primary NRM which, as 
was shown in §5, is probably thermo-remanent in origin. It may be 
noted in fig. 6 that the IRM induced in a steady field of 300 oersteds 
(M,-399)) is removed by alternating magnetic fields of the same strength. 
Moreover, if Mi 499, is induced in a specimen previously possessing a 
TRM (M,.0.6)) which has been treated in alternating fields in excess of 
300 oersteds (say 500 oersteds), the underlying TRM remains unaffected 
after the removal of Miso. The two components (Mis), and Myo.) 
after treatment in 500oersteds H,~) are independent of one another 
and are probably related to two different ferromagnetic components in 
these basalts, the former perhaps to the multi-domain ferromagnetic 
grains and the latter to high stability regions within these. 

An estimate of the fractional content of grains represented by high 
stability regions may be obtained from considering the demagnetization 
of Mio.) (fig. 6). The softer component which makes up about 90% 
of the intensity is largely removed in 500 oersteds, and it is this component 
which is likely to become demagnetized in the natural rock over a long 
period of time, or which has acquired a secondary magnetization. In 
fields above this, value a hard component, which is 6% of M o.¢. in the 
highly unstable Robertson basalt (site 3) and 15% of M yo.g) in the stable 
basalt at site 8, still remains. It seems that these high stability regions 
occupy approximately 10% of the ferromagnetic grains. 

Verhoogen (1959) has recently suggested in general terms that such 
high stability regions within ferromagnetic grains are due to structural 
dislovations. Another possibility is that intergrowths are responsible ; 
the fine scale intergrowths contributing the hard magnetic component 
and the coarser-grained matrix providing the softer component. 
Dislocations are more likely to occur in rapidly cooled rocks whereas 
slow cooling will favour the development of suitable intergrowths. It is 
noteworthy that in extensive work on hypabyssal dolerite intrusives 
(Almond et al. 1956, Irving 1956, Graham and Hales 1957, Bull and 
Irving 1960, Creer et al. 1959) and plutonic intrusions of gabbro 
(Gough and Niekerk 1959, Khan 1960) few cases of instability have been 
found, whereas instability in basalts is very common (Brynjélfsson 1957 
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Creer 1958, Khan 1960, this paper). In dolerite and gabbro the ferro- 
magnetic grains are larger than those in basalt and it is clear that there 
is no simple correlation between stability and the overall size of the iron 
mineral grains. It seems more likely that the question of their internal 
subdivision is of more importance. It seems not unlikely that the 
explanation for the hard and soft magnetic components which are found 
not only in these basalts discussed in this paper but also in most other 
igneous rocks studied is to be sought in the exsolved phases of differing 
grain size in ferromagnetic grains. 
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‘Conservative Climb’ of a Dislocation Loop due to its 
Interaction with an Edge Dislocation} 
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ABSTRACT 


A new type of dislocation motion, called ‘ conservative climb >, has been: 
observed in zine platelets deformed in tension inside an electron microscope. 
Circular, prismatic dislocation loops are often forced to climb out of their 
glide cylinders by their repulsive interaction with edge dislocations. It is 
observed that the area inside a loop is conserved during its movement, from 
which it is concluded that the climb occurs, not by self-diffusion, but by 
the transfer of vacancies around the loop by pipe diffusion. Calculations 
of the interaction energy, of the total forces on the loop, and of the forces 
per unit length on the elements of the loop show that conservative climb 
can be expected to occur when the temperature is not too low and the 
separation of loop and edge dislocation is small. 


DURING a transmission electron microscope study of the formation and 
climb of dislocation loops in deformed zine crystal platelets (Price 1960), 
an unusual type of dislocation motion was observed, which has not pre- 
viously been considered. A number of crystals were deformed on a 
tensile straining stage (Fisher 1959) inside the electron microscope, and it 
was often observed that during pyramidal glide dislocation loops were 
generated by moving screw dislocations with a Burgers vector in the 
direction c+a=1[1123]. Figure 1(q@)|| is a transmission micrograph 
showing loops with a [0001] Burgers vector which are lying in the basal 
plane in a row along the path taken by the screw dislocation (not shown 
in the figure) from which they were formed. The basal plane is parallel 
to the plane of the micrograph. Also in the micrograph is an edge dis- 
location with a Burgers vector a=}[1120]. (For the method of deter- 
mining the Burgers vectors, see Price (1960).) One end of the dislocation, 
which is situated mainly on the basal plane, ends at the crystal surface at 
A. During straining it was observed that the dislocation, in attempting 
to glide past the loops, forced one of them to move slowly sideways as 
shown in the sequence of fig. 1(a), (b) and (c). In fig. 1 (c) it has also 
begun to displace another one of the loops. The area inside the loops 
appeared to be conserved during the motion, so that there was no net flow 


of vacancies to or from the loops. This type of loop movement as a result 
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of interaction with a moving dislocation was observed several times and 
will be called ‘conservative climb’ to distinguish it from ordinary climb, 
which involves vacancy diffusion. The purpose of this note is to give a 
qualitative explanation for conservative climb in terms of classical 
elasticity theory. 

Let us consider the interaction between a circular edge dislocation loop 
and an infinite, straight dislocation whose glide plane is parallel to the 
plane of the loop. We assume an infinite isotropic medium characterized 
by the shear modulus », and the Poisson constant v. The loop lies in the 
plane z=0 (fig. 2) with its centre at the origin of coordinates, with radius 
R and with Burgers vector b® along the z axis. The straight dislocation, 
1, is located at y,, z,, extends along the x axis and has a Burgers vector b’. 


Fig. 2 


Circular edge dislocation loop and infinite straight edge dislocation. 


If b! is along the wx axis (i.e. if the dislocation is pure screw), the inter- 
action energy of loop and line is zero and the total force on the loop is zero. 
We now consider the case of a pure edge dislocation with b! along the y 
axis (fig. 2). The components of the stress tensor, o;;’, for the dislocation 
I are given in the book by Cottrell (1953), with his coordinates 2, y, 2 
replaced by y-41,2—%,. The interaction energy H, was obtained by 
integrating the expression H°!=b°{/o,1(x, y, 0)dady over the area, A, of 
the loop. The result, = 
Hl — (yb%t/(1— v))RZ{1 —{2(V2+ 22-1) + 2y/[(1+ Y24+ 222-42] We 

x [E+ AP F)\(¥ = 1) Ze) + (Ya Zh YY Zana 
is plotted in fig. 3 as a function of Y, for Z=1 and Z=0-1. wh = 
and Z=z,/R. Reapers ce 

The total forces on the loop were found from the ex i 

pressions F’,°=0, 
P= 0b" /0y,, F,°=0H°'/dz,. F,°, the force tending to cause climb, is 
plotted in fig. 3 as a function of Y for Z=1 and Z=0-1. F,°, the force 


: 
Conservative Climb’ of a Dislocation Loop 245. 


tending to drive the loop out of the crystal by glide (not plotted in the 


figure), has a maximum value which is smaller than that of F ,° and occurs 
at a different value of Y. 


Fig. 3 
2 


dS 


W4 
Interaction energy, ¢®!, and total force, ® ,°, on the loop as a function of Y, with 
69 = 2 /[ub°b!R/2(1—v)], 0,°=F' ,°/[ub°b'/2(1—v)]. The full lines are 
for Z=1 and the broken lines for Z=0-1. 


The force per unit length, f°, on every element of the loop was obtained 
by applying the Peach—Koehler formula (Peach and Koehler 1950). At 
the point (Rcos6, Rsin 8, 0) the components of f° are 

epee \sind= 7) 2? | 


Oe eee ee Se  CONY, 
ver ay A(T Ee So cee 


b°D! Z[(sin6 — Y)?—Z2?] 


c= ee ee. BIND, eee 
lie aRaienee Veer e ?) 


b°D! (sin@— Y)[(sin@— Y)?—Z?] . 
S em ae 
im i. R [(sin fx Yy2+ 227 sin | 
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where D! = .b!/27(1—v). The radial forces, given by fr? = [(f,°)? + (fy)? Ph? 
depend only on o,,1, which is zero along a plane at 45° to the glide plane 
of dislocation 1 (see fig. 2) and changes sign across this plane. For example, 
when Y=Z=1, that is, when the dislocation is a distance & above one 
edge of the loop, the plane intersects the loop at y=0 and the forces are 
directed toward the centre on one side and away from the centre on the 
other side, as shown in fig. 4. The lengths of the arrows in the drawing 
are proportional to the magnitude of the forces. 


Fig. 4 


Radial forces, ¢r°, on the elements of the loop for Y=Z=1, with 
gr°=fr°/[ub°b*/27(1 — v) RP}. 


We now apply these results to the electron microscope observations, 
assuming that they are qualitatively the same when the dislocation 1 is 
finite and the medium is crystalline, anisotropic, finite and bounded by a 
thin oxide film. The Burgers vectors b° and b! are directed along the ¢ 
axis and an a axisand have the values 4-94 A and 2-66 4 respectively. Thez, 
y plane is the basal plane of the platelet. If the glide plane of the disloca- 
tion is far from the plane of the loop (Z large), the total force on 1, — ae 
may be small enough for it to glide past the loop. If Z is not large, 
the dislocation will be unable to surmount the energy barrier and will be 
held up near the loop and exert a repulsive force on it. 

Since /’,°< F',°, it is not surprising that prismatic glide of the loop was 
not the predominating process, especially when the Peierls force and the 
repulsive force of the oxide layer at small distances are considered. The 
total forces on the loop favour conservative climb, and if the temperature 
is not too low it should be quite easy for this to occur at a reasonable rate 
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by the transfer of vacancies from one end of the loop to the other by pipe 
diffusion along the core of the dislocation. Since the Burgers vector of 
the loop is quite large—almost 5 &—the activation energy for this process 
should be much smaller than for the movement of vacancies through a 
perfect lattice. 

The large total force F',° and the conservative climb of the loop are due 
to the large gradient of the stress component o,,1._ (The average gradient 
of o,, in the region of the loop is ~105kg/mm?/mm for Z=1 and 
~10%kg/mm?/mm for Z=0-1.) Although we have made observations 
and calculations for only one special case of conservative climb, there is 
no reason to doubt that this type of dislocation motion is a general one as 
a result of the interaction of closed prismatic loops with each other and 
with other dislocations, provided a suitable stress gradient is available. 
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ABSTRACT 


The concept of the stability of a wavefunction under a perturbation is 
defined and the properties which distinguish stable wavefunctions are deduced. 
The stability of approximate wavefunctions is particularly examined and a 
number of general results stated. Self-consistent wavefunctions are shown 
to be stable under one-electron perturbations. The concept enables many 
hitherto isolated results to be unified and extended. In particular a new 
method of calculating nuclear magnetic shielding constants for atoms is 
given. The concept also suggests means of improving an approximate 
wavefunction so that the derived estimate of some property is made more 
accurate. 


§ 1. INTRODUCTION 


THE application of quantum mechanics to the study of atoms and 
molecules is made difficult by two related obstacles. The first is that 
the number of problems, for which exact solutions in terms of familiar 
functions can be found. is very severely limited and the second is that 
there are few theorems which can be used to simplify the evaluation of 
properties or to estimate the accuracy of approximations. Almost all 
calculations have to be made using approximate wavefunctions and the 
theoretical values of properties involve inaccuracies. It is very difficult 
to estimate these inaccuracies theoretically and their behaviour, judged 
from agreement with experiment, is often erratic. There have been 
several examples published in which the accuracy of the wavefunctions, 
as judged by the minimum energy criterion, bears no relation to the 
accuracy of the estimate of some derived property. Sometimes a change 
in a wavefunction, which produces no appreciable improvement in the 
energy, may make a large difference to the estimate of a property. 

In this paper the concept of the stability of a wavefunction under a 
perturbation is introduced and used to illuminate some of these difficulties. 
Much of the erratic behaviour is associated with the use of instable 
wavefunctions. A number of different methods of making wavefunctions 


stable is also discussed. 


§ 2. DEFINITION OF STABILITY 


The estimation of many atomic and molecular properties can be 
considered in terms of perturbation theory. The unperturbed equation 


+ Communicated by the Author. 
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is usually taken as the non-relativistic Schrédinger equation and the 
perturbation may be a relativistic term, some additional internal operator, 
a change in one of the terms already present or a term representing the 
interaction of the system with some external field. The perturbed 
Schrédinger equation is written 
(HEAP) b= Wb. 3. oe eee tha 
where A is a parameter governing the strength of the perturbation P, 
and the unperturbed equation has eigenfunctions ¢,, satisfying 
Hd.= Wd)... = OL a Oe 
In most applications only the ground-state solution of (2.01) is needed. 
This equation can be generalized by considering a Hamiltonian 
Ha), a8 (2.03) 
whose eigenfunctions are known for some non-zero Ay. This generalization 
is easily carried out and need not be given here. The quantities of interest 
are usually the coefficients of A and A* in W(A). It will be assumed that 
the operators H and P are of such a kind that the energy levels W(A) 
are differentiable functions of A at least atA=0. A discussion of conditions 
sufficient to ensure this has been given by Kato (1951). 
If the eigenfunction (A) in (2.01) is normalized, so that 


[vara tie (ee, Amara 


for all A, then the energy can be expressed as 
WO)= [YUL +AP dr. rete rs 
In general, therefore, since % depends on A, 
Ce ira adifs* ety ob 
+= | Part [Fi Hbar+ | Ho dr, . . (2.08) 


but for certain wavefunctions 


dlaj* d 

Te Hides | yrHh dr=0, cd OU ROD) 
and so 

dW 

x = | v"Plde= (P). aa 


Such wavefunctions will be called stable. More generally, if % is only 
an approximate wavefunction, (2.05) is the corresponding approximate 
energy and this % also will be called stable at A, if it satisfies (2.07) for 


some A=A,). If % is not normalized the corresponding equations to 
(2.07) and (2.08) are 


d * 
Te Wynd + [rit W) War =o, (2.09) 


y= [etPbar| fey drs (P). = ao eme (200) 
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§ 3. Tue Stapiiry or Exact EIGENFUNCTIONS 
The discussion of stability is easier when % is an exact solution of 
(2.01) but there are three different cases to be considered. Tf, inthe 
first place, the lowest eigenvalue W, is non-degenerate then % is stable 
for any perturbation which gives rise to an energy differentiable at 
A=0. This can be proved from the definition by using (2.02), (2.04) 
and the limiting form of % at A=0, i.e. do, for 


dif yy ee difs* y Us 
la Heyie+ {de AF dr=Wo| bode + Wo| Jar 


ON ya. 
=Wo | prbdr 
=!0; OM eth Ghee. pee ER 
where the derivatives are taken at A=0. 

When W, is degenerate and the degeneracy persists in the presence of 
the perturbation then, even though the energy varies smoothly, the 
individual wavefunctions can be discontinuous since the choice of basis 
functions in the degenerate sub-space may vary from point to point. 
To overcome this difficulty it is usual to introduce the operator or 
operators necessary to label the degenerate wavefunctions uniquely and 
to consider only those wavefunctions in which these labels are kept 
constant. Thus the simultaneous eigenfunctions of the energy and these 
additional operators are uniquely defined and so, according to the first 
part, can be differentiated with respect to A and will satisfy (3.01). 

When W, is degenerate and the degeneracy is split by the perturbation 
then the wavefunctions can be discontinuous but only at A=0. Thus 
the derivative dys/dA exists except possibly at A=0 which is the point 
involved in (3.01). If the degeneracy at A=0 is resolved by defining % 
there by 


HD) stl), 2p Cannes fen ee O42) 
then this difficulty is avoided. Thus the perturbation is regarded as 
virtually present even whenA=90. The usual modification of perturbation 
theory consists in choosing the zero-order degenerate wavefunctions so 
that the matrix elements of the perturbation are diagonalized. These 
are already the functions necessary for (3.02) and consequently are the 
stable wavefunctions. If the splitting is only partial it will be necessary 
to combine the techniques above. The significant feature of this case is 
that the stable zero-order wavefunctions depend on the perturbation 
and, in general, will differ according to the perturbation. When two 
perturbations are applied simultaneously the stable zero-order wave- 
functions may even vary as the magnitudes of the perturbations vary. 

The calculation of the dipole moment of the hydrogen atom in its 
first excited state gives an illustration of the influence of the degeneracy. 
In the non-relativistic theory the first excited state has 4-fold degeneracy 
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with wavefunctions which can be labelled 2s, 2p,, 2p,, 2p, Each of 
these functions gives a symmetrical distribution of charge and so has no 
dipole moment. This does not imply, however, that the atom has no 
dipole moment since these are not the stable zero-order wavefunctions 
in the presence of an electric field and an electric field is required for a 
dipole moment to be observed. The electric field H in the z direction 
produces an additional term in the Hamiltonian so that the Schrédinger 
equation becomes 


(H—E2)b= Whee yn ee a) 
and the dipole moment is defined as 
ow 
= —{ ~— aaMecage ior 20) ok Mae e 
be (Sr) E=0 Z ( ) 


The stable wavefunctions are, then, those linear combinations of 2s; 
2p,, 2p,, 2p,, which diagonalize z and these are the digonal hybrids 


JEEP) Vee ee oe 


and the degenerate pair p, and p,. For these stable wavefunctions 
equation (2.08) leads to the simpler formula 


= eee 


The hybrid functions (3.05) no longer have symmetrical distributions of 
charge but have dipole moments of +7-6 Debyes. (For further details 
see Van Vleck (1932), p. 204.) 


§ 4. THE StaBiyiry or APPROXIMATE WAVEFUNCTIONS 


The definition of the stability of a wavefunction under a perturbation 
by means of eqn. (2.07) applies equally well to approximate wavefunctions. 
Indeed its importance is even greater since it gives a new method of 
discussing the probable accuracy of the properties deduced from an 
approximate wavefunction. 

There are a small number of examples in which it is possible to find 
approximate wavefunctions which are completely stable in the sense of 
satisfying (2.07) exactly. Thus both the exact and the approximate 
wavefunctions may be eigenfunctions of the perturbing operator for 
all A and consequently will satisfy 


Op 

ais 0. 
which ensures (2.07). Practical examples of this include perturbations 
depending only on angular momentum operators. 

In general, for most approximate wavefunctions and most perturbations, 
only approximate stability can be attained even when there is no 
degeneracy present. The basic method of finding wavefunctions which 
are stable to a good approximation is to introduce into the approximate 


(4.01) 


Stability of a Wavefunction under a Perturbation 253 


wavefunction parameters which enable it to adjust to the presence of 
the perturbation. Thus if the approximate wavefunction is 


PA)=P(oy, Mar s--> Mn)» - - » + . (4,02) 
where the parameters 4; depend on A, and if these parameters are 
determined by minimizing the unperturbed energy so that they satisfy 


a * fad 
[Zo auars oe hen) ea (4.03) 
Op; On; 
then by using this and 
djs Op dy; 
Dee (1704) 


eqn. (2.07) is satisfied. The degree of stability will depend on the 
success with which these parameters 1; can be used to represent the change 
in % due to the perturbation even though it is the wnperturbed energy 
which is minimized. The method may be described as determining the 
wavefunction by treating the perturbation as latent and the resulting 
stable wavefunction has the significance that it enables the effects due 
to the perturbation and those due to the approximations to be isolated 
at least to zero order in the wavefunction and first order in the energy. 
Similar effects in higher order terms have been discussed by Dalgarno 
and Lewis (1956) and Schwartz (1959). If there are degeneracies present 
the procedures above need to be combined and this can be done by 
minimizing the perturbed energy so that 
* 

SP (H+ Pybde+ [UML +AP) SE dr=0 ee HCAS) 
He i 
since this implies stability at any A and hence also at the limit A-+0. 
This expression should also be used for a non-degenerate wavefunction 
if there is a degeneracy induced by the perturbation, i.e. if the energy 
becomes degenerate for some particular 4,40. Sometimes the degeneracy 
is prevented by the non-crossing rule but still the form of the wavefunction 
changes very rapidly with A and (4.05) must be used. 

In one important particular instance of this basic method the para- 
meters are linear and the trial wavefunction is a linear combination with 
arbitrary coefficients of some fixed functions (the Ritz method). For 
many perturbations a suitable set of functions is the set suggested by 
Hassé (1930) and Lennard-Jones (1930), viz. 

ee EaG ae AE Wig enn eat a ae clot ema (2.00) 
where ;, is a good approximation to the unperturbed wavefunction. 
The use of this set is analogous to the power method of finding matrix 


eigenvectors. Corresponding to this another useful criterion for stability 


can be written as 
| (Piig)* Hy dr + [vor Py) demo. eee 07) 


Another linear set which is often useful consists of the functions 
Ui(Ao),  dyi{Ag)/da, MIN dn eae. B. .  w (4.08) 
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It may be more convenient to use this set than to permit A to vary freely 
particularly when the dependence on A is non-linear. With a sufficient 
number of terms the set is equivalent to a Taylor expansion in A and so 
leads to a stable wavefunction. 


§ 5. SIGNIFICANCE OF SELF-CONSISTENT WAVEFUNCTIONS 


A self-consistent wavefunction has the significance of being stable 
under certain perturbations. If the Hamiltonian for an atom or a 
molecule had the form of a sum of one-electron operators then, in the 
absence of degeneracy, a single determinant of orbitals would be an 
exact wavefunction. Such a wavefunction would also be stable under a 
perturbation which was the sum of one-electron operators. Since the 
true Hamiltonian also contains electron repulsion terms the single deter- 
minant is only an approximation. The best single determinant wave- 
function 


Y= Det{e (1) 2.) Oo.) Leer os 
where #,; are spin-orbitals, functions of both space and spin, is defined by 
the equation 


8[ ¥*(H — WF dr=0, 13° ecg sepaamareiay 


in which W is the approximate energy and § acts only on the various Wy. 
By the turnover rule or some equivalent procedure, this equation simplifies 
to 


5 | *(UMa"(2)... (HW) de=0 Re ks 


and, by treating the variation of individual functions, the set of 
equations 


[¥s80)... 80)... WYP dr=0, Sela ty (5.04) 


is obtained. The coefficient of the arbitrary variation 5, must then 
vanish and this condition is the usual differential equation for the self- 
consistent functions 4; Thus, for the self-consistent w;, eqn. (5.04) is 


satisfied for arbitrary 5%;. If the system is now subject to a perturbation 
of one-electron form 


P= Spl) ae) 2 ee ee 


the first term in eqn. (4.07), after modification as in (2.09) for lack of 
normalization, can be similarly simplified to become 


| PY*(H—W)¥ dr= ¥ | da*(1)... [p(¥*(6)]...(H-W)¥dr (5.06) 


and since each term of this has the form of eqn. (5.04) the sum vanishes. 
The second term in (4.07) is the Hermitean conjugate of this and so the 
whole equation is satisfied. Thus the self-consistent wavefunction is 
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stable under any one-electron perturbation provided that there is no 
degeneracy present. 

The degeneracy due to the spin restricts the states which can be described 
by a single determinant which is an eigenfunction of S2, to those in 
which all the orbitals are doubly occupied and those whose singly -occupied 
orbitals all have the same spin factor. In addition, any degeneracy due 
to the spatial symmetry has to be resolved in a similar way. 


§ 6. APPLICATION TO ENERGY CALCULATIONS 


The idea of using stable wavefunctions to discuss various perturbations 
is of special value in energy calculations. It enables some familiar 
results to be generalized and related to a few basic theorems. 

The simplest of these results is the quantum mechanical virial theorem 
(for full references see Léwdin (1959)). It states that the mean values 
of the total potential energy and the total kinetic energy of any atom 
or molecule are in the ratio 2: —1. An approximate wavefunction, even 
if its total energy is close to the true eigenvalue, does not necessarily 
give this same ratio. The theorem is proved, in the usual way, by 
introducing a variable scale-factor into the trial wavefunction. The 
equation for the optimal scale-factor is of the form (2.07) and consequently 
any function in the optimal scale is stable with respect to scale and so 
gives the correct ratio. Another method of finding a stable wavefunction 
is by the Ritz method using the functions (4.06) or (4.08). In the helium 
isoelectronic sequence, for example, with the wavefunction a simple 
product of orbitals of the form 


=exp (—1) Stet: oe Bes ae ee sni6, 03) 
the scale can be adjusted exactly by scaling or approximately by using 
Pca ( ES e7 N04 OM ne) en ote Sd, ERO 


which is a Ritz-type function with 6 as linear parameter and the second 


function is 
£ [exp (=A)] aan 5 kee) cept etn 6.05) 


in accordance with (4.08) when A is a scale factor. The equivalence 
between these can be demonstrated by attempting to apply both. The 
energy found using ?}, and the optimal scale for each value of Z, the nuclear 
charge, is, in reduced Hartree units, 


Was (Z—2 — 3.5? 258 4 2h t+ PUL + $644 ...). (6.06) 


This is insensitive to 5 and is least for 5=0. Thus, when the scale is 
optimized, the extra freedom introduced by the parameter 5 in fy does 
not improve the wavefunction since this type of improvement 
has already been used in the scaling. A similar situation is found in 
many more elaborate calculations in which scaling is used and its symptom 
is extreme insensitivity to a parameter or the presence of an over-complete 


set (Léwdin 1956). 
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Another well-known result is the theorem due to Hellmann and Feynman 
(1939), according to which the force on a nucleus can be expressed 
as the sum of forces of classical electrostatic type from the electronic 
cloud and the remaining nuclei. This theorem is derived by considering 
a displacement of the nucleus as a perturbation. Many approximate 
wavefunctions do not satisfy the theorem since they are not stable under 
this perturbation. Hurley (1954) has shown that the theorem is satisfied 
if ‘floating’ trial wavefunctions are used. This is not the only type of 
stable wavefunction. Since the perturbation is of the one-electron 
type self-consistent wavefunctions will also be stable in the circumstances 
discussed in §4. Approximate stability can also be attained by including 
in a Ritz-type trial wavefunction terms which correspond to an 
infinitesimal change of nuclear position. 

Any other parameter in the. Hamiltonian can be treated as variable 
and will yield corresponding results. Thus if one nuclear charge Z in 
an atom or molecule is allowed to vary then the theorem is deduced from 
(2.08) that the total potential energy (Z) due to the electrons and that 
nucleus is given by 

(Ly=ZeW eB 2 eee ee a, 
where W is the energy of the atom or the electronic energy of the molecule 
in the Born—Oppenheimer approximation. Equation (6.07) will also 
apply to approximate wavefunctions as long as they are stable under the 
perturbation of a change in Z. Since this is another one-electron per- 
turbation the self-consistent wavefunctions willbe stable. The major 
effect of a change in Z is to change the scale of the wavefunction and 
consequently for an atom the use of the optimal scale will make the 
wavefunction closely stable. The difficulties of convergence that arise 
when (ZL) is calculated using unscaled functions have been noted by 
Hylleraas and Skavlem (1950). 

For an atom in a 1S state the nuclear magnetic shielding constant o 
is very closely related to this potential energy. Lamb (1941) has shown 
that 


o=FX(>).- eee AIO 


when « is the fine structure constant. Thus by analogy with (6.06) 
the simpler relation 


a2 
==> 3 OW /0Z . . . . . . . (6.08) 


can be given. This equation needs some correction if W is the experi- 
mental energy and relativistic corrections are appreciable. Its importance 
is that o can be deduced from W without any quantum mechanical 
calculation. For the two electron series, Hylleraas and Midtdal (1956) 
have given the energy formula, in H R 
W = —Z?45Z/8— 0-157 655 + 0-008 535Z-1 — 0-000 34Z-2 
— 0-000 82Z-3—0-0024457-4 . . . . (6.09) 
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and this yields 
o = $a?(2Z — 5/8 + 0-008 535Z-2 — 0-000 68Z-3 — 0-002 46Z-4 — 0-009 787-5. 
(6.10) 
The probable error of c should be the same as that for W which is about 


0-00001 for Z>2. The values given by this formula are shown in table 1, 


along with those calculated by Ormand and Matsen (1959) and Hylleraas 
and Skavlem (1950). 


Table 1. Calculated nuclear magnetic shielding constants (o x 10°) 


HH He Lit Bet+ 
Hylleraas and Skavlem 5-9938 
Ormand and Matsen 2°3475 5:9856 9-5396 13-091 
This paper 2-432811 | 5:993347 | 9-541983 | 13-091305 


The kinetic energy associated with the different nuclei and with the 
electrons can also be deduced from the total energy by treating the masses 
as variable parameters. The electronic kinetic energy is given by 


CE=mdWlo(m)=—moWlom .-. . -. (6.11) 
and the kinetic energy of the nucleus of mass M by 
= LOW (Oi = ee eee (6.12) 


The hydrogen-like atoms, for example, have a total energy proportional 
to the reduced mass p and these equations lead immediately to the 
result 

GI LS en, eat seer (OL LS) 


§7. APPLICATION TO DrpoLE MoMENT CALCULATIONS 


When a molecule is placed in a uniform electric field of strength H in 
the z direction its Hamiltonian will contain the perturbing operator 


Pate kya oe mee Ee! eA (DOL) 


Since this operator is unbounded for large values of the variables the 
system is not strictly quantized but nevertheless, as Kato (1951) has 
shown, psuedo-eigenvalues may be defined and calculated by perturbation 
methods. The dipole moment of the molecule is defined, as in (3.04), by 
differentiating the energy with respect to # since this corresponds to the 
process needed to measure the dipole moment. ‘his derivative reduces 
to (z) only for a stable wavefunction. For a trial wavefunction to yield 
a good estimate of the dipole moment from <z) it is necessary, as in §4, 
that it should contain parameters which would enable it to adjust the 
molecule continuously to the presence of the field. When the wave- 
function is built from atomic orbitals the approximate effect of the field 


R 
P.M. 
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is to mix each orbital with one which is z times itself (cf. 4.06). Thus 
the adjustment of the 2s can be represented by including 2p, among the 
orbitals but the 1s function requires what has been called a Kp, function, 
ie. ap, function whose radial dependence makes it largest in the K shell. 
Similarly the 2p functions will require Ld functions. 

The importance of using stable wavefunctions can be illustrated from 
some recent calculations on LiH. These are summarized in table 2. 
It can be seen that the changes in the dipole moment as the wavefunction 
becomes more stable are considerable and out of proportion to the 
corresponding improvements in the energy. For more complete stability 
it may be necessary to include Ld functions on the Li. 


Table 2. Dipole moment of lithium hydride 


Atomic functions used <Z> im energy in 


Reference 
Debyes ev 


—0-87 | (—206-98)+ | Hurst et al. 
—6-04 | (—208-01)+ (1957) 


— 6-05 — 217-51 Karo and Olsen (1959) 
— 5:86 — 217-07 Platas and Matsen 
(1958) 


+ Energy value calculated from data in reference. 
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ABSTRACT 


Creep of vacuum annealed polycrystalline 65/35 alpha-brass was studied 
at constant stress in the range 77—300°K, close to the yield point, at elongations 
not exceeding 3%. The tensile strain ¢ increased logarithmically with time, 
and the characteristic strain sey =d«/d log, t, given by the slope of the creep 
curves in semi-logarithmiec representation, decreased approximately linearly 
with increasing temperature, from 29x 10-4 at 77°K to 11x 10-4 at 300°x. 
The creep was correlated with logarithmic stress relaxation, and the relation 
Ser=—Srel/x Was established, where y=da/de is the coefficient of work- 
hardening, and syej=de/dlog,,¢ the slope of the corresponding curve 
obtained with stress relaxation at constant elongation. Both the creep and 
the relaxation are ascribed to a process of residual slip propagation in which 
the migration of edge dislocations impeded by the drag of conservatively 
moving jogs is rate determining. The driving force on the jog is the 
component of the line tension of the dislocation, which impels the jog in the 
direction of its Burgers vector; the effective stress on the dislocation is the 
difference between that necessary to tear it from its locked site and the stress 
due to the intragranular sub-structure, which opposes its movement 
through the lattice in the unlocked state. At temperatures exceeding about 
300°K, when dislocation locking is weak, both ser and sye] tend to very small 
values. 


§ 1. INTRODUCTION 


In a recent paper on the low-temperature stress relaxation of polycrystal- 
line copper and alpha-brasses containing 10-35% Zn, at plastic strains. 
of up to 2% (Feltham 1961), the tensile stress was shown to relax from 
its initial value o, according to the logarithmic law 
GeO = SrellO0 i, (Layton gs dee ee ely 
where v is independent of the time, ¢, and 8re1, the most characteristic 
time-independent parameter of the relaxation, given by 
Sql OG ONOL 15 ba sie Me potas = ae ete st anne) 
increased approximately linearly with decreasing temperature. At any 
given temperature in the range investigated (77-300°K) sre: was found 
to have a minimum at a composition close to Cu;Zn, which suggested 
that the rate determining process was influenced by short-range order 
of the AuCu, type, occurring in the immediate vicinity of the cores ot 
dislocations. 
The main object of the present work was to study the creep of alpha- 
brass, again at plastic strains close to the yield point, and in the same 
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temperature range as in the case of stress relaxation, for the functional 
form of the creep law and the temperature dependence of the charac- 
teristic parameters were expected to give a more specific indication of 
the nature of the rate process operative in the creep and in the associated 
relaxation. Also, existing theoretical treatments of low-temperature 
creep in metals (Nabarro 1952, Mott 1953, Taubert 1959), and experimental 
studies (Davis and Thompson 1950, Wyatt 1953, Conrad 1958, Thornton 
and Hirsch 1958, Michelitsch 1959) were anticipated to be of considerable 
assistance in the interpretation of the results. 


§ 2. MaTERIAL AND METHOD 


Several rods of 65/35 brass, each 13cm long and of 2mm diameter, 
were together wrapped in thin brass foil of the same composition and then 
annealed at 730°c for $ hour in an evacuated silica tube, in which they 
were also allowed to cool after withdrawal from the furnace. The resulting 
grain size was approximately 4:5 x 10-?cm, as in the rods of brass from 
the same stock previously used in the work on stress relaxation. Zine 
losses during annealing were insignificant. 

Fig. 1 
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Stress-strain curve of 63/35 alpha-brass at room temperature. Grain size 
approximately 45x 10% cm. Lower curve: annealing temperature 
730°C; upper curve: annealing temperature 800°c. (Extensometer 
drum chart.) ; 
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After hard soldering the specimens into self-aligning grips the gauge 
length was 95cm. An account of the precautions taken in making the 
test pieces, as well as the analysis of the brass have been given previously 
(Feltham 1961). 

The stress-strain curve of an annealed specimen extended at room 
temperature is shown by the lower curve in fig. 1, and was obtained by 
direct transcription from the drum chart of the extensometer. The 
tensile yield stresses, cy, were close to the corresponding ones determined 
by Feltham and Copley (1960) with the same type of brass annealed at 
800°c for 12min, given in table 1. However, the initial part of the 


Table 1. Creep and relaxation data for 65/35 alpha-brass 


°K 77 200 273 | 289 | 299 
oy X 10-3 kg/cm? 14:3 10-0 ano 0 a 
Gx 10-5 kg/cm? 4:50 4-35 | 4:21] 4-20] 4-18 
Sre) kg/cm? 26 19 12 1m 10 
Ser X 10+ 29 23 15 12 11 
(Ser/Sre1) x 10+ kg-! em? 1-1 1:2 1-2 1-1 1-1 
G/(Sre1/Ser) 50 53 52 46 46 


work-hardening curve was flatter than obtained by them, shown by open 
circles in fig. 1. The steep portion found on annealing at 800°c is absent ; 
it appears to occur oniy in conjunction with a polygonized substructure, 
which does not develop readily below about 750°c (Doo 1960). The 
flat type of curve was also observed by Feltham and Copley in specimens 
annealed to the same grain size as at 800°c, but at a lower temperature 
~ (600°C). 

With the exception of the first 0-03 % of plastic strain, when the harden- 
ing rate is anomalously low due to a slight yield-point elongation, a 
constant work-hardening coefficient, y, was found up to about 2% of 
strain at all temperatures, with 

40 < G/y < 50, Rg ack Tie ENN Brae cod we Ad 
where @ is the shear modulus given in table 1. All tabulated data 
were obtained from creep observed at total plastic strains between 0-03 
and 2%. 

Specimens, completely submerged in liquid nitrogen, solid carbon 
dioxide in aleohol, or melting ice-plus-water, maintained at a constant 
level in a 35cm long tubular cryostat, were mounted in a Hounsfield 
Tensometer. The ‘125lb’ beam was used. If, at stresses exceeding 
the yield point, a strain increment (generally between 0-03 and 0-6 %) 
was applied to the specimen by the manual operation of the drive of 
the machine, the stress tended to relax, and the level of the mercury in 
the capillary of the load indicator diminished with time unless the stress 
was maintained constant at the initial value by the continuation of the 
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straining of the specimen. This was done by means of the manual 
control of the ‘fine’ drive of the reduction gear of the machine. A 360° 
protractor was mounted concentrically with the coarse drive, and was used 
for the measurement of the strain as function of the time; a 1° rotation 
was equivalent to an extension of the specimen by about 9. 
Some difficulty was experienced in holding the stress constant during 
the first. few seconds after the sudden application of a new stress 
increment, particularly at 77°K, when the initial rates of relaxation 
were high. Again, after about 200 sec the creep rate required to maintain 
the stress at its initial value was so small that the manual control of the 
stress was no longer fine enough to give a smooth compensation for the 
residual relaxation. Apart from these limitations the method was 
found suitable for determining the time-extension relations, and results 
were reproducible. 


Fig. 2 
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Creep curves of 65/35 alpha-brass annealed at 730°C, at plastic strains within 
about 3% of the yield point. Semi-circular symbols denote anomalous 
curves obtained in the yield-point elongation range up to about 0-03% 
total plastic strain. Curves denoted by other empty symbols refer to 
creep at plastic strains of 0-03-2%, and the full symbols refer to curves 
obtained above 2%, of total plastic strain. 


Creep curves are shown in fig. 2. All, at any given temperature, were 
obtained with the same specimen. After the determination of a curve 
due to one stress increment a new small load was added, the strain-time 
curve was again recorded, and so on. The extension was in each case 
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measured with the zero at the point on the stress-strain curve from which 
the new increment was applied. 


§ 3. EXPERIMENTAL RESULTS 
The principal features of the creep, apparent from fig. 2, were its 
logarithmic form, and the effect of temperature on the time-independent 
parameter Scr in the equation 
e=Scrlogij(1+vt), . . . . . . . =. (4) 
which was found suitable for the representation of the tensile creep 
strain. In fact, within the ranges of time amenable to observation by 
means of the equipment used, the simpler form 


CN SelOligviee ss os we () 
was adequate. Values of the strain soy obtained from the ae 
Ser =Ge/dlogigt,. 9. we ee wt eee (8) 


representing the slopes of the parallel lines drawn through some of the 
representative curves in fig. 2, are shown in table 1. Large values of 
Ser occur at low temperatures and vice versa. Reproducibility was best 
at total plastic strains of 0-03-2%. Anomalous curves were obtained at 
strains of less than about 0-03%, presumably because yielding was not 
taking place homogeneously throughout the entire length of the specimen 
and, as will be apparent later, because of the low coefficient of work- 
hardening near the yield point. Such curves are indicated by semi- 
circular points in fig. 2. Somewhat high values of scr were also found at 
total plastic strains exceeding 2%, again probably because of a decrease 
in the work-hardening rate below the level in the range 0-03-2% of 
strain. These anomalies were not pronounced; the strains at which 
they occurred exceeded those of the ‘normal’ range by at most 1$%, 
i.e. giving a total of 34%. 

The occurrence of esas times, if the strain increments were less than 
about 0:03°%%, apparent for example in the lowest curves at all tempera- 
tures, appears to be a general feature associated with “source-looking’ 
in stochastic processes of deformation. An interpretation of the 
phenomenon has recently been given by Yokobori (1959). 

The time-independent frequency v in eqn. (5) (and also in eqn. (15)) 
can be determined from the intercepts of the creep curves with the strain 
axis in fig. 2. At any given temperature its magnitude increases with the 
stress increment; a detailed study of the effect of changes of stress and 
temperature on v was not. however made. 


§ 4. THEORY 
4.1. The Relation between Creep and Stress Relaxation 


In order to be able to discuss both the creep and the corresponding 
low-temperature stress relaxation of alpha-brasses jointly, we shall 
first establish a general relation between them. Now, the rate of change 
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of the applied stress with time may be expressed in terms of the relaxation 
at constant strain and the variation due to work hardening: 


do 0c  Ocde 
— mk 
dt oF dedi ) 
which, for stress relaxation at constant strain, reduces to 
die ae 
dt oat dt 
while for creep under constant stress 
OS ea! 95 opener = SS eA eet 


Ga. digas 
x is the coefficient of work-hardening, which we shall assume to be inde- 
pendent of time and strain in the part of the stress strain curve used, so 
that 


=— =— =y= const. le 


On eliminating dc/0t between eqns. (8) and (9), and multiplying both sides 
by t, one obtains 
do de 
dlogyt ~*dlogit’ (11) 
so that, with a constant value of y, logarithmic creep and logarithmic 
stress relaxation are conjugates, in accord with our observations. 
Further, on substituting for the differential coefficients from eqns. (2) 
and (6), eqn. (11) yields 
Sre1/Ser = x: Shh eg ee Sc ee (12) 
The ratios G/y (table 1), where G is the shear modulus, were determined 
by the use of eqn. (12) with sp values measured in the present work, | 
and sre, as given for the same type of brass by Feltham (1961). They 
agree well with the corresponding ones obtained with x's derived directly 
from the stress-strain curves (eqn. (3)). 


4.2. The Equation of Logarithmic Creep 
Application of the theory of rate processes to creep gives for the shear 
rate 
dy 
= NAby, exp (= Gk) aa eee ee (13) 
where WN is the number of dislocation elements per unit volume contribut- 
ing to the creep, A is the area Swept out by one such element after 
activation, 6 the Burgers vector, and ¥ an atomic frequency. The 
tensile creep rate would be given by 


(14) 


where « is a constant fraction of the order 4-4, 
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In logarithmic creep the apparent activation energy may be expressed 
with sufficient accuracy, within the intervals of time readily amenable 
to experimental study of the flow, by 


U(t)=kTIn(vt), with v=NAbv/2-3sey, . .. (15) 
and if the assumption is made that it depends linearly upon the effective 


shear stress acting on the dislocation (e.g. Thornton and Hirsch 1958) 
one obtains 


> 


kT In ey rs) ee ete etc (16) 
where U, is independent of stress and time. The effective shear stress 
is equal to the difference between the applied shear stress 7 and the 
average opposing internal shear stress 7,, both resolved in the slip 
direction. The product v(+—7,) is a function of the time; v is generally 
termed the activation volume. If the shear stresses are replaced by 
corresponding tensile stresses, then in eqn. (16) 


OT 7.) —9(C-—Ce) eee os oe Rees) 
with 
2>q/U>3. 
From eqns. (16) and (17): 
deg nG kT (18) 
rinse Fry ee ere 
which, with eqns. (2) and (12), yields 
23hL 
Srel = (19) 
q 
and 
23hE 
Ser = . . . . . . . : : (20) 
qx 


As may be seen from eqns. (19) and (20), the assumption that the activation 
energy is linear in the effective stress (which we shall show to be consistent 
with the results obtained), and the relation between scr and sre) facilitate 
the evaluation of g, and hence of the activation volume v, as well as 
of the stress-dependent part of U(t), i.e. 

O= UG.) = days 4); ee et ee (et) 
where we have replaced o by oy, as the former was always within a few 
per cent of the yield stress in the experiments. Values of @ calculated 
from the relationt 

= (2:3kT'|Sre1)(ay — Fg); . . . . . . (22) 
obtained by substituting for g from eqn. (19) into eqn. (21), are given 
in table 2, together with sre; and oy taken respectively from the data of 


eee 

i ite ti lso Q> U 

As 7’ 0 creep will be observable at a finite time only if also @> Uy 

an (16)). ‘bie e therefore the energy supplied by the line tension in 

moving a jog without thermal activation (see §4.4). If it is assumed that 

o (1) =oy(400°K)G(T)/G(400), the values of Q in table 2 are decreased by 

g 

Oo 


= te) 
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Table 2. Activation energies derived from eqns. (22) and (23); given in 


keal/g atom 
Brass Cu } 90/10 80/20 70/30 65/35 
ly 5-3 8-0 5-6 5-0 a9 
Q; | 5-0 7-6 53 4-7 8-5 
200°K | Q ts 9-9 5-9 6-3 10-6 
Q, | 6-2 8-4 5-0 53 8-8 
299° K ir@ 7-0 oo fel) 4-6 11-0 
Q, | 6:5 ey 5-9 3°5 8-3 


Feltham (1961) and Feltham and Copley (1960). The stresses o, are 
shown in fig. 3, and were found by extrapolating the oy/T' curves of Feltham 
and Copley from 300°K to the virtually temperature-independent part at 
400°K; hence og=(oy)499.,-. A similar method of deriving o, in low- 
temperature creep has also been employed by Thornton and Hirsch 
(1958). We shall return to the discussion of the probable significance 
of a, below. 


Fig. 3 
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The intragranular elastic stress ¢,. Values indicated by points are equal to the 
tensile yield stresses of the brasses at 400°x, and were obtained by extra- 
polation from the data of Feltham and Copley (1960). 
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a“ 


4.3. The Rate Process 


The values of Q in table 2 may be in error by as much as + RUAN 
particularly those referring to room temperature, where o, approaches 
oy closely (eqn. (22)). It is nevertheless apparent that Q is not appreciably 
temperature dependent, which implies (eqns. (16) and (21)) that with ¢ 
equal to a few seconds In (vt) in eqn. (16) must be close to zero. Values 
of v found by Feltham (1961) in the work on stress relaxation referred to 
above were of the order of 10/min, which, on substituting into eqn. (16) 
shows that for small t, i.e. with o~oy We can write 

OD as = tas Peat ae eit (23) 
irrespective of temperature within the range studied. The small values 
of v are probably due to the small number of dislocation elements per 
unit volume contributing to the creep, and to the small area swept out 
by each such element after activation (eqn. (15)). Similar low values of 
the corresponding parameter were found by Davis and Thompson (1950) 
in precipitation hardened alloys. 

The actual magnitude of @ in copper, for example, is close to the activa 
tion energy determining the temperature dependence of the activation 
volume gq in high temperature creep, i.e. 6-6 keal/g atom, which Feltham 
(1957) associates with the conservative migration of jogs in dislocations. 
The lattice defect in the re-entrant corner of a jog in an undissociated 
edge dislocation resembles a small vacancy, which the dislocation must 
drag along on moving (Maddin and Cottrell 1955). A somewhat more 
complex defect would exist in an extended jog. Feltham (1957) suggested 
that the activation energy for jog migration should take the form 

Q=0,)(1—W Tah. 2) es oo (24) 
where the constant @; is equal to Q at 0°K, and 7p is the melting point 
of the metal. The values of Q,; obtained by means of eqn. (24) are also 
shown in table 2; 7m was in each case taken to be solidus temperature 
of the particular brass. For any given brass @, can be seen to be constant 
within the limits of experimental error. ; 

The mean values of Q;, shown in fig. 4, show a definite minimum close 
to the composition Cu,Zn. If jog migration is rate determining this 
could be interpreted as signifying a reduction of the misfit energy 
associated with the jog defect due to the interaction of the defect with 
solute atoms, enhanced by short range order localized in its immediate 
vicinity. It should be noted that @; does not appear to depend in an 
obvious manner on the stacking-fault energy, and it cannot therefore be 
attributed readily to processes such as the mutual intersection of dis- 
locations as visualized by Thornton and Hirsch (1958), or to cross-slip, 
the break-down of Cottrell-Lomer barriers, the formation of constrictions 
at the points of intersection of extended dislocations, and other processes 
requiring the recombination of extended dislocations. 

In relation to the short-range order it is interesting to note that 
Masumoto et al. (1952) obtained evidence of short-range order in 
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alpha-brasses at low temperatures by specific-heat measurements, and 
found that ordering ceased to take place at temperatures above about 
360°K, i.e. at about the same temperature at which the s;e; versus 7’ curves 
of alpha-brasses appear to converge to a common, low, value (Feltham 
1961). 


Fig. 4 
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Mean values of the activation energies given in table 2. 


Feltham and Copley (1960) suggest that the strong temperature 
dependence of cy below about 300°K results from a special form of 
Cottrell locking of dislocations, in which the tendency to the formation 
of short-range order near dislocation cores becomes instrumental in 
confining partial dislocations to close-packed crystallographic directions 
and that ordering is most pronounced in brasses with compositions 
close to CusZn. They algo show that o q (eqn. (22)) does depend somewhat 
upon the intragranular sub-structure resulting from any given heat 
treatment. This, and the observation that co, increases uniformly with 
the zine content, and hence with decreasing stacking -fault energy (fig. 3) 
as well as the small temperature dependence of cy above about 300°K, 
suggest that a significant part of the internal stress o, at temperatures 
es 2, Panes fe 300-400°K may be due to the interaction between stacking 
aults and solute atoms (and vacancies , of the t i i 
(1957) and Flinn (1958). ieee 

We believe that the foregoing observations, particularly the derivation 
of constant Q; values from eqn. (24) are consistent with the assumption 
that the driving force of the activated process in the low-temperature 
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creep and relaxation considered is the effective shear stress acting on the 
dislocation after it has broken away from its locked position, i.e. the 
strongly temperature sensitive part of the yield stress, represented by 
the bracketed term in eqn. (22). The locking is virtually ineffective 
in brasses above about 360°K. In copper the creep and relaxation are 
very small at all temperatures in the range, and the effective stress 
probably derives from the pinning of dislocations by residual impurities 
and vacancies. 
4.4. The Activation Volume 


An estimate of v (eqn. (16)) can be obtained by considering the migra- 
tion of a jog in an edge dislocation. Neglecting the dissociation of the 
dislocation, the force F on the jog J (fig. 5) is 


BF =27T" sin «, 


where « is the angle between the tangent to the circular dislocation arc 


Fig. 5 
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Schematic representation of the resultant force F acting on the jog J. The 
dislocation is bowed out into arcs by the effective shear stress 7-7 ,. 


at J and the radius of curvature, p, passing through J, / is the mean 
spacing between adjacent jogs, and 7” the line tension of the dislocation : 


TY = p(t 7). 
As sina=1/2p, we have 

F =(r—7,)Obl, 
and hence 

O= (G7, )074, 
or with 

ci, = 4(a—@,) 

and g=oy: 


ly/b =3Q/(cy—0,)b°, 


which, on substituting measured values, gives for 65/35 alpha-brass 
-lyb/~ 200 at 77°K, and approximately 500 at room temperature. 
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ABSTRACT 


A detailed analysis is given of examples of the combination, in crystals 
of silver bromide, of two dislocations on intersecting glide planes with Burgers 
vectors at 120°. The observations provide direct evidence for the occurrence 
of pencil glide in which slip can occur on any plane containing a <110> lattice 
vector. 


§ 1. IyTRODUCTION 


WE have recently described a method for sensitizing crystals of silver 
bromide which allows surface terraces and dislocation lines introduced 
during the earliest stages of plastic deformation to be made visible simply by 
exposing the crystals to radiation from a medium pressure mercury vapour 
lamp (Bartlett and Mitchell 1960 a). This permits the glide planes of the 
dislocations to be determined with reasonable accuracy. The Burgers 
vectors of the dislocations can then be deduced if the direction of the 
effective shear stress can be inferred. This method has already been 
used to show that two dislocations with Burgers vectors of 4a¢110) 
can combine in a silver bromide crystal to form a single dislocation with 
a Burgers vector of a{100) (Bartlett and Mitchell 1960 b). 

In the present work we have made a detailed study of two cases of 
the interaction between dislocations on intersecting glide planes with 
Burgers vectors at 120°. Many cases of such interactions have been 
observed near the surfaces of silver bromide crystals. 

This interaction, in which the intersecting dislocations are unable to 
form Lomer-—Cottrell dislocations, has been discussed by Whelan (1959) 
in terms of the Thompson nomenclature (Thompson 1953) with particular 
reference to stainless steel which deforms by glide on {111} planes. 
He observed frequent examples of the interaction when dislocations 
piled up on one glide plane were intersected by a dislocation lying on 
another glide plane. The position with silver chloride and silver bromide 
which deform by pencil glide (Nye 1949) is more complicated. Here, 
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as we have established by direct observations both in our previous work 
and in the present work, any plane containing a (110) lattice vector, 
the pole of which does not lie in the immediate neighbourhood of the 
pole of a {100} plane can act as a glide plane. To emphasize this, we 
present two examples of the particular interaction to be discussed here 
in which different pairs of glide planes are involved. For the description 
of the interactions, we prefer to use indices to define the Burgers vectors 
and the glide planes of the dislocations as the Thompson notation is less 
convenient when the glide planes are not {111} planes. We use an arrow 
to show the positive sense of the dislocation line and write the Burgers 
vector beside the arrow. 


Stereographic projection of the crystal of fig. 1. 


§ 2. EXPERIMENTAL OBSERVATIONS 


Full details of the experimental methods have been given in two 
previous papers (Bartlett and Mitchell 1960 a,b) and nothing further 
need be added here. The systems of dislocations form networks in 
three dimensions and considerable difficulty has been experienced in 
producing a satisfactory record at the necessary high magnification. 
This problem has been partially resolved, where necessary, by taking 
a series of photomicrographs with the microscope focused successively 
on planes at an increasing depth below the surface. The sharply focused 
sections were then cut out from enlargements of the negatives and 
assembled to form a mosaic. 
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The first interaction which we will describe is shown in the photo- 
micrograph of fig. 1¢. The orientation of the crystal was determined 
from a back-reflection Laue x-ray diffraction photograph and the pole 
of the surface, P, is plotted in the stereographic projection of fig, 2. 

We have found that the glide planes which we determine experi- 
mentally are always near planes containing a (110) lattice vector. This 
can be seen on the stereogram where we have plotted the poles of the 
glide planes. The orientations of the glide planes were deduced from 
observations with the microscope of the relative positions of the surface 
terraces and of the corresponding dislocation lines with respect to the 
known orientation of the crystal. The lengths of the arcs represent 
the error in the observations. We therefore conclude that the two 
approximately parallel dislocations were introduced by glide on a pair 
of planes with poles very near the pole, E, of the (323) planes. As the 
specimen was deformed by the application of very small tensile stress, 
the dislocations must have a Burgers vector of }a[101]. They are crossed 
by a third dislocation on a (011) plane with a pole at F on the stereogram. 
This dislocation has a Burgers vector of 4a[0I1]. It combines with the 
other two dislocations lowering the elastic strain energy to form a 
dislocation with a Burgers vector of }a[110] along the line of intersection 
of the glide planes ; this is a [533] direction which is represented by 
the point G on the stereogram. This deduction is consistent with the 
experimental observations, the error in which is indicated by the length 
of the are at G. The glide plane of the resultant dislocation would be 
a (338) plane but the dislocations meeting at the nodes could not move 
simultaneously on their three intersecting glide planes so that the nodes 
are effectively immobile. The results of the analysis of this interaction 
are summarized in the diagram of fig. 3. 

A second interaction of the same type is shown in the photomicrograph 
and diagrams reproduced as figs. 4, 5 and 6. In this crystal, dislocations 
with Burgers vectors of }a[011] and 4$a[101] on (111) and (212) glide 
planes, respectively, are involved. The poles of these glide planes are 
indicated, as before, by E and F in the stereogram. The dislocations 
combine to form a dislocation with a Burgers vector of 4a[110] along 
the line of intersection of the two glide planes which is a [143] direction. 
This is shown at G on the stereogram. The glide plane of the combined 
dislocation would therefore be a (335) plane. The validity of this analysis 
of the interaction was confirmed by the analysis of further interactions 
between these dislocations and other intersecting dislocations in the 
surrounding area. At A in fig. 6, the dislocation on the (212) plane 
interacts and annihilates with a dislocation with the same Burgers 
vector on a (111) plane. The pole of this plane is plotted on _the 
stereogram at H. At B and C in fig. 6, the dislocations on the (111) 
and (212) planes are intersected by a third dislocation on a (111) plane 
ty pe ae 


+ Figures 1 and 4 are shown as plates. 
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Fig. 3 
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Diagram showing the Burgers vectors and glide planes of the dislocations in 
fig. 1. The dotted lines are the projections on the plane of the diagram 
of the traces of the intersections of the glide planes with the surface 
of the crystal. 


Fig. 5 


Stereographic projection of the crystal of fig. 4. 
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represented by J in the stereogram with a Burgers vector of Za{101]. 
A further example of an interaction between two dislocations with Burgers 
vectors at 120° therefore occurs at B, while the two dislocations mecting 
at C combine to form a segment with a Burgers vector of a[100] (see 
Bartlett and Mitchell 1960 b). 

The one possible alternative assignment of Burgers vectors was excluded 
as being inconsistent with deformation by the application of a tensile stress. 


Fig. 6 
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Diagram corresponding to fig. 4. 
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ABSTRACT 


The thermal expansion of aluminium has been determined by a sensitive 
optical method between 100°K and 15°x; that is, down to temperatures below 
©p/20. At the lowest temperatures the Griineisen parameter shows a 
significant rise. The results are discussed from the point of view of the 
theory of lattice vibrations. 


§ 1. INTRODUCTION 


In the past few years there has been renewed interest in the thermal 
expansion of crystals at low temperatures (Rubin ef al. 1954, Bijl and 
Pullan 1954, 1955, Barron 1955, Domb and Zucker 1956, Figgins ef al. 
1956, Simmons and Balluffi 1957, Blackman 1958, Dheer and Surange 
1958, Gibbons 1958). The importance of the coefficient of thermal 
expansion 


is that it is one of the three independent derivative thermodynamic 
parameters which are experimentally accessible. In the dynamical 
theory of crystal lattices its importance is that it leads to information 
about the variation with volume of the frequency distribution of lattice 
vibrations. 

The expansion coefficient « is best studied via the Griineisen parameter 
y= aV/c,«, which is of course constant for simple isotropic solids over 
a wide range of temperature. On the Debye model this implies the 
Griineisen assumption that for all frequencies (v;) the quantities 
y;(= —dInv,/dIn V) are equal and given by y=—dIn@p/dinV. There 
have been few attempts to predict theoretically the variation of « at low 
temperatures and only that by Barron makes use of a proper lattice 
theory. 

Experimentally, the situation seems to be that for simple isotropic 
solids y remains nearly constant at least down to @p/10. Below these 
temperatures the facts are still uncertain because the values of w to be 
measured (~10-® for aluminium, say) are so small. In a previous 
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paper (Figgins et al. 1956) measurements were reported on aluminium 
using x-ray diffraction for the measurements of lattice parameter. ‘This 
method gave good results down to ©)/10(#40°K), and it is, of course, 
an absolute method. In the present paper we describe measurements 
using a new optical method giving good results below ©)/20 and capable 
of use at still lower temperatures. 


’ 


§ 2. EXPERIMENTAL MreTHop 


Previous methods for the measurement of « at low temperatures have 
included the use of x-ray diffraction (Figgins et al. 1956, Simmons and 
Balluffi 1957), an optical lever (McLennan e¢ al. 1931), interferometric 
techniques (Rubin et al. 1954) and capacitative dilatometers (Bijl and 
Pullan 1955, Dheer and Surange 1958). 

The present method makes use of the optical lever principle, but 
employs a microptic autocollimator (Hilger and Watts), a kinematically 
designed twin-strip mechanism (Geary 1954), and an optical system 
designed to eliminate all errors due to incidental movement or to thermal 
expansion of parts of the apparatus other than the test specimen. The 
ultimate sensitivity in measuring the linear thermal expansion of a 
5em specimen corresponds for the optical lever used to a displacement 
of about 304, that is, defines an « of 1 x 10-® with an accuracy of about 
20% with a temperature interval of 1°K. 

The apparatus is shown schematically in fig. 1. Expansion of the 
test specimen relative to a control specimen causes the rotation of a mirror 
attached to a twin-strip mechanism held between them. The two 
specimens are mounted on a glass plate, and a second mirror, inclined 
at about 90° to the first, is carried by the control specimen. The auto- 
collimator, having an ultimate sensitivity of 0-3” of arc, indicates the 
angle between the mirrors through the separation between the images 
formed by the two beams of light passing between the mirrors in opposite 
directions. To the first order the separation between the images depends 
only on the angle between the mirrors and is not affected by incidental 
rotations or movements in the system. The field of view of the instru- 
ment is only 10’ of are so that to keep the images suitably placed it is 
necessary to make adjustments from outside, using control rods which 
can be disengaged. In separate experiments at room temperature the 
device was checked for linearity. 

The assembly is set up in an evacuated cryostat with anti-vibration 
mountings. By means of helium as exchange gas the comparator system 
can be brought into thermal equilibrium with the cryostat, whose tem- 
perature is controlled by balancing the Joule—Kelvin cooling of hydrogen 
against electrical heating. 

The comparator is absolute in the sense that it measures the relative 
expansion of two specimens of the same metal at different temperatures 
and does not depend upon the expansion coefficient of another material. 
There is however a temperature gradient near the specimen so that it is 
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Schematic diagram of apparatus. 


essential to use a dynamic method—somewhat similar to that used in 
the ‘ discrete heating ’’ method of determining specific heats—in which 
temperatures and deflections are continuously watched. 

Temperatures were determined using copper—constantan thermo- 
couples calibrated at the steam, ice, oxygen and hydrogen points and at 
the hydrogen triple point. At the lowest temperatures helium gas 
thermometers, arranged differentially, may also be used. 

The sensitivity of the optical lever depends critically upon the separation 
of the strips and so was determined in situ by making use of a value of 
a at room temperature of 7-05 x 10->deg-c taken from the tabulated 
data of Nix and McNair (1941). Cooling the mechanism from room 
temperature to 20°K cannot alter the sensitivity by more than 4%. 
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The aluminium was stated to be at least 99-99% pure (British 
Aluminium Company), and was poly-crystalline. 


§ 3. RESULTS 


The results are shown in fig. 2, where each circle represents a separate 
experimental determination. Because of the difficulty of measuring 


Fig. 2 
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Coefficient of thermal expansion of aluminium as a function ot temperature: 
results of the present investigation. 


small values of « with accuracy a greater number of experiments were 
performed at the lowest temperatures. We believe that the values 
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represented by the curve as drawn are accurate to within about + 15% 
at 100°K, and +30% at 15°x, a 


Fig. 3 
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Coefficient of thermal expansion of aluminium: comparison of present results 
with those of other workers. 


In fig. 3 the same curve is drawn, with values plotted from the results 
of Bijl and Pullan (1956, private communication), Figgins et al. (1956) 
and Gibbons (1958). The only other published results on aluminium 
near @p/20 are those of Gibbons and it will be seen that our results are 


‘in disagreement with his at low temperatures. 
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In order to obtain a value of the Griineisen parameter characteristic 
of the lattice only, we have used the formula 


aV 2 
n= (— = 30) [teen 
Kp 


(Bijl and Pullan 1955) where c, is the electronic specific heat. As can 
be seen in fig. 4, we find a marked rise in y, below @p/10 which we believe 
to be well outside experimental uncertainty. Rounded values of « and 
y, at intervals of 5° or 10° are shown in the table, together with the data 
used in the computations, and their sources. 


Fig. 4 
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Variation with temperature of the Griineisen parameter (y;) of aluminium. 


§ 4. Discussion 


The theory of Barron computes y as a weighted average of the y, taken 
over the different modes of vibration, and using a central force model 
predicts a drop in y near 0:20). A general treatment making use of 
the Debye—Brillouin model and leading to restricted conclusions is that 
of Bijl and Pullan (1954). They adopt the method of assuming the same 
cut-off wavelength for longitudinal and transverse waves and thus two 
different cut-off frequencies. This leads to the existence of a high- 
temperature and a low-temperature limit for y, and they conclude that 
the low-temperature limit will be the greater}. An attempt is made 


} The definitions of y employed by Bijl and Pull 
identical, but their conclusion must ee. both ae Sn Ly ails 
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to compute the effective contribution y, of the lattice only by assuming 
that the electronic contribution corresponds to the free-electron value 
of % for y,. We have used the same procedure, and it is worth noting 
that at the lowest temperatures the differences between y and y, are 
significant, as may be seen in the table. The main result of the present 
work is that y, rises at the lowest temperatures, as suggested by Bijl 
and Pullan. 

The deficiencies in the present theoretical situation are particularly 
obvious for the case of metals and it is worth considering whether any 
improvement can be made by empirical means. It does not seem possible 
at present to improve on the somewhat over-simplified procedure of 
extracting y, from the data by assuming a value for y,. It is however 
possible to insert certain data on the variation of elastic constants with 
pressure due to Lazarus (1949). On the assumption that the y,; take 
constant values y,,,, and y;,,.,, for the longitudinal and transverse 
modes of vibration respectively in the Debye—Brillouin model this leads 
to a variation of y, with temperature given by 


= 2Y trans. Cl Oren) + Yiong. od (Set £8) 
26,(@mantL) + 6(Ormeae) , 


where 
) a Cr 3Caa UCuy 
trans. 2 Cay d p ? 
1) Cua $Caa Fey | 
Yiong. 2 { Cr d ip > 
and 
Op?* rang Oe : hong. C11 


Using Lazarus’s figures, this leads to an increase of y, with falling 
temperature, the low-temperature limit being about 8% higher than 
the high-temperature limit. The main rise in y, takes place below 
0-30p but above 0:10, where the rise in our experimental values begins to 
appear. There is thus no more than a rough qualitative agreement 
between this semi-empirical conclusion based on the Debye-Brillouin 
model and our own results. 
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Note added in proof.—Our results are in good agreement with those of 
D. B. Fraser as reported to the VIIth International Conference on Low 
Temperature Physics (September 1960). The question of the effective 
value of y, has recently been discussed by Klemens (1960) and Gross and 
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Olsen (1960). The latter have determined the pressure shift in the 
superconducting transition temperature of aluminium and estimate 
y-+10. Our results at the lowest temperatures can be interpreted in 
terms of a constant y, +6 rather than 3. 
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ABSTRACT 


An account is given of the theory that a crystalline or non-crystalline 
array of atoms will make a transition from the metallic to the non-metallic 
state as the interatomic distance is varied. Experimental evidence relevant 
to the theory is summarized. Wigner’s (1938) suggestion that an electron 
gas might crystallize at sufficiently low densities is examined, and an experi- 
ment suggested by which it might be tested. 


§ 1. INTRODUCTION 


Ty anumber of papers (Mott 1949, 1952, 1956, 1958, Mott and Stevens 1957), 
the present author has put forward the hypothesis that a crystalline array 
of hydrogen-like atoms, or more generally atoms in which there is an 
incomplete shell, may not necessarily show metallic conduction, by which is 
meant a conductivity which does not vanish at the absolute zero of tempera- 
ture. The hypothesis states, moreover, that as the lattice parameter of the 
crystalline array is increased, there will at the absolute zero of temperature 
be a sharp transition from a metallic state, with a finite or infinite 
conductivity at the absolute zero, to a non-metallic state in which it 
will not conduct at the absolute zero and in which there is a finite activation 
energy for conduction. Much earlier Wigner (1938) suggested that an 
electron gas of low concentration might ‘crystallize’, the electrons 
becoming localized. These hypotheses have proved controversial and are 
not, I believe, universally accepted. There is however a good deal of 
new experimental evidence in the favour at any rate of the first, and the 
purpose of this paper is to restate this hypothesis, to summarize the 
evidence and to suggest an experiment by which Wigner’s proposal 
could be tested. 


§ 2. THEORY OF THE TRANSITION TO THE METALLIC STATE 
(DivALENT METALS) 


Since A. H. Wilson’s papers of 1931, an insulator has been described by a 
Slater determinant formed from Bloch wave-functions y;,(r) such that the 
wave vectors k are all within totally occupied Brillouin zones. Thus in a 
simple case, for instance a divalent cubic metal, one might envisage a full 


zone (f) and an empty zone («) as shown in fig. 1, with an activation energy 
PAs Staak ie 
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E needed to excite an electron from the full to the empty zone. As one 
compresses the crystal, decreasing the lattice parameter from its equili- 
brium value b,, one would expect both bands to widen so that # would tend 
to zero at a value of 6 equal to by. For values of b smaller still the bands 
overlap and the model predicts metallic conduction ; thereafter the number 
of free electrons increases uniformly. 

The first point that we must make about this model is that the predicted 
continuous increase in the number of free electrons and holes from the 
value zero is not possible. An electron and a positive hole will attract 
each other with a force behaving at large distances r like e2/«r?, where x is a 
dielectric constant. With a Coulomb force of this kind, the electron and 
hole will always, when in the state of lowest energy, form pairs (excitons, 
like positronium), with binding energy of order 
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Showing a full zone («) and an empty zone (8) varying with the interatomic 
distance b. 


where m* is their reduced mass. They will not therefore be able to carry a 
current at the absolute zero of temperature. It is only possible to escape 
from this conclusion if there is a considerable concentration, n per unit 
volume, of free carriers; the Coulomb field will then be replaced by one 
derived from the potential energy 


“=(¢*/er) oxp (Gr) eee eT) 


where q is a screening constant given in the Thomas—Fermi approximation 
by (Mott and Jones 1936, Chap. IT): 


Q? = 4men' 3 /ch?, Me es re ee) 


A screened potential of this type does not necessarily lead to any bound 
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states. We therefore predict that the number of free carriers (electron and 
holes) varies discontinuously} with 6, as shown in fig. 2. 


§ 3. THEORY OF THE TRANSITION TO THE METALLIC STATE 
(MonovaLent METALs) 


We turn to the consideration of materials which, according to the Wilson 
model, ought to be conductors for all values of b; the simplest is an arrange- 
ment of one-electron centres (hydrogen atoms) arranged on a cubic close- 
packed cubic or similar lattice. According to the Wilson model, the zone 
will be half full and the material a metallic conductor for all values of b, as 
in fig. 3, though the effective mass will increase with b. 


Fig. 2 


\/b 


Plot of the predicted number of free electrons of (say) a divalent metal against 
the reciprocal of the distance b between the atoms. 


Now for large values of b this is against common experience and, one 
might say, common sense. Suppose that the overlap between the centres 
is small. Then clearly one would expect a current at low temperatures 
only if one or more of the electrons were missing. Into the ‘hole’ so formed 
an electron from the neighbouring site could move and the process could 
continue ; the hole would behave like a carrier of high effective mass. But 
if all centres are occupied by an electron, then for a current to flow it is 
necessary to ionize at least one atom and to put the electron back on a 
distant atom; two free carriers are thereby obtained. The work necessary 
for this is to a first approximation I — #, where J is the ionization energy and 
FE the electron affinity. Thus an activation energy of this order is necessary 
for conduction. We want to emphasize that the presence of polar states in 
the wave function, such as one finds to some extent in the hydrogen molecule 
however small the overlap between atomic orbitals, does not lead to con- 
ductivity ; this point will be considered further in § 4. 


+ The assumption that 7 changes discontinuously may break down if « tends 
to infinity at the transition value of b. This can only occur if the energy of the 
optically allowed transition from the full to the empty zone tends to zero. 
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These findings accord with what we know about the conductivity of 
some of the transition metal oxides (e.g. NiO, Morin 1954). This particular 
oxide has a cubic structure; the nickel ions (Ni?+) should have the con- 
figuration 3d8; according to the Wilson model, the 3d band should be 
partly filled. Nevertheless this oxide when pure is an insulator and 
transparent; it conducts only when activated, say by Li* ions replacing 
some Ni2+ ions, with consequent formation of Ni** ions to preserve neutrality 
ofcharge. This shows that the nickel ions overlap sufficiently to allow some 
current to pass, but that the oxide only shows conductivity when Ni** ions 
are present. 


Fig. 3 
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Showing a half full zone for a monovalent metal varying with the interatomic 
distance 0. 


We consider therefore that in the array of one-electron centres, as in the 
array considered in § 2, free carriers are only possible if they are present in 
such numbers as to give a screening constant g big enough to prevent the 
formation of pairs. We come therefore again to the conclusion that, as the 
interatomic distance b is varied, there must be a sharp transition between a 
non-metallic state which can carry no current at the absolute zero and one in 
which a large number of electrons are free. Since the latter is the well- 
known metallic state, it is reasonable to expect all the electrons to be free. 
A very approximate estimate can then be made of the value of 6 at which 
the transition will occur. 

If the potential energy of a positive charge and an electron is 

— (e*/«r) exp (—qr), 
the condition that pairs will not form is (Friedel 1952, Blatt and Weisskopf 
1952) | 
q > 1-0/(h7«/me?). ree eee ee 
We do not know the effective masses of the ‘hole’ or the electron under 


discussion ; m should be the reduced mass of these which we shall take equal 
to the electronic mass. Using formula (3) for q, (4) then gives 


Vedas Oa. snes 
Shea H . . e . . (5) 


yy = h?«/me?. 
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If we are dealing with an array of hydrogen atoms, the ‘background’ 
dielectric constant is of course unity and we see that free electrons do not 
form until the interatomic distance is two or three multiples of the hydro- 
genic radius. If the centres are impurity centres in semiconductors, « 
should be the dielectric constant of the semiconductor, and m should be an 
effective mass of an electron (or hole) in its conduction band. 

Formula (5) is doubtless dimensionally correct but the derivation of the 
numerical constant contains too many crude assumptions to be taken very 
seriously. It is, however, interesting to apply it to bismuth, where the 
number of free electrons and holes is certainly small, about 10-5 per atom 
(Smith 1959). On the other hand, the effective masses are small, and the 
background dielectric constant x, due to the very nearly full bands, is 
likely to be large. There can be little doubt that (5) is satisfied. 


§ 4. A More Detattep Discussion or THE Non-conpuctinG WAVE 
FUNCTION 


A detailed discussion of the non-conducting wave function and a proof 
that a system described by it does not conduct may be given as follows 
(Mott 1956, 1958). We consider a crystalline array of N one-electron 


Fig. 4 
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(hydrogen) atoms. For an electron moving in a periodic field due to N 
atoms, there are N Bloch wave functions y¢,(r) possible in the first band for 
each spin direction. The usual description of the wave function of a 
monovalent metal is a determinant or sum of determinants using half 
these wave functions, the half with the lowest energy. These wave func- 
tions describe a state which can carry a current because the excited states 
of the system lie very close to the ground state, being separated from it by an 
energy which tends to zero as N tends to infinity. It is thus possible, for 
energies of the whole system very near to that of the ground state, to set up 
complex wave functions, which can describe a state with a current. 

It is, however, possible from the N wave functions in the first band to set 
up a set of orthogonal wave functions, the Wannier (1937) or Lowdin (1950) 
functions w,(r) which are localized, one on each atomic site ; here the 
suffix » denotes the site. Such Wannier functions are shown in fig. 4. 


T2 
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They decay exponentially as one moves away from the site in question, 
but there can be considerable overlap onto neighbouring sites. 

It must be emphasized that to set up Wannier functions for each lattice 
point, one must make use of all the Bloch wave functions in the band, 
We may use these Wannier functions to describe a state of the system which 
is antiferromagnetic or ferromagnetic and is non-conducting. If the 
material is antiferromagnetic we may set up for each site a wave function 

W,(t) X,(0) OF Wy(r) Xe(), 
the functions y,, y, being spin wave functions for the two spin directions, and 
from these we may form a single determinant to describe the state of the 
whole system. The properties of this wave function are: 


(a) It is real, and therefore does not describe a state carrying a current. 


(>) The only state with the same or nearly the same energy is that in 
which all the spin directions are interchanged, except states formed by the 
excitation of spin waves, which are believed not to carry a current. 


(c) Ifit is treated as an approximation to an exact wave function, we see 
that the first excited state with zero spin will have an energy of order J— EH 
(several electron volts) above the ground state. 


(d) Since the wave functions w,(r) overlap, polar states are present, in 
the sense that the wave function allows the presence of two electrons on one 
atom; but since the wave function is real this cannot lead to a current. 
The polar states therefore do not lead to mobile charge carriers so long as the 
state of the crystal is correctly described by this kind of wave function. 
This is because the holes and electrons in the lowest state wave function are 
bound to each other in space. 

For crystals containing two electrons per atom, it is well known that the 
determinant formed from Bloch wave functions and that formed from local- 
ized (e.g. Wannier) functions are identical, both describing a non-conduct- 
ing state ; but for crystals with one electron per atom they are not the same. 
The above considerations show that the two approximate wave functions 
then correspond to two different states of the crystal with different exact 
wave functions, one conducting (and with an electronic specific heat linear 
in 7 and with other metallic properties), and the other non-conducting. 
They are thus not alternative zero-order descriptions of the exact function ; 
they correspond to different states of the system, the true ground state (at 
the absolute zero) in each case being that with the lower energy. The wave 
function formed from Wannier functions will clearly give the lowest energy 
for large interatomic distances; the existence of metals leads us to believe 
that the other gives the lower energy for small interatomic distances. 

We have assumed up to this point that the non-conducting wave function 
describes an antiferromagnetic arrangement of spins. There are other 
possibilities—that it may be ferromagnetic or paramagnetic without any 
long-range order in the directions of the spins. We do not know of any 
non-conducting material corresponding to the second case, nor whether 
such a case can exist. For ferromagnetic non-conducting materials the 
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argument would be the same as that given above; one can set up a single 
determinant as a first approximation to the wave function, so that the 
wave function is real and cannot describe a state in which a current flows. 

In either the ferromagnetic or antiferromagnetic case, by reversing the 
spin of one or more electrons and setting up the wave functions correspon- 
ding to a spin wave one can obtain excited states with excitation energies 
tending to zero as N tends to infinity. One can therefore set up complex 
wave functions. But itis assumed} that aspin wave cannot carry a current, 
the current vector corresponding to a wave function describing a spin wave 
being zero. This theoretical prediction corresponds to the observed fact 
that typical antiferromagnetic materials show an exponential increase in 
conductivity (as exp(—W/kT)) as the temperature is raised, while the 
energy of spin waves excited is believed to be proportional to 74 (Anderson 
1950). 

Finally, we want to emphasizé that the sharp transition described here is 
only expected in an infinite lattice. It goes without saying that for a 
finite number of atoms there will be a gradual decrease in the weight of the 
ionized states in the wave function as the interatomic distance is increased, 
or in other words a gradual transition from the LCAO model to the London- 
Heitler model. 


§ 5. SLATER’S MODEL FOR SUBSTANCES OF THE NICKEL—OXIDE 
TYPE 


Slater (1951) has suggested that antiferromagnetic substances of the 
nickel—oxide type should be described by Bloch wave functions and their 
non-conducting properties described as follows. Any electron described 
by a Bloch wave function will, in terms of the Hartree-Fock equation, be 
repelled considerably less by electrons with spins in the parallel direction 
than by electrons with spins in the antiparallel direction. Therefore in a 
crystal in the antiferromagnetic state each electron will move in a crystalline 
field on which is superimposed that of a superlattice. This field of a super- 
lattice will split the first zone of a cubic lattice into two, each containing 3N 
states for each spin direction. With one electron per atom, the first zone 
could be completely full and the second empty. The material would then 
be an insulator. 

The suggestion that the band model can lead to a non-conducting state 
in this way is certainly correct. It can also lead to localized magnetic 
moments, both antiferromagnetic or ferromagnetic, and with integral or 
non-integral spin. This seems first to have been emphasized by Lidiard 
(1954). A recent paper by Friedel et al. (1961) discusses the transition 
metals in terms of these concepts, and shows that, always retaining the 
band model and if necessary non-integral moments, localized moments 
may survive above the transition point. This suggests that as one 
increases the interatomic distance of an array of one-electron atoms, the 


+ If this assumption is false, our argument falls to the ground. 
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first thing that happens may be the appearance of localized moments 
without the disappearance of conductivity. 

What is misleading in Slater’s model is the following. It suggests 
incorrectly that the crystal will show metallic conductivity above the Néel 
temperature or in other words that the activation energy required to 
destroy the long-range antiferromagnetic order and that to produce carriers 
of current are of the same order of magnitude. In fact the former will 
depend on some overlap integral giving the strength of the exchange or 
superexchange interaction, which will tend to zero as 6 tends to infinity. 
The activation energy for conduction, on the other hand, tends to J—#. 


Fig. 5 


A Brillouin zone split in two by a superlattice, giving sub-zone boundaries 
shown. by the dotted lines. 


‘ 


Nonetheless, Slater’s model of the non-conducting state must be the same 
as that set out in §4 at the absolute zero of temperature. It is interesting 
to note that it, too, predicts a discontinuous change in n, the number of 
free electrons. This may be seen as follows. 

Figure 5 shows schematically a Brillouin zone split into two by the super- 
lattice lines due to an antiferromagnetic arrangement of spins. Then 
suppose that the material is still a conductor and there is a small overlap rp, 
at B into the second zone and an equal number of holes in the first zone at A. 
What we have to show is that a small number is impossible. This is because 
the states within the dotted zone near its boundary are those in which the 
wave function of an electron has its maxima in those atoms with the same 
spin and nodes in the planes containing atoms with opposite spin. The np 


electrons at B are electrons which spend most of their time in atoms with 
opposite spin. 
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Thus if the ng electrons are transferred to A, the moment on each atom 
will be increased by about npu, where pz is the Bohr magneton, and the 
width of the energy gap along the dotted lines will be increased y a term 
proportional to nz. Thus the energy of the whole system will be lowered 
by a term proportional to ny. On the other hand, the work done to transfer 
the electrons is proportional to n,°!*, since their Fermi energy is propor- 
tional to n;?". Thus, if , is small enough, energy must be gained in 
transferring electrons from B to A. This model, therefore, shows that a 
discontinuous change in the number of free electrons is to be expected at 
the value of b for which the non-conducting state occurs. 

Other discussions in the literature on the transition to the metallic state 
are given by Des Cloiseaux (1959), Frood (1960) and Coldswell-Horsfall 
and Maradudin (1960). These will not be discussed here. 


§ 6. INTERATOMIC DISTANCE AT WHICH THE TRANSITION OCCURS 


A rough estimate of this has been given in § 4, based on the concentration 
of electrons at which screening is strong enough to prevent pair formation. 
Although this gives satisfactory results, it should not perhaps be taken too 
seriously. Neither should a treatment based on the Slater model given by 
the author (Mott 1956). 

Ideally, it should be possible both to calculate the energy of a degenerate 
electron gas in the field of an array of ions, and the energy of an array of 
overlapping hydrogen-like atoms with an antiferromagnetic spin arrange- 
ment. With the positive charge smeared out this has been done by Carr 
et al. (1961); the accuracy is however doubtful, since they use the method 
of Gell-Mann and Brueckner (1957) which converges badly for a low value 
of the electron density. One would expect curves of the form shown in 
fig. 6. ABC is the normal curve for the energy of metallic electrons in an 
array of hydrogen-like atoms, calculated for instance by the Wigner—Seitz 
method with full allowance for correlation ; ED is for the antiferromagnetic 
array of hydrogen-like atoms. Since the curve is meant to represent the 
lowest state of the whole system, the dotted parts have no physical meeting, 
but are just extrapolations of hypothetical calculations. The transition 
should occur when the two curves cross. It is an assumption that dH/db 
is discontinuous there, but if there is a discontinuity in the number of free 
electrons it would be surprising if it were not. Moreover, we give in §7.1 
some evidence that it is discontinuous. 

If this is so a sharp transition should also occur at finite temperatures at 
the value of 6 for which free energy curves similar to those in fig. 6 cross. 
Since the entropy of the electrons in a metal increases as 7’ and that due to 
spin waves among localized electrons as 7°, a rise in temperature will favour 
the metallic state and increase the value of b for which the transition occurs. 
A suitable analysis would be the following: let 


Ey (V— Vo) Vos Ew (E= VilVo 
be the energies per electron at 7’=0 in the metallic and non-metallic states 
at a volume V per atom differing slightly from the transition volume Vy. 
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Fig. 6 


Energy of the metallic state (ABC) and the non-metallic (ED) as a function of 
interatomic distance b (schematic). 


metallic non-metallic 


A b 


The full line shows the boundary in a 6, 7 diagram between the non-metallic 


and metallic states. .A crystal with a value of b denoted by A in the 
figure will make the transition at the temperature 7',. 
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Let H, be the Fermi energy in the metallic form, so that the electronic 
specific heat is 

tn ke TLE. 
Then, neglecting the spin wave term in 74, the difference in free energy 
between the two forms is 


(Ey — Eymu)(V — Vol Vo— a7 kT? |B. te eo) 
We therefore expect a boundary between the two forms as shown in fig. 7. 
These ideas will be applied to certain oxides in the next section. 


§ 7. APPLICATIONS 
7.1. Oxides and Sulphides of the Transition Metals 


A recent paper by Morin (1961) classifies these materials into metals and 
non-metals according to the degree of overlap between the 3d functions. 
He considers not only that the direct overlap between 3d functions with e, 
symmetry is important, but also that partly filled orbitals are formed in 
which electrons are shared between these and the 2p — z orbitals of oxygen, 
and that overlap between these may have a comparable effect. Figure 8 
shows his division into metals and non-metals plotted according to his 
estimate of the overlap integrals of the two types. Some of the oxides listed 
as metals only show metallic conductivity at reasonably high temperatures, 
there being a sharp transition to the metallic state as the temperature is 
raised with arise in conductivity which may be as large as 10° (Morin 1959 a). 
Figure 9, taken from Morin’s paper. shows the behaviour of the vanadium 
and titanium oxides. We consider that such oxides must lie near the 
transition line in fig. 7, e.g. at A, and that the transition occurs at 7’, where 
the isotherm (dotted line) crosses the dividing line between the two phases. 
Essentially this explanation is given by Morin (1959 b), though he expresses 
it in terms of the splitting of a band by an antiferromagnetic arrangement of 
spins, which is Slater’s (1951) model. 

There are, however, difficulties in too simple an interpretation of the 
transition. There is first the low conductivity in the ‘metallic’ region, 
smaller by about 1000 than in true metals. The conductivity of a metal 
should be given by 

o=Nne*T / Mett 
where 7 is the time of relaxation and Merz the effective mass. If we write 
+ =1/v where / is the mean free path and v the velocity at the Fermi surface, 
then 
o=ne?l/Mertv. 
Now merrv/h, the wave vector at the Fermi surface, must be of the same 
order as for a metal; the low conductivity can thus only be explained by the 
assumption that /is of the order of the lattice parameter. Several tentative 
suggestions may be made about this. The coupling between the electrons 
and an ionic lattice may be greater than in normal metals, though we should 
be surprised if this made a very large difference to the mean free path. 
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There may be localized spins of the kind suggested. by Friedel e¢ al. (1961), 
which would give a high electrical resistivity nearly independent of tem- 
perature (de Gennes and Friedel 1958). Perhaps this is the most probable 
explanation. Morin (1959b) ascribes it simply to a very narrow band. 
We cannot accept this, because we do not think there could be a transition 
to the metallic state under these conditions. 


Fig. 8 
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Comparison of the overlap integrals S(de, de) and S(de, pz) for a number of 
compounds (Morin 1960). The dotted line is an arbitrary one which, 
in fact, separates metals from insulators. 


Then we must ask how near the point A in fig. 7 must be to B. We suppose 
that all the energies Hy, Hy,, and Ey are likely to be of order me*/2h2x2, 
Unless « is particularly large, and there is no reason to believe it to be for an 
oxide, this should be of the order of an electron volt. Ifso, (V—V,)/V, 
must be of the order 10-3, It seems unlikely that a transition would be 
observed for a number of oxides if A had to be as near to Bas this. On the 
other hand, for VO, a phase change has been observed at the transition 


Transition to the Metallic State 299 


point, (Paoletti and Pickart 1959) and we suggest that the forces holding 
the metal ions in position in the lattice of oxide ions may change somewhat 
at the transition point and that there may be a considerable vhange in 
vibrational entropy. If the vibrational frequency is of order k7'/h, and 
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Conductivity of certain oxides as a function of temperature (Morin 1959 a). 


changes by a factor of order 2, we expect an entropy change of order k 
and thus a free energy change of order k7’. Under these conditions 
(V—V,)/V_ should be of order a few per cent. A pressure of 100000 
atmospheres should, if this is the case, lower the transition temperature bya 


factor of order 2. 
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7.2. Highly-doped Semi-conductors 

The rapid drop in the resistance of most extrinsic semi-conductors that 
occurs when the concentration of impurity exceeds a certain value is well 
known, and has been described by the author (1956) as a transition to the . 
metallic state. Above a certain concentration of impurity, it is found 

(i) that the Hall constant indicates that the number of free electrons is 
equal to the number of centres ; 

(ii) that the conductivity at low temperatures is independent of the 
temperature, and tends to a finite value as 7' tends to zero. The mean free 
path, though of the order of the distance between centres, is roughly 
independent of concentration (Mott and Twose 1961) ; 

(iii) that the magnetic susceptibility, at any rate of n-type germanium, is 
that expected of a degenerate electron gas, the carriers having the effective 
mass appropriate to the semi-conductor. In particular, if V is the number 
of centres per unit volume, the susceptibility is proportional to N13 and is 
independent of temperature (Bowers 1957, Bowers and Yafet 1958, 1960, 
Yafet 1959). 

As a means of investigating the transition to the metallic state, doped 
semi-conductors suffer from two disadvantages. 

(a) The centres are not ina crystalline array, but in general are completely 
disordered (far more so than the metal atoms in a liquid metal). 

(6) At any rate in germanium and silicon and probably in other semi- 
conductors, when the mean distance b between centres exceeds the assumed 
critical value, one does not find a simple state of affairs where weakly 
interacting one-electron centres allow conductivity only with an activation 
energy J—H (which would be near to the energy J required to excite an 
electron into the conduction band). This is because the position is compli- 
cated by the presence of minority centres (acceptors in n-type conductors). 
These will receive electrons, and the resulting unoccupied donors can, by 
tunnel effect, receive electrons from neighbouring donors, thus allowing a 
current to pass. However as pointed out first by the present author 
(1956; see also Price 1957, Miller and Abrahams 1960 and Mott and Twose 
1961), this impurity conduction in compensated semi-conductors cannot be 
independent of temperature. As the temperature 7 tends to zero, the 
conductivity will tend to zero as exp (— W/kT) where W is an activation 
energy due to the charged minority centres. Mott (1956) first considered 
the case when the compensation ratio K , defined by 


K=N min/ N maj 
is very small; here Nmin and NV maj are the concentrations of minority and 
majority centres. He stated that in this case one should have 
W=e|cR 
where & is the mean distance from a minority to the nearest Majority 


centre. W is then simply the energy required to remove the ‘hole’ on a 
majority centre to an infinite distance from the charged minority centre. 
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As K increases, W decreases rapidly as predicted by Price and by Miller and 
Abrahams, and observed by Fritzsche and Cuevas (1960). 

The transition is fairly sharp, as shown in fig. 10. To the right of the drop 
the conductivity is independent of temperature. 


Fig. 10 
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We expect, according to the crude estimate of § 3, that the critical density 
N for semiconductors that are elements should be given by 


N43 a, = 0-25 


where 
Ay = Kh? |Mers e?. 


Fritzsche and Cuevas (1960) have shown that the critical concentration 
depends little on compensation between values of K of 0-04 and 0-4, not 
by more than 10%. An exact experimental value of V1? a, is difficult to 
obtain because the wave functions of the centres are not hydrogen-like ; 
approximately the results for Ge and Si suggest a value of about 0-2 
(Mott and Twose 1961). 

For semi-conductors which are not elements it is not quite clear what 
value of the dielectric constant should be taken. In the metallic state the 
force between the electrons will be e?/« yr? where x, is the high-frequency 
dielectric constant, but in the non-metallic state the charged centres 
however will polarize the ionic lattice round them. Further investigation 
of this problem will be necessaryt. 


7.3. Heavily Compensated Semi-conductors and Wigner’s Hypothesis 


If K, the compensation constant, is nearly equal to unity, an interesting 
case may arise which would enable one to investigate whether Wigner’s 
assumption that a free electron gas can ‘crystallize’ is true or not. One 
will then have many more majority centres than electrons. Ifthe distance 
between majority centres is large enough, then (Anderson 1958) the states 
for these electrons are localized. In other words, the wave functions for any 
one of these states decays exponentially with distance from one of the 
centres. As V, increases, however, and overlap between the wave functions 
of the centres becomes strong, the states of a single electron become 
delocalized, the wave functions extending throughout the lattice. We thus 
have a gas of n( = Ny —N,) electrons per unit volume moving in the field of © 
arandom distribution of positive and negative point charges, the mean 
distance between them being small compared with n-3, This seems to us 
the nearest approach available experimentally to the original Sommerfeld 
model of a metal, a gas of electrons moving in a uniform distribution of 
positive charge (‘jellium’). 

If Wigner (1938) is right, then at the absolute zero of temperature these 
electrons ought to crystallize, say in a body-centred cubic lattice, if the 
quantity na, is greater than some critical valuet. We shall make a 
very rough estimate of this critical value. We denote by 7, the radius such 
that 47nr,3/3=1; then in the non-metallic state each electron is thought to 


em ee ee 
a Compare McIrvine (1960) who has given an empirical formula relating N to 
the dielectric constant for several materials, 
{ For intermediate values it is always possible that an ordered antiferro- 
magnetic arrangement of spins would arise, as envisaged by Friedel et al. (1961). 
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move in the field of a sphere of this radius filled with a uniform distribution 
of positive charge, and thus with a potential energy 


where r is the distance from the centre of the sphere. Wigner ascribes to 
each electron a fixed potential energy —4e2/r, and a zero-point energy 
3hw/2, where 
w= /(e?/mrp). 

This will be correct for large interatomic distance, but near the transition 
point it may be better to set for the potential energy the average over the 
sphere 

To 

| V(r) 4rr? dr/(477,3/3) = — 2e?/5ro, 

0 
and for the zero point energy that of a particle in a spherical cavity of 
radius ry, namely 7?h?/2mr,?. 

For the metallic phase, since each sphere is occupied by all the electrons, 
the Hartree approximation will give half the above contribution to the 
potential energy, namely 

— 2e7/107. 


The zero point.energy is now the Fermi energy, namely 


SA Peri 8 Va 
5 2mry? \ 2a 


The transition will thus occur when 


Le? Cat 
a 3 
oo 2mro 
or 
1) ~ 20 aq. 


While this paper was in proof, Dr. Carr showed the author the manuscript 
of a more detailed calculation along these lines (Carr et al. 1961). They 
find energy curves cutting when ro/a, ~~ 3, though they appear to assume, 
contrary to our belief, that there is a continuous transition from one form 
to the other. 
If the phenomenon is to be observed, the density of centres must be big 
enough for bound states to disappear. According to an extension of 
Anderson’s (1958) work due to Twose (Mott and Twose 1961) this should 
occur when Ra, is of the order 3-6, where N =3/47K® is the density of 
majority centres. In germanium this would mean that the density of 
centres must be greater than about 10!’ cm~*. ; 
Actually it may be doubted whether a crystallization with long distance 
order could ever beobserved. The energy of melting might be expected to 
be of the order of 1/100 of e?/«ryand thus only a few hundredths of a degree in 
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germanium. Also the random field due to the charged centres might lead 
to a polycrystalline arrangement, even at the absolute zero. But a sharp 
transition between a normal liquid at low concentrations and the condensed 
electron gas at high concentration might occur. Perhaps the most straight- 
forward way to observe any such effect would be a change in the magnetic 
susceptibility. 

The transition can be compared to that between liquid helium I and IT; 
the fact that the repulsion between helium atoms falls off as a higher power 
of the distance than the zero point energy and a lower power for electrons 
means that the super fluid (condensed gas) is destroyed by compression in 
the former case, by expansion in the latter. 


FLAT-BAND POTENTIAL CVS CALOMEL) 


ELECTRON DENSITY 


The flat-band potential as a function of electron density for crystals freshly 
etched in H;PO, (Dewald 1960). 


7.4. Fermi Energy in Zinc Oxide Activated by Indium 


Dewald (1960) has measured the capacity of the interface between an 
electrolyte and zinc oxide activated by indium, and from these measure- 
ments he deduces what he calls the ‘flat band potential’, i.e. the potential 
at which there is no Schottky barrier in the surface layer. This gives a 
measure of the Fermi energy of the electrons in the semi-conductor for 
various concentrations of the indium, measured against a fixed level 
determined by the concentration of the electrolyte. His results are repro- 
duced in fig. 10, which shows—with sign reversed—the Fermi energy as a 
function of concentration. At about the concentration marked A metallic 
conductivity (i.e. conductivity without temperature-dependence) is known 
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to set in (Thomas 1958). We make the hypothesis that the discontinuity 
in the Fermi energy at A, its drop to a minimum at B and subsequent rise is 
due to the same cause as the corresponding behaviour of the total free 
energy shown in fig. 6. The minimum at B, which should not be sharp if 
this explanation is correct, should be the point where the free energy of the 
electron gas isa minimum. The curve suggests that the transition to the 
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Solubility gap (a) and limit of solubility (6) for sodium in ammonia, as & function 
of concentration c of ammonia (Bingel 1953). Temperatures in °c. 


metallic state occurs when the electron density is less by a factor 2 than at 

this minimum. If 7,/7, is of order 3 for the transition to the metallic state 

as for germanium, this would suggest that the minimum occurred for 

rolty~ 2-4. This is considerably larger than the values calculated ig 

metallic hydrogen by Wigner and Huntingdon (1935) and by Stern an 

Talley (1955), (1-5-1-6), but since the electron in indium is in a p state 
there is no reason to think these results applicable. 


U 
P.M. 
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7.5. Metal_-Ammonia Solutions 


If the centres were free to move about, fig. 6 shows that one would 
expect a two phase region} between the points marked P and Q. ‘This is 
observed for solutions of sodium in liquid ammonia (Kraus 1953, Bingel 
1953), for which the phase diagram is as in fig. 12. Solutions of about 10% 
sodium in ammonia are probably ‘metals’ of the type discussed here, the 
sodium atoms being expanded by the surrounding dielectric constant (its 
high-frequency value). The free energy of the solutions, including the 
terms representing the cohesive energy between these centres, must be ata 
minimum about here. 


§ 8. LocatizeD States IN METALS 


Several authors recently (e.g. Goodenough 1960) have referred to localized 
states in transition metals and the present author (Mott and Stevens 1957) 
and also Lomer and Marshall (1958) make use of this concept in discussing 
the magnetic electrons in iron. In this section we shall not discuss 
individual metals in detail, but shall state certain general principles. 

We have emphasized that, for an array of one-electron atoms, impurity 
states in a semi-conductor or inner shells of a transition metal oxide, the 
question whether the material is ‘ metallic’ or not is one capable of a precise 
answer. Whether it is a conductor or insulator at the absolute zero of 
temperature is a physical question and not a matter of the most convenient 
theoretical model. For inner shells of transition metals, the question can 

-only be answered to the approximation to which one writes down separate 
many-electron wave functions for the inner and outer electrons. Assoonas 
one takes into account the interaction between them, the question is no 
longer precise; one can only ask which is the best model. Nonetheless, 
the question may be worth asking. 

Let us then consider a metal in which the atoms have an incomplete 
inner shell, for instance the rare earth metals. The overlap between the 
4f shells is so small that it is most unlikely that a ‘metallic’ type of wave 
function would be appropriate or that the 4f shells contribute to the 
Fermi surface. On the other hand, one will always get a metallic type of 
wave function if the number of electrons in the inner shell is non-integral 
(as in nickel). Such a state of affairs is possible in, for instance, a rare 
earth metal. Suppose an extra electron is added to one of the 4f shells. 
The energies that this electron can have will lie in a narrow band. If the 
Fermi energy of the conduction electrons lies in this band, it will fill up until 
the Fermi energy of the narrow 4f band and of the conduction band are the 
same. Ifthe Fermi energy of the conduction band lies below it, the wave 
functions of the 4f electrons will remain non-conducting. If it lies suffici- 
ently far below it, however, it may overlap the band obtained by taking an 
electron out of one of the 4f shells. 


} Lam indebted to Dr. Dewald for pointing this out to me. 
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If then one were to plot the number N of 4f electrons against the Fermi 
energy of the conduction electrons, a curve similar to fig. 13 would be 
expected, being an integer. In the regions where JN is non-integral, one 
has a band of energy levels and a Fermi surface. A band of energy levels. 


Fig. 13 


Ferm: Energy 


Number of electrons expected in an inner band as a function of the Fermi energy 
of a metal. 


Fig. 14 


. 


Suggested splitting of an inner band into two (or more) sub-bands when the 
overlap between atomic orbitals is small. 


containing only one electron per atom isanew concept. If however forn and 
for n +1 electrons the states are non-conducting, it must be a valid one; 
and s band, for instance, that would normally have two electrons will split 


into two bands, a p band into six (fig. 14). 
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We may also emphasize that it is only in a metal that we may anticipate 
that, if a band is narrow enough to have a non-metallic wave function for 
an integral number of electrons, it will necessarily have a Fermi surface at 
T =0 when the number is non-integral. For non-metals the position is 
made clear by the behaviour of compensated semi-conductors (§ 7.2) ; at 
very low temperatures the carriers are held in position by the corresponding 
positive charges, and the material is not a conductor. 
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CORRESPONDENCE 


Hammer Tracks from the Photodisintegration 
of Light Emulsion Nuclei 


By W. T. Morton and T. G. Waker 
Department of Natural Philosophy, University of Glasgow 


[Received September 14, 1960] 


Ir *Li or *B is produced in a nuclear disintegration it will decay at rest 
into two alpha particles through an excited state of *Be, and so produce 
a characteristic T-shaped track. These are usually referred to as 
‘hammer tracks * (Occhialini and Powell 1947, Alvarez 1950). Events 
of this type in which the initial disintegration was produced by a particle 
have been investigated by a number of workers (for example, Titterton 
1951). In the present investigation a search has been made for hammer 
tracks originating in the photodisintegration of the light nuclei (carbon, 
nitrogen and oxygen) of nuclear research emulsions. 

Ilford C.2 nuclear research emulsions, 400 microns in thickness, were 
exposed perpendicular to a bremsstrahlung beam of maximum energy 
120 Mev from the Glasgow electron synchrotron. The emulsions were 
scanned for nuclear disintegrations in which more than two charged 
particles were emitted. 4300 events of this type were examined and 
eight hammer tracks identified. A small number of the three-prong 
events could have been due to the emission of a ®Li or §B particle from 
a heavy nucleus of the emulsion (silver or bromide) for which no other 
charged particle was emitted. These events might be confused with a 
low energy 3« photodisintegration of #C and so were neglected. In 
each of the accepted hammer track events the initial disintegration 
had also produced a singly charged, or doubly charged particle of energy 
less than 5 Mev, and it is considered unlikely that these would be 
emitted simultaneously with a hammer fragment from a heavy nucleus 
of the emulsion. Hence it seems reasonable to assume that all of the 
hammer tracks accepted originated in the light nuclei of the emulsion. 
In the investigation of hammer tracks emitted in reactions induced by 
neutrons with an energy of 150 Mev, Titterton (loc. cit.) reached a similar 
conclusion. It would be expected that the probability of emission of a 
hammer fragment from silver or bromine would decrease with energy 
of the incident particle. 

If a 8B particle was emitted from a carbon nucleus, by the conservation 
of charge, a singly charged particle would also be emitted. Hence, 
although the §B would be recorded as a three-prong event, its true nature 
would be indicated by the sharp change in grain density and possibly 
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also in direction at the position of the initial disintegration. No events 
of this type were found. This result is in agreement with Titterton who 
found that the probability of a hammer track resulting from ®Li is higher 
than that for °B. In the remainder of this communication all hammer 
tracks will be assumed to be due to SLi. 

No events were observed which could be assigned to °Li, which decays 
into *Be* followed by a decay at rest into a neutron and two alpha 
particles (Fry 1953). In this case the two alpha particles would not be 
collinear. 

The energies of the *Li tracks were found to be 2-2, 3-2, 3-2, 5-1, 6-5, 
8-3, 10-5 and 14-5 Mev. Three of the events in which the °Li particle 
originated had three prongs, one of which was singly and the other 
doubly charged. Hence by charge conservation these events are produced 
by the reaction 'C(y, *He p)8Li with a threshold value of 41-2 Mev. 
The remaining events had four prongs, three of which were due to singly 
charged particles. These events must also have originated in a carbon 
nucleus of the emulsion. Since it was not possible to differentiate between 
protons and deuterons these events were produced either by the reaction 
Cy, 3pn)$Li with a threshold value of 50-9 Mev or the reaction 
LC(y, 2pd)§Li with a threshold value of 48-7 ey. 

To obtain the incident y-ray flux the three-prong events were analysed, 
using the criteria of Goward and Wilkins (1952), for events produced 
by the photodisintegration of 12C into three alpha particles. Using the 
cross section given by these authors for that reaction the cross section 
for the photoproduction of SLi from "C averaged over y-ray energies 
between 50 mev and 120 Mev was found to be (4:9 + 2-0) x 10-89 cm?, 
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A Comparison of Vacancy and Interstitial Loops in Graphite 


By G. K. WitLramson 


Berkeley Nuclear Research Laboratories, Gloucestershire 


and C. BAKER 
Department of Metallurgy, Cambridge University 


[Received September 22, 1960] 


A recent letter by Amelinckx and Delavignette (1960) suggested the 
possibility of distinguishing between interstitial and vacancy loops in 
graphite because the vacancy loop has a component of its Burgers vector 
in the basal plane whereas the interstitial loop has a Burgers vector 
parallel to the C axis. Independent work carried out in the Department 
of Metallurgy, Cambridge, has proved that this is so and has also produced 
positive evidence for Amelinckx’s interpretation of the Burgers vector 
of the vacancy loops. 


Fig. 4 


Possible interaction of a loop and a gliding basal plane dislocation to produce 
the dislocation arrangement shown at F in figs. 1-3 (arrows indicate 
Burgers vectors). 


The photographs shown in figs. 1, 2 and 3} are of quenched and 
annealed graphite containing vacancy loops. The graphite was heated 
by direct current in a vacuum, and quenched by switching off. Figures 2 
and 3 are dark field photographs of the area shown in fig. 1, taken with 
1210 and 1120 reflections. The partial Burgers vector of the type 


+ Figures 1-3, 5 and 6 are shown as plates. 
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a/3[1010] will lie along one of the three projections of the {1120} planes 
on the basal plane and so must go out of contrast on one and only one 
of the three different (1120) dark field views (Williamson 1960 and 
Delavignette and Amelinckx 1960). It is necessary to observe the absence 
of contrast to prove that the component of the Burgers vector in the basal 
plane is in fact a partial. This does in fact occur for every loop. In the 
example shown loop A disappears on fig. 2 and loop B on fig. 3. The 
large loop C at first sight appears on all three dark fields, but this is due 
to the contrast of dislocations held to its periphery by long-range stress 
fields, the segment where there are no overlapping dislocations does go 
out of contrast on fig. 3. Because these loops contain a component in 
the basal plane, interaction is expected with basal plane dislocations and 
one type of interaction is shown at F. Figures 1-3 show that the loop 
has different partial Burgers vectors at opposite sides and a possible 
mechanism for producing this effect is shown in fig. 4. 

Examples of large interstitial loops produced by neutron irradiation 
(~3x 10” nvt) at 300°C are shown in fig. 5. In the case of interstitials 
a disc of atoms can be inserted as a C layer in an ABAB ete. sequence 
without producing a high-energy stacking fault. Thus the basal plane 
component which makes the vacancy loops so easy to observe is missing. 
Interstitial loops rarely show good contrast and that shown in fig. 5 
is produced by tilting the specimen so that the C axis is about 30° away 
from the electron beam. The contrast effects in this case are clearly 
very complex and a detailed study has been hindered by the fine scale 
radiation damage clearly visible in the broad extinction contours. The 
basal plane dislocations can interact with these loops through their 
long-range stress fields as in fig. 5. Figure 6 shows the contrast of 
interstitial loops when the C axis is parallel to the electron beam. This 
contrast has been discussed in detail by Howie and Whelan (1960) who 
show that the loops are visible because of the radial displacements 
perpendicular to the Burgers vector. Under these conditions two 
concentric loops on either side of the dislocation loop are seen with 
contrast most intense at those parts of the loop parallel to the reflecting 
plane and vanishing at points perpendicular to it. This is exactly the 
contrast observed and in all cases examined the diameter along which 
there is no contrast is parallel to the reflecting plane. 

A fuller account of this work will be published elsewhere. 
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Evidence of Sub-surface Cross-slip in Fatigued Copper 


By A. H. Metexa, W. Barr and A. A. BAKER 
Metal Physics Section, British Iron and Steel Research Association, London 


[Received October 5, 1960] 


Ir is now generally accepted that the formation of extrusions and 
intrusions on the free surface of fatigued metals leads to the initiation 
of the ultimate fatigue crack. Particular attention has therefore been 
given to the possible atomic movements during cyclic stressing responsible 
for extrusion-intrusion formations. Various mechanisms have been 
suggested including those due to Cottrell and Hull (1957), and Forsyth 
(1957) but the one proposed by Mott (1958) appears the most likely. 
Briefly, it is suggested that during the first half-eycle a dislocation line 
moves along the slip plane and then cross-slips to an adjacent plane. 
In the second half-cycle the dislocation glides backwards, then cross-slips 
to the original plane, and so on. In the presence of a screw component 
of the dislocation this continuous movement will lead to the formation 
of an extrusion above the free surface and a cavity below the extrusion. 
Cross-slip is an essential feature in this mechanism. Although the 
importance of cross-slip to the phenomenon of fatigue is now generally 
accepted by most workers in this field, Mott’s mechanism has not, as 
yet, been confirmed by direct observation. We have observed, by using 
a fine electropolishing technique and by penetrating progressively into 
the surface of fatigued copper specimens, the existence of cross-slip in 
those slip bands which eventually lead to crack formation. 

Specimens of OFHC copper were fatigued in reversed bending to about 
half their expected life. The wide slip bands characteristic of fatigue 
deformation were observed. Upon electropolishing, using the technique 
developed by Jacquet (1957), most of the surface markings were removed, 
leaving behind a number of dark bands, the so-called persistent slip 
bands. The present technique is superior here to standard electro- 
polishing methods for two reasons. Firstly, very thin layers, estimated 
at 100A, can be readily removed; this is particularly useful since early 
fatigue damage is concentrated in a thin surface layer. Secondly, the 
persistent slip bands produced by this technique are considerably finer 
than those obtained by standard electropolishing ‘methods. Persistent 
slip bands can, therefore, be observed at higher magnifications and fine 
structural details are revealed. ; 

Figure 1(a)+ shows the appearance of a grain after removal of 
approximately 3 microns ; persistent slip bands on two slip systems can 
be seen. Successive layers were then removed and fig. 1 (6) shows the 


+ Figures are shown as plates. 
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same field of view after the removal of a further 9 microns. The few 
persistent slip bands that remain, differ from those in fig. 1 (a) in that 
they deviate from linearity by the presence of short segments inclined 
to the main slip plane. In the region to the left of fig. 1 (b), for example, 
the segment indicated by the arrow A coincides with the direction of the 
conjugate slip plane, while that indicated by B coincides with the 
cross-slip plane. This was confirmed by constructing a stereographic 
projection based on the sideway shift of the primary slip bands produced 
by the removal of 9 microns. Other manifestations of cross-slip are 
indicated and more can be seen elsewhere in the micrograph. 

It is interesting to note that cross- and conjugate-slip segments are 
more in evidence upon deeper penetration into the surface layer. It is 
likely, as suggested by McEvily and Machlin (1959) that cross-slip 
provides the means for the easy transport of piled-up material beneath 
the intrusion to the region beneath the extrusion. 

It was further observed that, in certain cases, cross-slip segments 
were present below the extrusions. Examples of extrusions and the 
corresponding cross-slip regions are shown in fig. 2 (a) and (6) respectively. 
Although this one to one correspondence was observed in a few cases, 
it is by no means general. Numerous extrusions were observed which 
could not be related to regions of cross-slip. 

Further examples of the abundance of cross-slip in fatigued copper, 
together with a full account of this work, will be published elsewhere. 
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REVIEWS OF BOOKS 


Graphite and its Crystal Compounds. By A. R. UspeLonpe and F. A. Lewis. 
(Oxford: Clarendon Press, 1960.) [Pp. xii+217.] Price £1 15s. Od. 


THE use of graphite as a moderator in nuclear reactors has given an impetus to 
the study of this most interesting and versatile material. It is now clear that 
the properties of graphite are so diverse as to require a monograph for their 
description. Professor Ubbelohde and Dr. Lewis have now provided such an 
account. Many of the properties of graphite, such as its increased electrical 
conductivity when electrons are either added to it, or subtracted from it, depend 
on its almost unique electronic band structure. Other properties, such as its 
layer form with normal bond distances within a layer and much larger distances 
between layers, link it to the peculiar chemistry of aromatic molecules with 
theirz and celectrons. Yet other properties, such as the bizarre way in which a 
wide variety of atoms and ions, including such diverse entities as alkali metal 
atoms, halogen atoms and AICI, moleclues are able to prize apart the layers and 
intercalate themselves between them, thus increasing the interlayer distance 
from 3-5 A to as much as 11 A, depend on the relative weakness of van der 
Waals forces as compared with valence forces. Much of the reactivity of 
carbons depends critically on its purity, and so the nature of the residual 
foreign atoms around the edge is of prime importance. So also are the lattice 
defects—a matter only recently recognized. All this—and much more—is set 
out in this book. The explanations are clear, and (almost always) consistent. 
This, incidentally, is by way of being a minor triumph, since the range of altern- 
ative theories suggested from time to time is prodigious. In places the attempt 
to refer to almost every published paper on graphite gives the book the character 
of a card index, but anyone working in this field is sure to be grateful for the 
bibliography of no less than 28 pages. As a reference work on graphite, with 
all its vagaries and still largely uncomprehended behaviour, this book is unique. 
It is easy to read and beautifully produced. C. A.C, 


Frozen Free Radicals. By G. J. Miyxorr. (New York: Interscience Pub- 
lishers, 1960.) [Pp. 134.] Price $5.00 (£1 18s. Od.). 


THE possibility of freezing or trapping high concentrations of free radicals has 
been very much in the scientific news over the last year or so, as illustrated by 
the large-scale research programme recently undertaken by the National Bureau 
of Standards in Washington on this very subject. Dr. Minkoff was one of the 
Guest Scientists associated with this project and one suspects that much of this 
enthusiasm, and many of the ideas developed in the writing of this book, stem 
from this association. The book has the flavour of a subject that is full of 
interest and possible practical applications, and is written by an author who is 
not only an active worker in the field itself, but also has obviously taken a keen 
interest in all the other lines of development. When considered in ten years’ 
time some parts of the book will undoubtedly appear to have undue emphasis, 
and others may need modification or correction, but the fact that it is impossible 
to tell which sections these are at the moment, shows that the monograph is a 
fair summary of the present position in the development of the subject. 

On the other hand, it should not be understood from these remarks that the 
book only deals with present day work and the more recent advances in free- 
radical studies. After a short introduction to the subject, a whole chapter is 
devoted to the early history of free-radical investigations, and forms a very 
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interesting and helpful summary, which is not readily available in many 
similar monographs. The more recent work is then introduced by a considera- 
tion of the experimental techniques, first grouped under the methods of pro- 
duction of radicals, and then the methods of trapping and stabilizing at low 
temperatures, and finally the methods of detecting and investigating the frozen 
radicals. Owing to the space available, each of these topics is considered in a 
rather condensed fashion, but considerable detail is given on the design of low 
temperature Dewars and cells, including several detailed diagrams, which will 
prove very useful to any who are entering this field for the first time. 

A chapter on the basic principles of the ‘ trapping mechanism ’ then follows, 
together with a quantitative treatment of both trapping and chain processes. 
The rest of the book is then devoted to a summary of the results obtained with 
different types of radicals. This starts with a chapter devoted to the study 
of individual atoms; hydrogen, oxygen and nitrogen being considered in some 
detail. The different experimental results are summarized for each case, and it 
is interesting to see how large a part electron resonance spectra play in the 
identification of the atomic species. The following chapter deals with diatomic 
and triatomic radicals in a similar way, while the concluding section summarizes 
the results on complex organic radicals of both aromatic and aliphatic nature. 

The book can certainly be taken as a fair and comprehensive treatment of 
the whole field of ‘ frozen free radicals ’ as it exists at the moment, although a 
second edition will undoubtedly be required before too long. There are one or 
two minor critisisms of careless phrasing in the book, such as the statement that 
“ scattering is less important at shorter frequencies ”, when it is not immediately 
evident whether this means ‘shorter wavelengths’ or ‘smaller frequencies ’. 
These hardly detract from the main content of the book, however. 

Ds J. EI, 


Beam and Wave Electronics in Microwave Tubes. By Rupotr G. E. Hurrer. 
(New York: D. van Nostrand, 1960.) [Pp. 378.] £3 13s. 6d. 


In the last decade it has become clear that the best description of the operation 
of microwave valves depends on the application of Maxwellian electrodynamics 
to the highly complicated system comprising the electron beam in interaction 
with some kind of slow-wave structure. This model has been the subject of 
hundreds of articles and several books in several languages. The work under 
review is an extension of an article in Advances in Electronics which appeared 
in 1954. The new material includes chapters on circuits (about one-third 
of the book), extensive quotation from the work of Gould on electron waves 
in crossed electric and magnetic fields, an account of Pierce’s mode coupling 
concept and a chapter on noise. The preface states that the level is that of 
a university textbook. 

The author has decided to ignore many well-known physical facts. For 
example, most microwave valves in use today use a quasi-cylindrical electron 
beam formed by some type of magnetic focusing. The trajectories of the 
electrons depend on the focusing technique and the space-charge wavelength, 
radio-frequency currents, etc., depend on the trajectories in a very marked 
way. The behaviour of a beam using Brillouin focusing is thus quite different 
from that of a confined beam. Yet none of this is mentioned, nor are references 
given. Again, it is well known that even the simplest case of interaction between 
a cylindrical beam and a circuit, in which the beam just fills the gaps of a 
cavity, has to be described using a complete set of orthogonal modes on the 
beam. The existence of modes of higher order is never mentioned and, as a 
result, the basic work in Chapters 8 and-9 is never put on a firm theoretical 
basis. _For the same reason, the discussion of energy conversion becomes 
unrealistic immediately it goes beyond the derivation of an appropriately 
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modified Poynting theorem to consider effects in real tubes. The delicate 
question of the best description of phenomena on the surface of a rippled 
beam has recently been discussed, but this is not referred to. 

While the author could be excused for not including full treatments of 
these questions, to omit them entirely is inexcusable and the book cannot 
be recommended. AN SW 3B. 


Electromagnetic Wave Propagation. Edited by M. Diistrant and J. L. Mrcurets. 
(London and New York: Academic Press, 1960.) [Pp. xiv+730.] $22.00. 


THIs is a collection of 54 papers which were presented in 1958 at a conference 
in Brussels sponsored by the Postal and Telecommunications Group of the 
Brussels Universal Exhibition. Five of the papers are in German, four are 
in French and the rest are in English. Nearly all of them relate to some aspect 
of radio wave propagation. 

About 20 of the papers are theoretical and the topics include diffraction 
theory, antenna radiation theory, surface waves, ray tracing in the ionosphere, 
propagation of waves of low and very low frequency, scattering, and the 
statistics of irregular media. About 27 papers describe measurements or 
observations, including ground wave propagation, tropospheric scatter, 
diffraction beyond the horizon, ionospheric backscatter, scattering from 
meteor trails, the scintillation of radio stars, the reflection of low frequencies 
from the ionosphere, and ‘round the world’ echoes. The remaining papers 
are concerned with analysis of data and with radio engineering techniques. 

The contributors include many well known investigators from research 
centres all over the world. All who are accustomed to watch the current 
literature on radio wave propagation will need to make a careful study of the 
contents list of this book. Ke Gab: 


An Introduction to the Physics of Mass Length and Time. By NoRMAN FEATHER. 
(Edinburgh University Press, 1959.) [Pp. ix+358.] 18s. 


PROFESSOR FEATHER’S aim in this book is “to introduce the beginning 
specialist in physics to those parts of his subject which in any well-desighed 
course he should study first, and to effect that introduction by examining 
first of all the bases of measurement of the fundamental physical quantities 
mass, length and time’”’. The term ‘ beginning specialist’ is perhaps mis- 
leading, for Professor Feather is writing not only for the future scientist but 
for all students, whether scientists or not, who wish to acquire a worthwhile 
grasp of some of the elements of physics, and this of course influences his 
presentation ; in particular, the calculus is not introduced. 

After an introductory chapter the author discusses the measurement of 
length (in both theory and practice) and of time, and he is then in a position 
to discuss a number of topics in kinematics, including Kepler’s laws and 
Newton’s treatment of them. A key chapter on the concepts of force and 
mass follows, after which rather more than half the book remains in which to 
make the work “‘ not very different in scope.from any elementary Mechanics 
and Properties of Matter’. 

An historian of science cannot but welcome the frequent efforts to set concepts 
and laws into their historical context, and here Professor Feather is usually 
well informed. Most unfortunately he does not recognize that the very attitude 
to physics which he accepts almost without question was itself the product of 
intellectual struggles far more profound that those he has occasion to discuss. 
Physics “may be defined as the body of knowledge derived by making 
measurements of things ’—what pangs attended the birth of that definition !; 
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and did not the adoption of idealized mathematical quantities constitute a 
paradoxical method of studying nature? For an understanding of physics 
such as Professor Feather aims to give, the reader must be made aware of these 
fundamental historical issues; but it is good to see this attempt to teach 
physics as a great human achievement. M. A. H. 


Some Mathematical Methods of Physics. By Guratp GorrTzEL and Nunzio 
Trait. (McGraw-Hill Book Co., 1960.) [Pp. xiv+300.] £3 6s. Od. 


Tuts book shows how to apply the calculus of operators to problems of 
theoretical physics. Much of the mathematical technique used in quantum 
theory is presented here, including the Dirac bracket notation. Yet there is 
hardly any mention of quantum theory, and all the problems discussed are 
classical. 

The book is in three parts. Part I deals with systems having a finite number 
of degrees of freedom, such as the loaded string, and the operators are repre- 
sented by matrices. Part II is concerned with systems having an infinite 
number of degrees of freedom. The operators are differential operators with 
boundary conditions, and the problems used for illustration include the 
vibrating membrane, electrical lines, heat conduction and problems of 
radiation and scattering. Part III discusses methods of approximation, 
including perturbation methods and the Rayleigh—Ritz variational principle. 

Most students taking a final honours physics course will find this book 
difficult to understand. It should be more useful to teachers already familiar 
with the main theme, who want ways of presenting the ideas clearly to students. 
Readers who master it will be well equipped to tackle an advanced course on 
mathematical quantum theory. BO Gen. 


The Theory of Brillouin Zones and Electronic States in Crystals. By H. Jonzs. 
(North-Holland Publishing Company, 1960.) [Pp. 268.] £3 Os. Od. 


THE bulk of the book is devoted to a thorough discussion of how to specify 
states in a periodic lattice including topics such as character tables, compati- 
bility relations, screw axes and glide planes, spin-orbit coupling, but only a 
very incomplete treatment of time-reversal symmetry. The necessary group 
theory is introduced very well, almost imperceptibly and without tears, in 
the form directly required by the applications without mathematical frills. 
One chapter deals with practical methods of calculating band structures 
discussing three methods in detail. 

Chapter 5 is concerned with the large zones in extended k-space, which are 
often introduced for metals and alloys having a large number of valence 
electrons per unit cell. Unfortunately the analysis is all based on the nearly 
free electron approximation (n.f.e.a.) with weak potentials. Since in practice 
atomic potentials are strong and the n.fe.a., if carried out, would converge 
to the ls energy of the core orbitals, one is left wondering what relevance it 
all has to real metals. This could easily have been avoided by basing the 
analysis on the orthogonalized plane wave method, where the structure factors 
ee in exactly the same way but where the matrix elements really are 
small. 


The book is clearly written throughout, and fills an important place in the 


series of monographs into which solid-state theory has now expanded. V. H. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 
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A. N. GOLAND Phil. Mag. Ser. 8, Vol. 6, PI. 13. 


(0) 
¥lectron micrographs of the two copper samples shown in fig. 2. Magnification 
is x 20250. (a) The normal grain structure of the out-of-pile copper 


sample. (6) The normal grain structure of the in-pile copper sample. 


A. N. GOLAND Phil. Mag. Ser. 8. Vol. 6, Pl. 14. 


A. N. GOLAND Phil. Mag. Ser. 8, Vol. 6, Pl. 15. 


Fig. 4 (continued) 
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Electron micrographs of the precipitation region of the two copper samples. 
Magnification is x 20250. (a) The grain structure of the precipitation 
region in the out-of-pile sample. (6) The grain structure of the 
precipitation region in the in-pile sample, showing an active grain 
boundary and one which is not active in supplying vacancies. (c) The 


grain structure in another portion of the precipitation region in the 
in-pile sample. 


A. R. C. WESTWOOD Phil. Mag. Ser. 8, Vol. 6, Pl. 16. 


Cracks formed at grain boundary by piled-up groups of edge dislocations. 


A. R. C. WESTWOOD 


Crack propagating away from boundary on 
a {110} plane. 


Fig. 4 


Phil. Mag. Ser. 8, Vol. 6, Pl. 17. 


Fig. 6 


\ - 
Ce a 
Crack formed at U by cooperation of dislo. 

cations in band V and from source SN. 


Crack at P possibly formed by interaction of dislocations in slip bands Q and T. 


Phil. Mag. Ser. 8, Vol. 6, Pl. 18. 


A. R. C. WESTWOOD 
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R. W. DAVIDGE and R. W. WHITWORTH Phil. Mag. Ser. 8, Vol. 6, Pl. 19. 
Fig. 1 


Etch pits on cleaved surface of undeformed crystal, showing pits at annealed-in 
dislocations and cleavage steps. (x 260.) 


Fig. 2 


Part of rosette of etch pits on crystal indented, polished for 1 min and etched. 
(x 410.) 


Fig. 3 


Normal and shallow etch pits in an edge slip band. (x 650.) 


R. W. DAVIDGE and R. W. WHITWORTH Phil. Mag. Ser. 8, Vol. 6, Pl. 20. 


Fig. 4 


Pits showing intermittent etching in an edge slip band. (x 580.) 
Fig. 5 


(b) 
Etch pits on edge slip bands before and after annealing. (a) Deformed, cleaved, 
polished and etched. (6) Same region annealed 1 hr at 250°c, polished 
and re-etched. (x 190.) 


R. W. DAVIDGE and R. W. WHITWORTH Phil. Mag. Ser. 8, Vol. 6, Pl. 21. 


Fig. 6 


Ktch pits on part of a screw slip band. 
External surface: (a) oblique illumination, (b) phase-contrast. 
Polished 1 min: (c) oblique illumination, (d) phase-contrast. (x 540.) 


Fig. 7 


Streaks of pits formed by etching freshly cleaved surface of undeformed 
crystal. Phase-contrast.  ( x 390.) 


R. W. DAVIDGE and R. W. WHITWORTH Phil. Mag. Ser. 8, Vol. 6, Pl. 22. 


Fig. 8 


Part of rosette in crystal vibrated at 90 ke/s. (a) Polished, etched, indented, 
vibrated at strain amplitude of 2-6 x 10-4and re-etched. (b) Same region 
polished 3 min and re-etched. (x 230.) 


R. W. DAVIDGE and R. W. WHITWORTH Phil. Mag. Ser. 8, Vol. 6, Pl. 23. 
Rigeo 


Pits associated with edge dislocations in crystal indented, vibrated at strain 
amplitude of 2-0x 10-4, polished 3 min and etched. (x 560.) 


Fig. 10 


(a) 


0 Ey TE yy 


(b) 


Etch pits in slip bands on matched cleavage faces of deformed LiF. (a) Etch W, 


(0) Etch H. Note poorly defin, d pits fi : : 
(diagonal lines). (x530) > Ormed by H in edge slip bands 


F. KROUPA and P. B. PRICE Phil. Mag. Ser. 8, Vol. 6, Pl. 24. 
Fig. 1 


(¢) O-lp 

; ; eee aaa 
Sequence of transmission electron micrographs showing the ‘ conservative 
climb ’ motion of a dislocation loop with a [0001] Burgers vector due to 
its interaction with a moving edge dislocation with a [1120] Burgers 
vector, The plane of the micrograph i is parallel to the basal plane of the 


zine platelet. 


Joa BART LEM andilevy, MIMGhEL Phil, Mag.Ser..8, Vols6;,El225: 
Fig. 1 


Interactions between dislocations with Burgers vectors at 120°. (x 2100.) 


Fig. 4 


a) oa: ft heroes Pa , ee i Le 
A group of interactions between dislocations with different Bur 


gers vectors. 
( x 1580.) 


G. K. WILLIAMSON and C. BAKER Phil. Mag. Ser. 8, Vol. 6, Pl. 26. 


Rige 


Dark field view of fig. 1 taken with 1120 reflection. 
Loop A has disappeared from contrast. 


G. K. WILLIAMSON and C. BAKER Phil. Mag. Ser. 8, Vol. 6, Pl. 27. 


Fig. 3 


Dark field view of fig. 1 taken with another 1120 reflection. 
Loops B and C and one side of loop F are out of contrast. 


Fig. 5 


Interstitial loops in irradiated graphite with the specimen tilted about 30°. 


G. K. WILLIAMSON and C, BAKER Phil. Mag. Ser. 8, Vol, 6. Pl. 28. 


(a) (0) 


(a) and (b). Dark field view of interstitial loop at slightly different reflecting 
conditions showing zero contrast along the diameter parallel to the 
reflecting plane (dotted line). 


A. H. MELEKA et al. Phil. Mag. Ser. 8, Vol. 6, Pl. 29. 
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The occurrence of cross-slip in fatigued copper. (a) 3 microns removed by 
electropolishing. (b) A further 9 microns removed. (x 340.) 


A. H. MELEKA et al. Phil. Mag. Ser. 8, Vol. 6, Pl. 30. 


Fig, 2 


As 


(0) 


The correspondence between extrusions: and cross-slip. (a) As fatigued. 
(6) 6 microns removed by electropolishing. (x 1350.) 


